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ABSTRACT: The liquefaction model presented based solely observed pore 
pressure responses sands and was constructed employing 
assumptions. Pore pressure generation sands under uniform well irregular 
loadings can simulated this model. The model requires only limited number 
laboratory tests and few model parameters. All model parameters can determined 
from results conventional undrained cyclic tests employing uniform stress cycles. 
Assumptions and performance this liquefaction model were examined through case 
studies. This liquefaction model predicted with good accuracy the pore pressure 
response for both uniform and earthquake-type irregular stress loadings. The model can 
incorporated into nonlinear site response methods and soil-structures response 
methods, and may used simple alternative other equally effective, but 
sophisticated liquefaction models currently available. 
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ABSTRACT: Ground movements associated with new insitu earth reinforcement 
lateral support system are analyzed using two-dimensional plane strain finite element 
technique specially developed for reinforced soil 
problems. Incremental construction and nonlinear soil characterization are included. 
The calculated ground movements are represented terms angular and lateral 
distortions. These were selected indications potential building damage ajacent 
open excavation. Typical design charts angular and lateral distortions relation 
the depth excavation and reinforcemeut length are presented. The system 
described this paper unique and many instances offers better choice over the 
conventional temporary earth support systems currently being used. The working stress 
analysis when coupled with the limit analysis calculation provide not only adequate 
design the system, but also information concerning ground movements critical 
building damage adjacent open cuts. 
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ABSTRACT: Two case studies lateral earth support system are presented. This 
system uses the insitu earth reinforcement technique, which different from the 
conventional systems that serve retain soils behind vertical cut, strengthen the 
native soil. The field results were compared with predictions obtained from finite 
element method analysis. general, the predictions agree well with the field 
measurement demonstrating that the analytical procedures developed can correctly 
predict the field behavior. Furthermore, has shown that ground movements adjacent 
open excavation can determined and controlled properly designed system 
installed. Hence, the insitu earth reinforcement lateral support system has shown 
offer viable alternative for providing lateral support deep excavation. 


REFERENCE: Shen, (Prof., Dept. Civ. Engrg., Univ. California, Davis, 
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ABSTRACT: description provided the design steel sheet pile cellular 
cofferdam required for construction graving drydock. The cofferdam was 
constructed approximately 550 (168 offshore within the Hood Canal 
Washington State and was required retain (24 water after basin 
unwatering. Interlock tensions associated with the required cell diameter this deep 
water location resulted one the first United States cofferdam applications high 
strength steel sheet piling and extruded wye connections. Much the cofferdam was 
designed permanent provide in-service laydown areas adjacent the 
completed drylock. Items discussed include the need for two pumped dewatering 
systems, the need for vibratory probe compaction the cell and are fill, site soil 
conditions, dredging and steel sheet pile corrosion protection. 


REFERENCE: Sorota, Max (Vice. Pres., Fay, Spofford Thorndike, Inc., Boston, 
Mass), and Kinner, Edward B., “Cellular Cofferdam for Trident Drydock: Design,” 
Journal the Geotechnical Engineering Division, ASCE, Vol. 107, No. GT12, Proc. 
Paper 16758, December, 1981, pp. 1643-1655 
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KEY WORDS: Cellular cofferdams; Compaction (soils); Dewatering; 
Dredging; Drydocks; Instrumentation; Steel sheet piles 


ABSTRACT: The construction and performance permanent steel sheet pile 
cellular cofferdam designed accommodate maximum water depth (24 
described. The cofferdam among the first constructed the United States 
which, because expected interlock tensions, required the use high strength steel 
sheet piles and extruded wye connections. Vibratory probe compaction the cell and 
connecting are fill was performed prohibit earthquake liquefaction, minimize 
settlements and obtain necessary weight for cofferdam stability. pumped cell fill 


dewatering system was used limit interlock tensions during and following cofferdam 
basin unwatering. 


REFERENCE: Sorota, Max (Sr. Vice Pres., Fay, Spofford Thorndike, Inc., 
Boston, Mass.), Kinner, Edward B., and Haley, Mark X., “Cellular Cofferdam for 
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ABSTRACT: Small-scale model footing tests were carried out two frictional soils 
having identical strength characteristics and different deformation characteristics. The 
results emphasized the importance soil compressibility. Test results are compared 
with existing analyical methods dealing with the bearing capacity compressible soils. 
The suitability these methods examined light the test results. 


REFERENCE: Ismael, Nabil (Assoc. Prof., Coll. Engrg. and Petroleum, Kuwait 
Univ., Kuwait), and Vesic, Aleksandar S., and Bearing Capacity,” 
Journal the Geotechnical Engineering Division, ASCE, Vol. 107, No. GT12, Proc. 
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16764 PROBABILISTIC SOIL EXPLORATION: CASE HISTORY 


KEY WORDS: Clays; Clay structure; Shear strength (soils); Soft soils; Soil 
investigations; Soil stratification; Two layer soil system 


ABSTRACT: Probability concepts are applied the interpretation soil exploration 
data site where failure weak layer considered probable. The subsoil 
modeled two-class material with soft clay layer included within stiff clay. 
hypothetical case history constructed illustrate the interpretation data obtained 
different stages the soil exploration program. The analysis the soil exploration 
program considers detection and recognition the soft material and inference that the 
soft material exists unexplored locations, given that has been detected explored 
locations. The judgment practicing engineers regarding the site conditions are 
presented subjective probabilities. 


REFERENCE: Wu, Tien (Prof. Civ. Engrg., Ohio State Univ., Columbus, Ohio. 
43210), and Wong, Kinfun, “Probabilistic Soil Exploration: Case Journal 
the Geotechnical Engineering Division, ASCE, Vol. 107, No. GT12, Proc. Paper 
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16735 LATERAL PILE RESPONSE DURING EARTHQUAKES 


KEY WORDS: Earthquake loads; Earthquakes; Liquefaction; Nonlinear 
systems; Offshore structures; Pile structures; Pore water pressure; Sands 


presented. Primary features the proposed method include consistent approach 
determine seismic p-y relationships and liquefaction model that suited pore 
pressure evaluation for irregular loadings. p-y 
relationships are determined from nonlinear stress-strain relations soils, and include 
soil nonlinearity well pore pressure buildup effects around pile. The 
performance and validity the method were evaluated through analysis shaking 
table test soil-pile-structure systems. Using the proposed method, seismic 
response characteristics idealized offshore pile-supported structure sand site 
were analyzed. Results indicated that liquefaction well pore pressure bulidup 
around pile can have large impact the response pile-supported structures. 


REFERENCE: Kagawa, Takaaki (Engrg. Consultant, McClelland 
Houston, Tex.), and Kraft, Leland M., Jr., “Lateral Pile Response During 
Earthquakes,” Journal the Geotechnical Engineering Division, ASCE, Vol. 107, No. 
GT12, Proc. Paper 16735, December, 1981, pp. 1713-1731 
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MODELING THE LIQUEFACTION PROCESS 
Takaaki and Leland Kraft, Members, ASCE 


INTRODUCTION 


Liquefaction cohesionless soils has contributed the failure soil structure 

systems during earthquakes. The mechanisms liquefaction have, recent 
years, been studied extensively, and many liquefaction models have been 
proposed. 

Available liquefaction models are based either: (1) Experimentally observed 
undrained stress paths during pore pressure buildup (5,6); (2) correlation between 
pore pressure response and volume change tendency dry soils (2,9); (3) the 
formulation pore pressure response directly from observed data (3,12,13,14); 
(4) plasticity theory which the plastic volume change related pore 
pressure build-up (10,18); (5) treatment the soil two-phase medium 
(1,8). Except for the models that formulate pore pressure response directly 
from observed data, all the other models explicitly implicitly include the 
correspondence between the volume change tendency soil and the pore pressure 
increase. Use volume change-pore pressure relationship sometimes, however, 
complicates the model and requires relatively large number model parameters 
that must obtained from several types laboratory tests. the other 
hand, the available models that formulate the pore pressure response directly 
from observed data undrained tests require less effort determine model 
parameters. The capabilities most models have been confirmed primarily 
for uniform loading cases. 

pore pressure generation model, based solely observed pore pressure 
responses, presented this paper. The model conceptually simple and 
requires only limited number laboratory tests and few model parameters. 
All model parameters can determined from conventional stress-controlled, 
undrained tests employing uniform stress cycles. The model also works well 
for earthquake-type stress variations. 

important feature the model that can used with any nonlinear 
dynamic response method that accounts for the time-dependent stiffness degrading 
effect due the pore pressure build-up the soil, and provides convenient 

Consultant, McClelland Engrs., Inc., Houston, Tex. 

McClelland Engrs., Inc., Houston, Tex. 

Note.—Discussion open until May 1982. extend the closing date one month, 
written request must filed with the Manager Technical and Professional Publications, 
ASCE. Manuscript was submitted for review for possible publication June 1980. 
This paper part the Journal the Geotechnical Engineering Division, Proceedings 
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alternative sophisticated liquefaction models evaluate the effects pore 
pressure generation within soil media. Due its simplicity, the model 
numerically efficient. 


Background.—During drained cyclic loading, loose soils densify and dense 
soils dilate. Similar phenomena also are observed during undrained cyclic loading. 
This volume change tendency key influence the liquefaction cohesionless 
soils. 

When element relatively loose soil subjected cyclic stress reversals 
undrained simple shear condition, pore pressure develops the soil. 
early stage loading when the pore pressure ratio low, the effective stress 
system the soil nearly identical that before loading. Therefore, the 
pore pressure increment during undrained stress cycle may adequately 
predicted using the potential for drained volume change corresponding 
the initial stress system. the cyclic straining proceeds generate large pore 
pressures the soil, the effective confining pressure reduces and the soil may 
exhibit tendency dilate. The pore pressure response this stage can 
longer predicted from the volume change tendency the soil under the 
initial stress system. 

Pore pressure generation can related the drained volume change charac- 
teristics the soil, which are affected the cyclic shear strain and the change 
effective confining pressure based the effective stress point view (2,9). 
This approach, however, requires experimental data cyclic stress-strain 
behavior for the soil, relationship between the cyclic strain and the volume 
change behavior dry soil, and relationship between the volumetric expansion 
and the effective stress reduction for the soil. 

simple alternative the effective stress method may obtained 
empirically relating the observed pore pressure responses the drained volume 
change the soil (11,16) assuming that the volume change tendency 
the soil during undrained cyclic loading can described the volume change 
the dry soil under the initial stress system. The stress ratio chosen 
key parameter these methods. Due the simplifying assumptions involved 
these methods, predicted pore pressure responses may not adequate for 
any type loading and soil conditions. major advantage this approach, 
however, simplicity. 

this study, the writers constructed liquefaction model based the observed 
pore pressure responses laboratory undrained cyclic tests employing uniform 
stress cycles. The liquefaction process viewed result two distinct 
phenomena: (1) early stages loading, pore pressure generation dominated 
the tendency the soil decrease volume under the initial effective 
system; and (2) the pore pressure generation capability the soil also 
affected the change effective confining pressure during cyclic loading. 
The former characterized using the drained volume decrease behavior 
sands, and the latter represented empirical ‘‘pore pressure 
Although the key model parameter the stress ratio, the ‘‘pore pressure 


implicitly accounts for the strain dependent nature the liquefaction process, 
analysed below. 
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Assumptions.—Three fundamental assumptions were made construct this 
liquefaction model: 


Pore pressure generation due volume change characteristics, and this 
influenced by; (1) The magnitude applied shear stress for given confining 
pressure and relative density the soil; (2) the stress history (repeated loading); 
and (3) the effective increase the applied shear stress relative the shear 
strength, effective stresses decrease with increase pore pressure. 

The pore pressure generation sand depends the magnitude applied 
shear stress (3,11,13,14,16). Previous studies also have shown that pore pressure 
generation sand related the volume decrease tendency the soil 
under repeated loading. Under drained conditions, volume during uniform 
stress cycle decreases rather quickly with increasing number cycles. 
this volume decrease sand were directly converted the pore pressure 


SHEAR STRESS, 


TIME, 


FIG. Pressure Computation Scheme 


increase the sand undrained condition, the amount pore pressure 
build-up stress cycle would decrease with increase the number 
cycles. Large pore pressures, however, are generally developed high pore 
pressure ratios near liquefaction. This phenomenon may consequence 
increase cyclic shear strain amplitude near liquefaction, can attributed 
the third influence previously meationed. 

Soil exhibits identical pore pressure generation both positive and 
negative shear stress applications only when the stress histories are the same 
both directions and the pore pressures before the stress application are the 
same. The soil remembers the stress histories independently for the positive 
and negative shear stresses. 

Pore pressure assumed develop only the loading and reloading 
process. earthquake-type irregular stress-time history shown Fig. 
The portions which irrecoverable pore pressure develops are represented 
solid lines. For example, the stress changes from 0-1, 6-7, 8-9, and 
12-13 correspond loading reloading the positive shear stress. Thus, 
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pore pressure develops during these stress changes. the other hand, the 
stress changes from 1-2, 5-6, 7-8, 9-10, and 13-14 are unloading processes 
the positive shear stress. Thus, pore pressure, according the assumption, 
develops during these stress changes. Similarly, pore pressure develops also 
during loading and reloading processes the negative shear stress (from 2-3, 
10-11, and 16-17). 


Model Description.— When for pore pressure build-up, 
exists, can compute the pore pressure increase due shear stress 


Ir2/8,| t 


which pore pressure increase; initial effective overburden stress; 


x 
< 
x= 
a 


SHEAR STRESS, 


(b) IRREGULAR LOADING 


FIG. 2.—Pore Pressure Evaluation and Memory Function, 


and shear stress. considered uniquely determined for given soil 
density and confining pressure following the first assumption, and expressed 


which and model parameters; n(N) that depends 


Cc, 
STRESS FUNCTION 
RANGE 
TIME, 
UNIFORM LOADING 
RANGE 
0-a 
c-d 
h j-k 
+ 
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DEVIATOR STRESS, 


PORE PRESSURE, 


NUMBER OF CYCLES 


(AFTER ISHIHARA AND YASUDA, 1972) 


FIG. 3.—Typical Pore Pressure Response 


Nth CYCLE 


LOG 


FOR GIVEN 
SOIL DENSITY 


LOG 


FIG. 4.—Determination Parameter 
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the number the times soil element sheared; and function 
representing the change the pore pressure generation capability soil 
element caused increase the pore pressure ratio, The function 
n(N) reflects the tendency for volume decrease, and the function, 


TABLE 1.—Summary Model Parameters Case Studies 


gram 
force per 

square 

centi- 

Sand sample meters 


Niigata sand 
Fuji River sand 


MODEL 
PARAMETERS 


NIIGATA SAND 


1.5 


~ 


=0.2 


NUMBER CYCLES, 


FIG. 5.—Model Parameter Determination Procedure 


0.20 3.14 2.06 187.5 
0.27 5.85 1.94 76.2 
0.155 6.98 1.72 402.8 
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accounts for the effective increase the applied shear stress result 
decrease confining pressure during the liquefaction process. Alternatively, 


NIIGATA SAND 
45% 


~0.20 


PORE PRESSURE FUNCTION, 


0.5 
PORE PRESSURE RATIO, 


FIG. 6.—Pore Pressure Function for Niigata Sand 


0.4 


TEST DATA (ISHIHARA AND 
YASUDA, 1972) 


PREDICTED 


NIIGATA SAND 
45% 
1.5 
1.0 


STRESS RATIO, 


NUMBER CYCLES LIQUEFACTION 


FIG. and Predicted Stress Ratio Liquefaction for Niigata Sand 


may considered empirical function that represents the effects 
increase cyclic shear strain near liquefaction and the volume expansion 
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tendency the soil due the reduction effective confining pressure 
the pore pressure generation capability sand. 

The memory function, n(N), its definition, depends the number 
times soil element sheared some stress ratio, Fig. shows variations 
the memory function, n(N), when stress-time histories are applied soil 
element. According the third assumption, soil element sheared only 


1.0 
04/25, 0.20 


SAND 
0, =45% 
1.5 
© TEST DATA (ISHIHARA AND 
1972) 


@ PREDICTED 


PORE PRESSURE RATIO, 


NUMBER CYCLES, 


FIG. and Predicted Pore Pressure Response for Niigata Sand and 
Uniform Stress Cycles 


5 6 7 
TIME, SECONDS 


= 
= 
«< 


WHGATA SAND (0, 45%) 
} ISWIMARA AND 
YASUDA (1972)" 

PREDICTED 

os; = 1.5 


PORE PRESSURE. u. 


TIME, SECONDS 


FIG. and Predicted Pore Pressure for Niigata Sand and Irregular Loading, 
Niigata Earthquake W-E Component 


during the loading and reloading processes. For uniform loading case, Fig. 
2(a), the memory function during the stress changes (0-1 and 2-3) 
the soil element sheared for the first time. The memory function for the 
stress changes (4-5 and 6-7), however, the soil element sheared 
for the second time these stress ranges. Variation for irregular 
stress history also will defined similar manner (Fig. 2(b)). the other 
hand, the function will considered independent the stress history 
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and uniquely related the pore pressure ratio, for given soil density 
and confining pressure. 

Determination Model Parameters.—Most pore pressure data have been 
obtained from undrained cyclic tests employing uniform stress cycles. Therefore, 
convenient for construct model the parameters which can 
determined from such tests. For uniform stress conditions, Eq. can integrated 
using Eq. 


which 1); and the applied shear stress amplitude. 

Fig. shows typical pore pressure response undrained cyclic triaxial 
test employing uniform stress cycles. The pore pressure increment each stress 
cycle can determined from results such those Fig. When cyclic 


= 
a 
a 
s 


SECONDS 


—— CM-TEST)| (SHINARA AND | MNGATA SAND (0, 
vasuba (1972) 


prroicteo 


PORE PRESSURE, u. 


TIME, SECONDS 


FIG. 10.—Measured and Predicted Pore Pressures for Niigata Sand and Irregular 
Loading, Niigata Earthquake N-S Component 


undrained tests are performed samples with the same density and confining 
pressure but for several different stress ratios, the parameter can obtained 
function the number stress cycles. method determine the parameter 
for practical values typically ranges from 0.0-3.0, based data 
our files and results previous studies (3,11,13,14,16,17). 

With the value thus determined, and appropriate shear stress amplitude, 
(Fig. 3), Eq. can simplified 


Cc, Co, ( 


Functions and can determined from pore pressure response 
such that Fig. determine these functions using test results with 
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several stress ratios, must employ the results with stress ratio that 
representative the field loading condition. The pore pressure response 
n(N), recall: (1) That the function n(N) dominates for small values, 
this function represents the pore pressure increase due the potential for 
volume decrease sand under the initial effective stress system; and (2) 
the drained volume decrease sand can approximated hyperbolic 
function for practical values (2,9,11,16). Thus, the initial straight line 
portion the curve Fig. must represent the pore pressure increase due 
the volume decrease potential the sand under the initial effective stress 
system that given 


which parameters and can determined from Fig. Therefore, the 
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FIG. 11.—Summary Pore Pressure Function Case Studies 
pore pressure increase during the cycle, due the volume decrease potential 


The observed pore pressure increase the Nth cycle Fig. excess 


(Pore Pressure Increase during the Nth cycle) n(N) 


The function, obtained for the pore pressure response Fig. 
normalized n(1) and shown Fig. 
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Assumptions and performance the liquefaction model were examined 
analyzing pore pressure responses obtained laboratory tests for both uniform 
and earthquake-type loading conditions. 
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FIG. 12.—Measured and Predicted Stress Ratio Liquefaction for Fuji River Sand 
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FIG. 13.—Measured and Predicted Pore Pressure Response for Fuji River Sand with 


Case Study and Yasuda (7) performed undrained triaxial tests 
Niigata sand employing input both uniform and earthquake-type irregular 
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stress histories. All the parameters and the functions, n(N) and for 
the liquefaction model except were determined using their test data for uniform 
stress loading (Table 1). Available data were not sufficient determine the 
parameter, Therefore, was assumed, and value 2.0 was used the 
analysis. The function, was determined, shown Fig. 

Fig. shows the relationship obtained between the stress ratio and the number 
cycles liquefaction for undrained triaxial tests Ishihara and Yasuda 
(7) and the analytic model. The agreement between the computed and observed 
results excellent. The pore pressure generation process obtained the model, 
however, deviated slightly from that the test data, due the nonlinear 
relationship between and (Fig. 8). 

further examine the performance the model, the writers predicted the 
pore pressure response the Niigata sand under earthquake-type, irregular 
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FIG. 14.—Measured and Predicted Pore Pressure Response for Fuji River Sand with 


loadings. Fig. shows the results for the stress-time history that simulates 
the W-E component the Niigata earthquake. Ishihara and Yasuda (7) reported 
results from tests’’ which the maximum shear stress developed 
the compression side the triaxial test, and tests’’ which the maximum 
shear stress developed the extension side. The pore pressure response from 
the model excellent agreement with the test results. similar degree 
agreement was also achieved for the stress-time history simulating the N-S 
component the Niigata earthquake (Fig. 10). 

Case Study torsional shear testing device, Ishihara, Iwamoto, 
Yasuda, and Takatsu (4) studied the effect the earth pressure coefficient 
rest, the liquefaction potential Fuji River sand. Their tests included 
both uniform and irregular stress loadings. 

The parameters and the functions n(N) and were determined using 
the test results the same way Case Study but for equal 0.5 


GT12 


SHEAR STRESS, T. PSI 


PORE PRESSURE 


and 1.5 (Table and Fig. 11). The value was assumed 2.0. 

The relationships between the stress ratio and the number cycles 
liquefaction from the test data and the model, are summarized Fig. 12. Figs. 
and compare computed and measured pore pressure histories. Again, 
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FIG. 15.—Measured and Predicted Stress Ratio Liquefaction for Ottawa Sand 
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FIG. 16.—Measured and Predicted Pore Pressure Response for Ottawa Sand and 
Irregular Loading 


the model satisfactorily predicted the pore pressure response for both uniform 
and irregular loading conditions. 


Other Case Studies.—Beside the case studies described above, performance 
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the model was examined using the test results reported Ishibashi, Sherif, 
and Tsuchiya (3), and Sherif, Ishibashi, and Tsuchiya (13). their data were 
not sufficient determine all the model parameters and functions, assumed 
appropriate values. The results for uniform stress conditions and earthquake-type 
irregular loading are shown Figs. and 16. could reproduce their test 
results with excellent degree accuracy, although the level simulation 
this case lower than those the two aforementioned case studies. 

Performance the model also has been studied using data our files 
undrained cyclic triaxial tests several silty sands. Although the comparisons 
these cases were limited uniform stress conditions, the agreement between 
computed and observed results was quite favorable. 


liquefaction process was modeled employing some simplifying assumptions. 
The liquefaction model presented this paper predicted, with good accuracy, 
the pore pressure response for both uniform and irregular stress loadings. 
the case studies presented, the model parameters were determined using results 
from single, undrained test employing uniform stress cycles. The parameter 
was tactically assumed. Results from undrained cyclic tests for least three 
stress ratios, however, are needed define the parameter when applying 
the model practical problems. course, further studies are needed define 
the system parameters more detail and relate the system parameters 
the status soil, especially the relative density soil and initial confining 
pressure. 

The major advantage the liquefaction model presented this paper 
its simplicity. The model requires relatively small number parameters that 
can determined from conventional undrained cyclic tests employing uniform 
stress cycles. The model requires less laboratory testing than most liquefaction 
models currently available that can incorporated into effective stress seismic 
response methods. The pore pressure evaluation with this model computation- 
ally efficient. Thus, the model may used simple alternative other 
equally effective, but more sophisticated liquefaction models currently available. 
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The following symbols are used this paper: 


model parameter; 
pore pressure function; 


n(N) memory function; 
number stress cycles; 
excess pore-water pressure; 
model parameter; 
pore pressure increment; 
initial effective overburden stress; 
initial effective ambient stress; 
shear stress; and 
potential function for pore pressure generation. 
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GROUND MOVEMENT ANALYSIS 
EARTH SYSTEM 


and Herrman,’ ASCE 


INTRODUCTION 


recent years, relatively new situ earth reinforcement lateral support 
system for deep excavation has been introduced. Unlike the conventional systems 
that serve retain soil behind vertical cut, this procedure based the 
concept soil reinforcement, the native soil adjacent the excavation 
strengthened that can stand unsupported depths which would normally 
require the installation sheet pilings soldier piles and bracings. Briefly, 
this system composed array reinforcing members that are grouted 
into the soil mass, wire-mesh reinforced shotcrete panel facing, and rows 
re-bars which form horizontal wales each reinforcement level. The various 
components this system are shown Fig. 

Excavation begins from the ground level. After each layer excavation, 
reinforcement immediately applied both the native soil and the exposed 
cut. This system provides unique means temporary earth support with 
the advantages that it: (1) Requires pile driving; (2) prevents any loosening 
sloughing the soil; (3) provides obstruction-free site for subsequent 
foundation work; and (4) significantly reduces construction time. The system 
has far been successfully used for thousands square meters excavation, 


found elsewhere (10). 


Movements Open Excavation 


Ground movements associated with deep open excavation have often been 
observed the field. Unless permanent type often costly lateral earth 
support system (such slurry walls) used, ground movement adjacent 
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open excavation almost unavoidable. the present time, however, there 
rational method capable accurately predicting the extent well 
the magnitude both the horizontal and vertical movements. This due largely 
the fact that construction practice temporary lateral earth support systems 
varies from place place, and, frequently, its execution entirely the 
discretion the contractor. The lack standardized construction procedures 
together with other local influences make difficult, not impossible, quantify 
ground movements associated with deep excavation. Consequently, only limited 
field information available pertaining this subject. For instance, Peck (9) 


Existing Grade 


2'-0 ) Adjust Elev. and Length of 
to Top Reinf. to Suit Location 
3'-0 J of Utilities, etc. 


First Appl. 2” 4 Second Appi. 2” 


With Grout a 


#68 Re-bar (Typ) 


6" x 6" x 1/2” Plate 


4° & 8/8 Wire Mesh 


view A (Typ) 


FIG. 1.—Construction Detail—Earth Reinforcement Lateral Support System 
0.305 in. 25.4 mm) 


studied the field records many excavation sites using soldier piles sheet 
piles that were braced and classified the field performance open excavation 
systems into three categories Fig. Zone which gives the best performance, 
refers well-braced systems with slurry walls substantial amount berms 
left permanently place. Zone II, describing intermediate performance, indicates 
systems with temporary berms and raker support. Zone III includes systems 
with ground loss from caisson construction insufficient support. Similar 
information gathered Chicago and Washington D.C. areas were reported 
al. (8). These records have clearly indicated that the amount and 
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distribution ground movements depends heavily upon the soil type, the method 
and sequential procedure construction, and the depth excavation. Since 
control these factors extremely difficult using the current excavation 
methodology, unlikely that accurate prediction ground movements can 
made possible. Ground movements, however, could bear grave consequences 
the stability adjacent buildings, the operation utilities, and the use 
nearby streets. adopting particular temporary lateral earth support system, 
one should ideally able judge estimate the magnitude and distribution 
ground movements appraise the potential for damage and evaluate 
the need for proper design protective measures both. Unless the procedures 
construction practice can controlled and standardized, the question 
ground movement prediction will remain unresolved. 

the other hand, the construction procedures lateral earth support 
system can standardized and systematically carried out the field and 
physical model can devised properly describe the system behavior and 


DISTANCE FROM EXCAVATION DEPTH EXCAVATION 


SETTLEMENT DEPTH EXCAVATION (%) 


FIG. 2.—Field Performance Open Excavation System (Ref. 


its construction sequence, possible that analytical procedure may 
formulated predict the ground movements more accurately. The system 
described this paper has unique and well-defined construction procedure 
such that can expected achieve better and more uniform control 
field practice than any currently employed temporary system. For this reason, 
appears that more rigorous analytical treatment can applied order 
determine the magnitude and distribution ground movements associated 
with deep excavation. Such analyses and the interpretation the analytical 
results form the remaining parts this paper. 


Computer Program.—Analyses were performed using the finite element 
computer program for general two-dimensional soil and soil reinforcement 
problems developed Herrmann (5). The program assumes that plane strain 
conditions exist. Because the finite width and the horizontal spacing 
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reinforcement the system studied, did not exactly satisfy the plane strain 
criterion. However, could reasonably assumed that deviation from plane 
strain conditions was not substantial for most practical purposes (6). Several 
elements were incorporated into the analyses simulate the special features 
the system properly; these include the interface bond layer, the sequential 
procedure construction, and the composite reinforced soil model. more 
detailed description given the following sections. 

Bond Layer.—Slippage between the reinforcing members and the surrounding 
soil was accounted for the analyses (5). The maximum bond stress was assumed 


INITIALIZE THE 
IN-SITU SOIL CONDITION 
EXCAVATE FIRST LAYER SOIL 


INSTALL FIRST LAYER REINFORCEMENTS 


AND PANEL FACING. EXCAVATE SECOND 
LAYER SOIL. 


EXCAVATE LAYER SOIL 


INSTALL LAYER REINFORCEMENTS AND 
PANEL FACING. EXCAVATE SOIL 


EXCAVATE LAST LAYER SOIL 


INSTALL LAST LAYER REINFORCEMENTS 
AND PANEL FACING 


FIG. 3.—Construction Steps Used Finite Element Analysis 


was assumed that adhesion was fully maintained even though slippage might 
have taken place. yielding re-bars relative the grout was considered. 
Construction Sequence.—The analyses were carried out sequence 
incremental steps with iterations conducted each increment. The incremental 
procedure simulated the sequential excavation steps and the installation the 
reinforcing members. The number increments required was determined 
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such way that the excavation and the stress history the soil were adequately 
represented. Generally, the greater the number increments, the better the 
accuracy the analysis. flow chart the computational steps the analysis 

Composite Model.—A composite model was used describe the reinforced 
soil system. The composite representation adopted for this analysis was different 
from the conventional composite approach expressing simply the orthotropic 
composite material properties functions the properties each the 
constituent materials and their geometric arrangement. The procedure used 
this study accounted for the composite nature the reinforced soil system 
combining the element stiffness matrices the constituent materials (i.e. 
soil, reinforcing members, and the bond behavior) directly produce the 
composite element matrices. The approximate solution was then selected 
minimizing the incremental potential energy the system. Detailed description 


TABLE 1.—Properties Soils Used Analysis 


Soil number 
(2) 


Soil number 


0.3 (assumed) 
121.3 pef 

0.013 
0.045 
0.070 


Note: psi 6.9 and pef 157 N/m’. 


the development can found Ref. The soil model employed the 
analysis was the nonlinear, and the inelastic model was proposed Duncan 
and Chang (4). This characterization assumed that the soil behavior can always 
described the instantaneous values modulus, i.e., the soil does not 
experience stress strain induced anisotropy. hyperbola was assumed 
represent the stress-strain relationship for primary loading and straight line 
for unloading and reloading. Duncan and Chang also proposed similar expression 
for the Poisson’s ratio. They, however, failed include the effects unloading- 
reloading Poisson’s ratio. was found during this study that the expression 
for Poisson’s ratio did not represent the proper behavior soil during excavation, 
particularly the prediction ground surface movements. constant Poisson’s 
value 0.3 was, therefore, used the analyses. 

Analysis.—Based previous parametric study (1), was concluded that 
the construction sequence; the depth excavation, the soil type; the length 


100 280 
K,, 300 840 
0.84 0.6 
27° 18° 
7.5 psi 12.64 psi 
0.77 0.93 
0.3 (assumed) 
124.6 pef 
0.002 
0.04 
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the size, and the spacing the reinforcing members were important 
parameters considered design. was also indicated that the reinforcing 
members were best oriented horizontal inclination between 15-25°. 
examine the effects these parameters ground movements further, the 
following analyses were performed. 

The depths excavation studied were range approx 0.5-1.5 times the 
length the reinforcing members 15-40 (4.5-12 length. angle 


FIG. 4.—Typical Finite Element Grid 0.305 
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in. 25.4 mm) 


was used for the reinforcement grid pattern. The diameters the reinforcing 
members were 3.33 in. (85 mm) and in. (127 mm), corresponding 0.24 
and 0.54% the total cross-sectional area the reinforced soil mass, respectively. 
the plane strain analysis, equivalent reinforcing members that had tensile 
strength equal that steel and contact area per unit width 
wall were used, where the diameter the grout hole; the horizontal 
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spacing reinforcement; and the length reinforcement. Two different 
soils, sandy-silt (soil #1) and sandy clay (soil #2) were chosen. Their physical 
properties and the hyperbolic parameters were obtained from Wong and Duncan 

Fig. shows typical finite element mesh and its boundaries for system 
having 15-ft (4.5-m) reinforcement length and maximum depth excavation 
(7.5 m). The shaded portion the mesh represents the reinforced 
soil forming the situ lateral earth support system. the analyses, the reinforced 
soil mass was modeled composite elements, and the facing was represented 
one-dimensional elements. Six increments were required simulate the 
excavation depth (7.5 and the installation reinforcement; 
one more increment was added generate the initial geostatic stresses the 
soil mass. 

The successive displacement profiles the ground surface adjacent the 
excavation for soil (for the case in. (127 mm) and 


DISTANCE FROM EXCAVATION DEPTH EXCAVATION 


in. 25.4 mm) 


(4.5 m)) are shown Fig. can seen that small amount rebound 
the ground surface near the edge occurs during the first (1.5 
excavation. The maximum amount rebound approx 0.3 in. (7.5 mm) which 
corresponds approx 0.5% the depth excavation. The deformed shape 
the facing approximately linear. However, excavation proceeds, the 
pattern ground movement generally concave downward with the maximum 
movements occuring along the edge the ground surface. The maximum 
horizontal displacement increases rapidly with increasing depth excavation. 
For instance, maximum horizontal displacement 0.5 in. (12.7 mm) (0.75% 
(7.5 excavation, the displacement 5.1 in. (13 mm) (1.7% the 
depth excavation). This indicates that the effectiveness reinforcement 
decreases with decreasing ratio reinforcement length the depth excavation, 
L/H. Furthermore, results ground surface settlement for soil are plotted 
Fig. Curves and represent the settlement profiles with reinforcement 
lengths 40, 30, and (12, and 7.5 m), respectively. The depth 
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excavation (12 m), and the diameter the reinforcing members 
in. (127 mm) all cases. The ratios L/H for curves and are 
1.0, 0.77, and 0.63, respectively. These results show again that the effectiveness 
reinforcement decreases with decreasing ratio Comparing Figs. 
and suggested Peck, can seen that, properly designed, the 
anticipated ground movements associated with this system appear well 
within the movements normally identified with currently used temporary lateral 
earth support systems construction practice. 


Discussion 


Conventionally, the safety evaluation earth retaining structure based 
stability calculation performed against failure. This limit analysis approach 
design methodology does not provide any information concerning the movement 
displacement the ground adjacent open excavation. result, 
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different values safety factors are imposed indirectly limit the movement 
the structure maintaining the shear stresses everywhere within the 
soil-structure system well below its shear strength. Unfortunately, the margin 
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FIG. 9.—Angular Distortion Soil No. 0.305 in. 25.4 mm) 


safety against failure alone not sufficient judging the damage potential 
surrounding buildings other installations. additional design criterion 
needed, therefore, relate, the working stress level, the ground movements 
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excavation depth and reinforcement. These movements can then translated 
damage potential surrounding structures that the appropriateness 
the design can properly evaluated. 
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FIG. 10.—Lateral Distortion Soil No. 0.305 in. 25.4 mm) 
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FIG. 11.—Lateral Distortion Soil No. 0.305 in. 25.4 mm) 


Ground movements expressed angular and lateral distortions have been 
correlated building damages reported recently O’Rourke, Cording, and 
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Boscardin (8). The definitions angular and lateral distortions are shown 
Fig. The angular distortion, 5,, the differential settlement between two 
points divided the horizontal distance separating them. The term angular 
distortion has often been related architectural and structural damages 
buildings due settlement (3). However, buildings near open cuts are also 
affected lateral displacements. These movements can significant, and 
can cause substantial amount strain adjacent structures. The lateral 
distortion, 5,,, defined the difference lateral movement between two 
points divided the horizontal distance separating them. Correlations between 


damages suffered old brick buildings lateral distortion can found 
Ref. 
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FIG. 12.—Angular Distortions from the Excavation 0.305 in. 
25.4 mm) 


Using the data obtained from the finite element analysis, Figs. and are 
presented show the distribution angular distortions for various depths 
excavation and reinforcement sizes for soils and #2. The length the 
reinforcement (4.5 for all cases shown. The angular distortion between 
two surface node points first calculated and then plotted the midpoint 
between the two nodes. noted that the further the distance from the edge 
excavation, the smaller the magnitude angular distortion. The rate 
decrease angular distortion with respect the distance from the edge 
excavation substantial; only negligible amount distortion recorded 
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when that distance equals the depth excavation Fig. for soil #1. 
Soil experiences even smaller amounts angular distortion. The effect 
reinforcement size the ground movement relatively insignificant 
comparison the L/H ratio. This because the prime source resisting 
force this system derived from friction and adhesion between the reinforcing 
members and the soil. can shown that for the same amount reinforcement 
terms volume per unit width the increase frictional resistance per unit 
length reinforcement can achieved more readily increasing the contact 
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FIG. 13.—Angular Distortions from the Excavation 0.305 in. 
25.4 mm) 


area; i.e., instead increasing the size reinforcing members, more 
beneficial reduce their spacings the soil mass. 

The lateral distortions, which are measurements surface tensile strains 
for the same excavation system, are plotted Figs. and for soils 
and #2, respectively. For relatively shallow depths excavation, the maximum 
lateral distortion occurs near the edge, but the pattern changes the depth 
increases. The maximum lateral distortion, for instance, located distance 
approx 6.7 from the edge when the excavation reaches depth 
for soil and (7.5 for soil #2. This indicates that the 
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largest tensile strain the ground surface does not always develop along the 
edge excavation but rather some distance away from it. This distance, 
however, determined largely the soil type, the characteristics the 
reinforcement, and the depth excavation. The locations maximum lateral 
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Figure 14. Lateral Distortions ft. From the Excavation 


FIG. 14.—Lateral Distortions from the Excavation 0.305 in. 25.4 
mm) 


distortion predicted the analyses seem agree with field measurements 
braced wall system reported O’Rourke, al. (8). According their 
measurements open excavation system having five rows struts, the 
maximum lateral distortion occurred near the edge excavation when was 
deep; and approx (4.5 from the edge when the excavation 
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was (18 deep. The predicted maximum lateral distortion appears 
relate tension cracks often observed around open cuts. 

Useful design information concerning both the angular and lateral distortions 
the ground surface can derived from Figs. and for soil #1. The 
distribution angular distortion with respect the depth excavation and 
reinforcement length are shown Figs. and distances (1.5 
and away from the edge excavation, respectively. Likewise, 
the corresponding distributions lateral distortion are shown Figs. and 
15. For given set reinforcement size and spacing, the magnitude the 
angular distortion increases with the depth excavation, and decreases with 
the length the reinforcement, i.e., greater ratio. The variation lateral 
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FIG. 15.—Lateral Distortions from the Excavation 0.305 in. 
25.4 mm) 


distortion the ground surface not always monotonically decreasing function 
the case for angular distortion. decreases monotonically for relatively 
shallow depths excavation, but, excavation becomes deeper, the maximum 
distortion occurs some distance away from the edge. This indicates that 
the lateral distortion near the edge decreases the depth excavation increases 
shown the curves Fig. 14. This pattern, however, does not hold 
distances further away from the excavation Fig. 15. 

Having constructed the design charts and knowing the constraints terms 
permissible ground distortions possible building damage adjacent 
open cut, one may determine first, the proposed system suitable choice 
for particular job; or, second, the required reinforcement spacing and length 
for given depth excavation that ground distortion can limited within 
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its tolerance. This design methodology unique its approach, and applies 
only the situ lateral earth support system described above. may 
noted that the illustrations given here have not included any form surface 
loading (such building, fill, etc.) the analysis. However, provisions are 
available the computer program approximate such loading conditions. The 
tolerance ground distortion buildings and other structures subject 
great importance foundation engineering; however, its formulation beyond 
the scope this paper. 


Summary 


Ground movements associated with new situ earth reinforcement lateral 
support system were analyzed using two-dimensional plane strain finite element 
technique specially developed for reinforced soil and soil-structure interaction 
problems. Incremental construction and nonlinear soil characterization were 
included. The calculated ground movements were represented terms angular 
and lateral distortions. These were selected indications potential building 
damage adjacent open excavation. The maximum angular distortion 
the ground surface occurs the edge excavation, and its value decreases 
rapidly with respect the distance from the open cut. The maximum lateral 
distortion the ground surface, however, may occur some distance away 
from the edge; this distance depends the soil type, the characteristics 
reinforcement and the depth excavation. Typical design charts angular 
distortion and lateral distortion were developed from the finite element analysis 
results. Similar plots can established for different types soils, for different 
size and length reinforcement, and various distances from the edge 
excavation. 

believed that the system described above unique and attractive 
one, and, many instances, offers better choice over the conventional 
temporary earth support systems currently being used. The working stress 
analysis, coupled with the limit analysis calculation (2), provide not only 
adequate design the system but information concerning ground movements 
critical building damage adjacent open cuts. 
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FIELD MEASUREMENTS EARTH 
SYSTEM 


INTRODUCTION 


The lateral earth support system using the situ earth reinforcement technique 
has been described accompanying paper (8). For this system recognized 
and accepted the geotechnical and construction profession, necessary 
that documented field performance records are made available. This paper 
presents the field measurement records two sites—the Portland study and 
the Davis study—and their comparison with results obtained from the finite 
element method analysis. 

The first application this system the United States was the foundation 
excavation site for the expansion the Good Samaritan Hospital and Medical 
Center Portland, Oregon, the summer 1976. Since the system was relatively 
new and unfamiliar the local engineers and contractors, limited amount 
field instrumentation was installed monitor its performance and safety. 
This instrumentation was not meant, however, provide complete set 
data for field evaluation and analysis. The Davis study was carried out the 
summer 1979. This was field investigation with extensive instrumentation 
and data collection. The field work was part comprehensive study 
the situ earth reinforcement lateral support system funded the National 
Science Foundation. Descriptions and comparisons the field measurements 
together with analyses the two sites are reported herein. 
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General.—The construction site was located immediately west the existing 
hospital shown the plot plan Fig. The general street level elevation 
around the site was about 110 (33.6 m). basement with finished floor 
grade established elevation (21.4 m), extended under the entire new 
addition. The earth reinforcement lateral support system was installed along 
the north (Marshall Street), the west (NW 23rd Avenue), and the south (NW 
Lovejoy Street) boundaries depths varying from 35-45 (10.7-13.7 m). 
total 683 reinforcing members were installed and lateral surface area 
about 23,000 (2,140 was covered 4-in. (100-mm) thick wire-mesh 


EL. WORTHRUP ST 


WARSHALL ST. (TOBE VACATED) 
EL. 1093 ft 


bd LOVEJOY ST 
EL 


FIG. Excavation Good Samaritan Hospital, Portland, Oregon 
0.305 


reinforced shotcrete facing. The westerly boundary along 23rd Avenue 
was the longest, approx 250 (76.3 m), with depths excavation ranging 
from (12.8 the southern end approx (10.7 near the 
north. The deepest cut (45 ft) was along Marshall Street, where the boundary 
excavation was parallel and about away from four-story brick 
building (the Marshall Street Annex). Fig. shows general view the excavation 
site near its completion; the brick building and its proximity the edge 
excavation shown Fig. 

Subsurface the record nearby well log, was 
indicated that the general area the site was underlain 70-90 (21.4-27.5 
thick formation silty, sandy lacustrine soil. Below this was layer 
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gravel, cobble, and boulder formation approx (27.5 thick. The basaltic 
bedrock was found elevations between (—16.8-—18.3 m). 
subsurface soil investigation (4) for this site was carried out drilling four 


FIG. 3.—Marshall Street Annex 


boreholes depths about (29 m). typical borehole log shown 
Fig. the top was brownish fine sandy silt which was followed 
silty fine sand with randomly interbedded layers fine sand and silt 


FIG. View Site 
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depth about (10 m). Below that depth was layer predominantly 
fine sands with occasional thin layers silt. The transition zone between the 
lacustrine soil formation and the underlying gravel formation was observed 
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FIG. 5.—Horizontal Deflection Profiles From inclinometer Readings 0.305 


elevation (6.7 m). ground water table was encountered within the 
depth explored. 

Field Instrumentation and Measurements.—This being the first job using the 
situ earth reinforcement lateral support system for deep excavation the 
United States, limited field instrumentation program was initiated monitor 
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and safeguard the stability the excavation during construction. Among the 
measurements obtained, two sets information were selected for comparison 
with analytical results: the lateral ground movement profiles from borehole 
inclinometer readings and the surface tension cracks the ground floor slab 
the Marshall Street Annex. 

The casing-lined borehole for inclinometer readings was (13.7 deep 
located from and approximately halfway along the length the 
westerly boundary. The depth excavation that location was about 
(11.3 m). continuous record the ground movements was taken during and 
right after the completion excavation Fig. however, detailed 
information relating depth excavation movement available. The maximum 
lateral movement the ground surface was 1.3 in. (33 mm), which was only 
about 0.3% the maximum depth excavation. 

The Marshall Street Annex, four-story brick building, being close 
the excavation, experienced severe tensile cracks the ground floor concrete 
slab. This building had additional half-story basement below the street level. 


=1 to 7 Cracks in Floorsiab 
8 Crack in Wall inside 
=9, 


° 
@ =9,10 Cracks in Window Sill First Floor 
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EXCAVATION 


FIG. 6.—Crack Survey—Marshall Street Annex 


The ‘‘surcharge top the northerly lateral support system was, 
therefore, quite substantial. This was the most critical area the excavation, 
insofar safety was concerned. The depth excavation along the northerly 
boundary was approx (11.3 m); however, large footing placed 
below the general floor level required the deepening the excavation along 
from the street level. Since the half basement and the concrete floor slab extended 
about (2.75 below the ground surface, the actual depth excavation 
the deepest area was about (11.3 m). During the period excavation, 
tensile cracks were noticed develop. The first set major cracks developed 
parallel to, and distance approx (2.75 from the edge excavation. 
Subsequently, the second and third sets major cracks were noticed distances 
approx and (5.2 and 7.6 from the edge excavation, respectively. 
The widths those crack openings the completion excavation were measured 
10/32, 8/32, and 2/32 in. (8.0, 6.4, and 1.6 mm), respectively. schematic 
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plot indicating the cracking pattern shown Fig. 

Finite Element Analysis and Comparison.—The two records described above 
were analyzed using the finite element method discussed the accompanying 
paper (8). Since the depth excavation was general deeper than 
(12.2 m), [except the 45-ft (13.7-m) depth near the Marshall Street Annex], 
was deemed appropriate that the soil strata considered the analyses 
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FIG. 7.—Direct Shear Test Result (ML) psf 0.0479 
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FIG. 8.—Direct Shear Test Result (SM) psf 0.0479 


confined the top formation consisting primarily approx 80-ft (24.4-m) 
thick layer fine silts and silty sands Fig. 

obtain the strength characteristics the soils the site, direct shear 
tests were performed undisturbed samples under normal pressures equal 
the magnitudes the existing overburden. The strength test results typical 
silty and sandy soils were plotted Figs. and respectively, for both 
the peak and residual values. The values cohesion and friction angle shown 
the figures were obtained from linear regression analysis. can noted 
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Soil Type 
(1) 

elevation 

cohesion peak 393.8 psf 

residual 239.6 psf 

friction angle peak 36.2° 
residual 34.4° 
0.06 
0.08 
103 pef 


Depth of Excavation = 37 ft. 
Length of Reinforcement = 28 ft. 
Diameter of Reinforcement = 4 in. 
Horizontal Spacing = 6 ft. 
Vertical Spacing = 5 ft. 

Angle of Inclination = 15° 


FIG. 9.—Finite Element Grid—Excavation Along 23rd Avenue 0.305 
in. 25.4 mm) 
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FIG. 10.—Comparison Horizontal Movement—23rd Avenue 0.305 
in. 25.4 mm) 
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that the two soil groups have similar strength parameters, and, because the 
lack available triaxial test data determine the hyperbolic parameters, 
single representative soil was chosen from Ref. reflect the average physical 
properties and the hyperbolic parameters the entire soil formation the 
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FIG. 11.—Finite Element Grid—Excavation Along Marshall Street 0.305 
in. 25.4 mm; psf 0.0479 
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site Table The corresponding hyperbolic parameters are 270; 

Lateral Ground Movement.—The finite element grid the cross section 
containing the location the borehole for measuring lateral ground movement 
shown Fig. other pertinent information necessary complete the computer 
input are also listed. comparison the predicted and measured horizontal 
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movements along the length the borehole Fig. 10. The predicted values 
were interpolated from movements the adjacent nodes. Because the inclinome- 
ter measures the ground movement relative the bottom the borehole, 
the actual ground movement can only obtained the reference point 
immobile. The actual lateral ground movement could, therefore, larger than 
those Fig. since the bottom the hole was only about deeper 
than the depth excavation. that level, ground movement result 
nearby open excavation was likely take place. Fig. 10, correction has 
been made shifting the measured curve the right that the reading 
the bottom the hole matches the movement predicted the finite element 
analysis that depth. can seen, both the shape and the magnitude 
the predicted soil movement profiles agree with the measured values. all 
cases, the maximum movement less than 0.32% the depth excavation. 
Tensile Cracks.—The distribution and magnitude lateral distortion the 
ground surface was chosen the basis for predicting the development tensile 
cracks. Lateral distortion, defined the differential horizontal movement 
between two reference points divided the length separating them, measure 
tensile strain the ground surface near excavation. Due the limitations 
the employed soil characterization describing the near-failure soil behavior, 
and the absence realistic knowledge dealing with the frictional behavior 
the soil-slab interface, beyond the present analytical capability predict 
the generation and propagation tensile cracks slabs adjacent open 
excavation. Consequently, the interpretation the analytical results and their 


comparison with the field behavior can only viewed qualitative manner. 

The finite element grid and other information related the analysis are 
presented Fig. 11. surcharge simulating the existing building load was 
applied the ground surface prior the excavation. Fig. shows the distribution 
lateral distortion the ground surface after completion the excavation. 
The maximum lateral distortion occurred distance approx (2.75 
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FIG. 13.—Orientation Reinforcing Members 0.305 
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from the edge excavation; this agrees well with the location the first 
set major tensile cracks recorded. The computer solution, however, not 
capable predicting the locations all the tensile cracks, simultaneously. 
The foundation excavation and the construction the expansion facility 
the hospital were successfully completed. The movements the ground were 
relatively small. While tensile cracks did develop the Marshall Street Annex, 
the building was fully occupied and functional throughout the construction period. 
According the master plan, the brick building was demolished, making 
room for parking facility. Though the field measurement records and the 
soil data were relatively incomplete, the study has shown that the analytical 


methodology employed can adequately capture the fieid behavior the lateral 
earth support system. 
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FIG. 15.—Typical Boring Log—Davis Site 0.305 
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General.—This was field prototype study conducted part the 
comprehensive investigation the situ earth reinforcement lateral support 
system. During the last three years, analytical procedure for evaluating the 
stability the system was formulated (3) and later verified means centrifuge 
model studies (7). addition, method for predicting the field behavior 
the working stress level was also completed reported the accompanying 
paper (8) and Ref. Due the lack available field measurements 
prototype behavior, was concluded that instrumented field study would 
beneficial further substantiate the capability the analytical solutions 
for the prediction field system behavior. The construction and instrumentation 
this sytem were carried out the summer 1979. The study site was 
located the northwest corner the sanitary landfill facility the Davis 
campus the University California. The 30-ft (9.2-m) deep excavation was 
constructed five lifts with each lift approx (1.83 thick. The geometry 
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FIG. 16.—Undrained Strength Results psi 6.9 


the system and the layout the reinforcing members are shown Fig. 
13. The remaining three sides the excavation were left untreated with slopes 
about 45°. The vertical cut had surface area 2,100 (195 and 
total 5,500 (4,215 soil was removed from the open excavation. 
this was small job with limited available funds, the construction methodology 
applied this site was different from the Portland site (9). The drilling and 
grouting reinforcing members were done manually. Instead using special 
fast setting, early strength grout and shotcrete mixtures, concrete mixture 
water-cement-aggregate (sand and pea gravel with maximum size 3/8 
in.) having 1:2.5:8 ratio weight was used grout the holes and gunnnite 
the vertical cut. general view the system near completion shown 
Fig. 14. 

Site Conditions.—The terrain the site was flat, and the soil deposit was 
composed primarily interbedded layers sandy silts and silty clays. Boring 
logs were recorded the time the inclinometer casings were installed. total 
four boreholes were drilled along the center line perpendicular the vertical 


DECEMBER 1981 


: Excavation Depth = 18 ft. 
@ : Finite Element Prediction 
Corresponding to Profile ® 
© : Excavation Depth = 30 tt. 
Finite Elemant Prediction 
Corresponding 


DEPTH BELOW GROUND SURFACE (ft.) 


02 03 


HORIZONTAL MOVEMENT TOWARD EXCAVATION (in.) 


o4 


FIG. 17.—Horizontal Deflection Profile Behind the Excavation) 
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FIG. 18.—Horizontal Deflection Profile (15 Behind the Excavation) 
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cut. These holes were spaced distances 15, 30, and (1.5, 4.6, 
9.2, and 13.7 from the edge excavation. The two holes closest the 
excavation were each (18.3 deep, while the remaining two holes were 
each (13.7) deep. Undisturbed cylindrical samples were taken appropriate 
depths using pitcher sample tube. typical soil profile extending depth 
(18.3 shown Fig. 15. gravelly sand layer and the ground-water 
table were encountered this depth. Later during the excavation, was 
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FIG. 19.—Horizontal Movement Ground Surface Excavation Centerline (By Field 
Transit Survey) 0.305 in. 25.4 mm) 


discovered that the soil deposit the site was extremely variable both the 
vertical and horizontal directions, typical the Putah Creek deposits this 

Undrained triaxial compression tests were performed order determine 
the shear strengths the various soil types encountered. The undrained strength 
associated with these test results are plotted Fig. 16. Further discussion 
the test results will given later. 

Field Instrumentation and Measurements.—Extensive instrumentation was 
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installed the site that the field performance the system could monitored 
and evaluated. Four boreholes with special casings for inclinometer readings 
were installed along the center line the excavation order obtain set 
horizontal displacement profiles. network surface movement markers 
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FIG. 20.—Measured Axial Force Distribution Second-, Third-, and Fourth-Level 
Anchor Rods 0.305 kip 4.45 kN) 


was installed prior the start excavation that the horizontal and vertical 
displacements the ground surface could monitored transit survey. Finally, 
the centermost column reinforcing members were mounted with resistance-type 
strain gages. Four sets gages were spot welded each anchor rod (#8 
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tt. 
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re-bar) that the magnitude and the distribution the axial force each 
member could measured. complete set borehole inclinometer and survey 
readings were taken before and after each lift excavation. Strain gage readings 
were made after each anchor rod had been grouted into place, and after the 
completion each subsequent lift. 

The horizontal deflection profiles yielded inclinometers and are 
shown Figs. and 18, respectively. may observed from these records 
that: (1) The horizontal deflection decreases with increasing depth from the 
ground surface; (2) the horizontal deflection decreases with increasing distance 
behind the vertical cut; (3) the horizontal deflection increases the excavation 
depth increases; and (4) the horizontal deflections tend relatively small 
depths greater than the current excavation depth. 

also interest note that profiles and (Figs. and 18) were 
both taken after the final lift excavation had been completed. Profiles 
were obtained immediately after completion, whereas profiles were recorded 
two weeks later. Thus, there was slight amount short-term movement. 
Additional readings were continued throughout the wet winter and spring months. 
During this period, several earthquakes centered approx miles (120 km) 
southeast the site shook the area. However, movements greater than 
those reflected profile were recorded. 

The record ground surface movements obtained transit survey shown 
Fig. 19. may seen, the survey readings are close agreement with 
the inclinometer measurements. both cases the maximum displacements are 
quite small—approx 0.8 in. (20 mm) the center the uppermost edge 
the excavation. 

The axial force distributions the second, third, and fourth level anchor 
rods are shown Fig. 20. may noted that, for given excavation depth, 
the maximum axial force the various levels approximately the same. The 
shapes force distribution correspond closely those measured Mitchell 
(5) reinforced earth models, and predicted Bang (1), Shen, al. (6), 
using the finite element method analysis, i.e., the force distribution roughly 
triangular, with maximum near the panel facing, and decreasing zero 
the free end. 


Finrre ano Comparisons 


The finite element grid and the locations the boreholes for inclinometer 
readings are shown Fig. 21. Also listed are the dimensions and spacings 
the reinforcement. stated previously, the soil condition the site was 
extremely variable, and would impossible and impractical delineate 
its variations with depth accordance with the boring log information. Further- 
more, there were not enough undisturbed samples determine the hyperbolic 
parameters for each sublayer the soil deposit encountered. average set 
soil parameter was, therefore, used for analytical purposes. choosing the 
parameters, considerations were given the possible sample disturbance and 
strength loss due sampling, handling, and testing. Consequently, more weight 
was placed samples which exhibit stiffer soil properties. The hyperbolic 
soil parameters were determined 548; 0.05; 0.79; and 
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The predicted horizontal deflection profiles inclinometers and 
and behind the vertical cut, respectively) are shown dots Figs. 
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FIG. 21.—Finite Element Grid—Excavation Davis Site 0.305 in. 25.4 
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FIG. 22.—Axial Force Comparison—Second Level Anchor Rod 0.305 
kip 4.45 kN) 


and 18, respectively. For each these locations, the predicted profiles corres- 
ponding excavation depths both and (7.3 and 9.2 are presented. 
all instances, the agreement between the measured and the predicted movements 
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good. may seen, for depths less than about the agreement 
remarkably close. However, for greater depths the finite element predictions 
tend exceed the actual measured ground movements. One possible explanation 
that the actual soil may have been stiffer than indicated the strength 
tests and modeled the selected parameters. 

may observed Fig. 19, the predicted surface movements are likewise 
excellent agreement with the results the field transit surveys. The predicted 
axial forces the second level anchor rod are shown Fig. 22. general, 
the predicted values are about twice large the measured; however, the 
shapes force distribution are similar. The sources discrepancy between 
the measured and the predicted axial forces the reinforcing members cannot 
positively identified. Three possible sources, however, can listed follows: 


The strain gages might have been damaged during the installation the 
re-bars and grouting the holes. 
Strain gage measurements might error caused the drifting 
reference readings under extremely hot weather conditions. 
The bond behavior between the reinforcing members and the surrounding 


soil might not fully captured simplified soil profile assumed the 
analysis. 


Summary 


This paper presents two field performance records lateral earth support 
system. This system uses the situ earth reinforcement technique strengthen 
the native soil, and different from the conventional systems that serve 
retain soils behind vertical cut. The field results were compared with predictions 
obtained from finite element method analysis. general, the predictions 
agree with the field measurement demonstrating that the analytical procedures 
developed can correctly predict the field behavior. Furthermore, has shown 
that ground movements adjacent open excavation can predicted and 
controlled properly designed situ earth reinforcement lateral support system 
installed. The system also economical saving time (30% less time than 
the conventional soldier piling and bracing system for the Portland job) and 
money (the cost approx $10/sq the summer 1979 for the Davis 
job, and 15% below the low bid for conventional system for the Portland 
job). Although the system intended serve temporary support only, 
the fact that additional movements the Davis site were recorded during 
12-month period following the completion construction—a period which 
included wet winter and spring, and several earthquakes—indicates that the 
system also has great deal long-term stability. Thus, the situ earth 
reinforcement lateral support system has shown offer viable alternative 
for providing temporary support deep excavations. 
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CELLULAR COFFERDAM FOR TRIDENT 
DESIGN 


INTRODUCTION 


Extensive waterfront and shore-based facilities are being constructed the 
Naval Submarine Base, Bangor, Bremerton, Washington connection with 
the development new submarine refit and crew training facility for support 
the Trident class submarine. The submarine base located Hood Canal, 
the state Washington, the location shown Fig. The base development 
includes graving drydock which involved the deep water, offshore construction 
steel sheet pile cellular cofferdam. extreme high water (EHW), the 
cofferdam retained (24 water. Interlock tensions associated with 
the size cells needed resist this water depth resulted the use one 
the first United States cellular cofferdam applications high strength steel 
sheet piling and extruded wye connections. the writers’ knowledge, 
one the deepest cellular cofferdams which has been successfully unwatered 
without the use internal earth stability berms. 

The purpose this paper describe the cofferdam and summarize 
important design considerations, many which are unique and represent 
advancement the state-of-the-art for cofferdam design and construction. 
Included are descriptions the cofferdam geometry, site constraints which 
influenced the design, foundation soil conditions, cell fill dewatering used 
limit interlock tensions, seismic design, cell fill densification, and artesian pressure 
relief systems. Additional information concerning cofferdam construction and 
performance are provided Ref. 


The facility was located offshore order permit the continued shallow 
water migration salmon fingerlings. The drydock centerline located the 
Hood Canal, approximately 550 (168 offshore from the mean higher high 
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water (MHHW) shoreline, shown Fig. The drydock and two immediately 
adjacent, pile supported piers form the three-sided Trident Refit Delta. 

Refit Pier No. and the South Trestle were constructed first, and the cofferdam 
was then constructed concurrently with Refit Pier No. The drydock was 
constructed immediately following completion the cofferdam. This scheduling 
accelerated the construction but imposed severe constraints the cofferdam 
design limiting the area available for positioning the structure relative 
the drydock and piers. addition, the offshore location resulted very deep 
cofferdam water depths. 

The following brief description the drydock design and construction 
which provided the basis for the cofferdam geometry and design requirements: 
The drydock basin (27.4 wide 690 (210.3 long. Permanent 
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PACIFIC 
FIG. 1.—Vicinity Map 


laydown areas (22.3 minimum width are required provide for 
work areas, roadways, and service buildings. The laydown areas are El. 
(6.1 provide adequate freeboard above extreme high water (EHW) 
El. 14.6 (4.5 m). Project datum El. mean lower low water 
The top the drydock floor El. —43 (—13.1 which 49.4 (15.1 
below mean tide level (MTL). 

Extensive drydock concept studies were performed evaluate alternate design 
and construction methods. Drydock designs considered included fully relieved 
type, which granular fill and drains would used intercept the seepage 
flow relieve hydrostatic pressures the walls and floor; relieved floor and 
non-relieved wall type; and gravity type, which the walls and floor would 
sized resist the full hydrostatic pressure. Designs considered for the laydown 
areas included pile-supported slabs, ballasted precast concrete cofferdam, and 
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FIG. 2.—Site Plan 0.305 


earth-filled steel sheet pile cellular cofferdam. Drydock construction methods 
considered ranged from tremie underwater floor and wall construction, tremie 


placement the floor, with subsequent in-the-dry wall construction 
unwatered cofferdam, both floor and wall construction in-the-dry 
unwatered cofferdam. 


The decision was made construct full gravity drydock, and have all 
construction performed in-the-dry earth-filled steel sheet pile cellular 
cofferdam. was further established that the cofferdam would permanent, 
except the drydock entrance, provide major portions the required drydock 
laydown areas. The gravity design the drydock required floor thickness 
from (4.9 (5.8 m), depending location. This thickness, 
when combined with in. (460 mm) gravel subbase and in. (150 mm) 
concrete working mat below the floor, resulted final excavation subgrade 
El. —61 (—18.6 El. —64 (—19.5 and maximum water depth 
78.6 (24.0 below EHW. 

When unwatered during construction, the cellular cofferdam had retain 
laydown area for the completed drydock. plan the cofferdam shown 
Fig. and typical section through the structure and future drydock 
shown Fig. Fig. photograph the cofferdam after its pump down. 

The cofferdam basin approximately 774 (236 long and has minimum 
width 130 (40 m). The cofferdam consists circular cells, each 75.83 
(23.11 diameter, and connecting arcs, hereafter termed arcs. The 
cells are spaced from 83.55 (25.47 93.57 (28.52 center, 
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FIG. 3.—Cofferdam Plan 0.305 
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FIG. 4.—Typical Section 0.305 


FIG. 5.—Cofferdam After Pump Down 
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obtain the required geometry and achieve closure. Cells and the in-between 
arcs are permanent. The remaining cells and arcs were removed the end 
drydock construction. All sheet piles are PSX-32 and the connecting elements 
are 40° extruded wyes from the United States Steel Corporation. The plan 
representative cell and arc shown Fig. and the details the 40° 
extruded wye are shown Fig. The manufacturer’s rated interlock strength 
kips per lineal in. (4,900 kN/m). Each cell contains 172 sheets and four 
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FIG. 6.—Representative Cell and Arc 0.305 
ARC 
RADius 


SECTION, 
TYPICAL 


EXTRUDED 
40° WYE 


FIG. Degree Extruded Wye in. 25.4 mm) 


wye sections. The arcs each contain from sheets, depending the 
cell location and spacing. The sheet lengths vary from minimum 
one splice was permitted per sheet, with this splice being the upper portion 
the sheet. During basin unwatering and drydock floor construction the 
cofferdam fill surface was El. (5.5 m). Thereafter, was raised 
required achieve the in-service level El. (6.1 for the drydock 
laydown areas. 


shown Fig. the presence the adjacent refit piers necessitated 
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construction the west side cells immediately adjacent the future drydock 
walls and floor slab, thus precluding the use interior earth berms. Similar 
constraints were imposed the South Trestle the head end the future 
drydock. Cost and environmental considerations directed construction the 
east side cells immediately adjacent the future drydock walls. 


generalized profile subsurface soil conditions, direction perpendicular 
the shoreline, shown Fig. alluvial deposit loose medium 
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FIG. Soil Profile Perpendicular Shoreline 0.305 
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FIG. 9.—Aquifer Piezometric Conditions 0.305 


compact sands and gravels forms the uppermost soil stratum underlying the 
Hood Canal. the cofferdam site, this deposit averages about 
thickness. glacial till stratum, consisting principally very compact 
gray, silty, coarse fine sand sandy silt, underlies the alluvium. The penetration 
resistance typically varied from over 200 blows per (0.3 using 
300 (136 kg) hammer falling in. (762 mm), and 2.5 in. (64 mm) inside 
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diameter split spoon sampler. The glacial till varies from about 
(12 thickness. The average natural water content, liquid, and 
plastic limits are 9.0%, 17.5%, and 13.5%, respectively. Underlying the till 
very compact, glacially overconsolidated deposit interbedded sands and 
gravels with silt layers, herein termed the aquifer sands and gravels. The 
penetration resistance ranged from over 400 blows per (0.3 m), and 
influenced the coarse gravel and cobbles present. Available data indicate 
that the aquifer varies from about 200 250 (61 thickness. 
From five fifteen percent weight the aquifer samples tested pass 
No. 200 sieve. Specific information soils greater depths not known. 
has been estimated (3) that bedrock the site area about El. —2,000 
(—610 m). 

unusual site condition encountered was high artesian pressures the 
aquifer sands and gravels. Shown Fig. are the piezometric elevations measured 
piezometers installed both onshore and offshore. the vicinity the 
cofferdam the piezometric level ranged from El. El. (12 
dramatic variation artesian pressure was observed near the apex the 
refit delta. approximately the end the delta, artesian conditions were 
encountered. believed that the artesian pressures are into 
the deep alluvial deposits this area which are shown the left side 


Fig. analysis the effect these artesian conditions design presented 
later. 


Dredging.—The drydock floor subgrade below the bottom 
grade the canal and much below the surface the dense 
glacial till. Positioning the cofferdam cells immediately adjacent the drydock 
floor subgrade required the installation the sheet piling the basin side 
the cofferdam below subgrade level. 

Because the density the glacial till, was believed that the sheet piles 
could not driven into this layer more than few feet without danger 
driving the piles out interlock. Accordingly, was established that dredging 
would performed within the limits the cells and arcs prior cofferdam 
construction, shown Fig. 10. order limit the total dredged volume, 
the glacial till was excavated in-the-dry from within the drydock floor subgrade 
area following basin unwatering. The typical relationship between the existing 
ground surface and the excavated depth shown Fig. 

Analysis Methods.—The cofferdam was designed general accordance with 
the procedures Design Manual DM-7 (2), except that the maximum interlock 
tensions caused the cell and arc were calculated accordance with the 
TVA secant formula (1). Because the very dense glacial till and the expected 
minimal sheet pile penetrations, the cofferdam was assumed act similarly 
cofferdam founded rock. The point maximum interlock tension was 
assumed occur one-quarter the cofferdam height (H) above the subgrade 
level, H/4. Comments the reasonableness this assumption are made 
Ref. 


The 75.83 (23.11 cell diameter was selected provide adequate safety 
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FIG. 10.—Dredging Plan 0.305 


against overturning, sliding, and shear the fill the cell centerline and 
between the cell fill and sheet piling. Cell interlock tensions obtained accordance 
with DM-7 assume sloped phreatic surface across the width the cells resulting 
from weep hole drainage the fill during and following cofferdam basin 
unwatering. Initial calculations indicated that this loading produced excessive 
interlock tensions. 

major effort was undertaken develop method for reducing the estimated 
tensions least the magnitude which would permit use high-strength 
sheet piling having rated strength kips/in. (4,900 kN/m). Cell diameters 
required for sliding and overturning, fill unit weights and lateral earth pressure 
coefficients, were all re-evaluated with little beneficial effect estimated 
tensions. The only adjustable factor which was found have significant 
effect calculated maximum tensions was the water level within the cells 
and arcs. Calculations disclosed that interlock tensions could reduced 
tolerable levels when the fill water levels were lowered below El. —45 
(—13.7 and maintained essentially flat across the full widths the cells 
and arcs. El. —45 (—13.7 only (4.9 5.8 above the 
excavation subgrade level. Thus, order for the cofferdam technically 
feasible, was necessary devise reliable procedure for obtaining and 
maintaining extremely low water levels within all cells and arcs. was also 
determined that during most the cofferdam basin unwatering phase, the cell 
and arc water levels would have lowered rates such that the water 
levels were maintained below the basin level. 

Fill Dewatering fill dewatering pumped wells was considered 
the only method which water levels could lowered and maintained 
the required elevations. The writers were not aware any similar previous 
application cell dewatering which could referenced guide for estimating 
leakage rates through the interlocks into the fill. The pumps were sized assuming 
leakage rate gpm (1.3 L/s) into the cells following fill dewatering. 
Assuming cell fill specific yield 0.2 and minimal interlock leakage, pumping 
rate gpm (1.3 L/s) also allowed the cells dewatered rate 
least (1.5 per day. 
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The system for dewatering the cell fill consisted two wells within each 
cell (three corner Cells 11, 21, and 24) and one within each arc, the 
wells being in. (150 mm) diameter slotted PVC pipe. The wells were located 
about (2.1 from the sheet piling, shown Fig. 11. Submersible 
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FIG. 12.—Location Plan: Cell Fill Dewatering Wells and Artesian Pressure Relief Wells 


1.5 (1.1 kW) pumps rated gpm (1.3 L/s) 100 (30.5 head were 
installed within in. (150 mm) the bottom each well, shown 
Fig. 11. The pumps were automatically operated on-off probes each well. 
The wells located adjacent the cofferdam basin were not required 
pumped but were provided function observation wells and also standby 
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wells the event excessive seepage breakdown the pumped well (Fig. 
12). Weep holes were also provided and below El. —49 (—14.9 
monitor the dewatering system and serve additional backup. The well 
design made imperative that dredge residue removed from the cells and 
arcs prior filling, and that filling performed such manner produce 
permeable soil within the lower levels near the dredge line. The method 
which this was accomplished discussed Ref. 

addition this pumping system, (4.6 minimum depth blanket 
semi-impervious fill was placed around the exterior the cofferdam, 
shown Fig. reduce underseepage and interlock leakage near the bottoms 
the cells and arcs. 

Seismic Design.—The submarine base located moderately active seismic 
area. Geologic studies disclosed that active faults were present within the 
site area which would require design for ground rupture. However, because 
most the cofferdam will permanent, consideration vibratory ground 
motion was required. Specifically, was necessary prohibit excessive interlock 
tensions due liquefaction the fill and limit cyclically-induced settlements. 
was established relative density approximately 75% would necessary 
prevent liquefaction. relative density could not achieved placement 
uncompacted fill in-the-wet. 

Placement and compaction the fill in-the-dry layers within the offshore 
cofferdam was not possible. Grouting chemical stabilization the soil 
substitute for achieving the 75% relative density was rejected because 
cost and Dynamic compaction with dropped weights was not 
possible, due the danger damaging the cells. Use nonliquefiable soil 
such peastone was rejected because construction seepage considerations 
and cost. 

Compaction vibratory probe procedures was deemed the only feasible 
method. Since (27.4 depth fill required compaction the 
deepest dredge area, consideration was given filling the cells partial depth 
followed probe compaction This would followed one more 
lifts additional fill and compaction. This procedure was rejected because 
interferences from the templates used construct the cells, difficulties 
controlling the work, and other factors. The compaction method selected consisted 
filling the cells and arcs above water level and then compacting the fill 
full depth, with equipment operating from the fill surface. The compaction was 
unusual that vibratory probe compaction depths from (21.3 
(27.4 was required. More importantly, the effects that compaction 
would have interlock tensions was major consideration. Compaction 
cellular bulkhead had been previously performed (4), but precedent existed 
for cofferdam which would subsequently unwatered. 

Choice Cell for use fill the cells and arcs had meet 
the following design requirements: 


High unit weight after compaction provide stability sliding and 
overturning after basin unwatering. 


High internal friction minimize lateral earth pressures and interlock 
tensions. 


Limited percentage gravel permit full depth penetration the probes 
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for vibratory probe compaction and facilitate subsequent augering through 
the fill for bearing pile installation during drydock construction. 

Relatively low percentage fines facilitate probe compaction and 
permit rapid fill dewatering during cofferdam basin unwatering. 


Adequately low permeability limit underseepage near the bottoms 
the cells and arcs. 


The fill specified was well-graded, gravelly, coarse fine sand having 
in. (50.8 mm) maximum stone size, from 15% 45% passing No. sieve 
and 10% maximum passing No. 200 sieve. 

Corrosion systems corrosion protection were provided 
for the sheet piles the permanent cells and arcs. two-coat, shop applied, 
coal tar epoxy-polyamide coating was provided the exterior surfaces all 
sheets exposed the Hood Canal after construction. automatic, impressed- 
current type, cathodic protection system which became operational the end 


drydock construction, was also provided for the exterior surfaces the 
exposed sheets. 


Pressure 


The artesian conditions, shown Fig. posed unique challenges the 
project. Studies disclosed that there was risk uplift failure the glacial 
till during dredging the aquifer piezometric level were not reduced. addition, 
major pressure reduction the aquifer was required permit cofferdam basin 
unwatering. Lastly, permanent partial reduction artesian pressure was required 
for drydock in-service conditions limit uplift pressures the drydock and 
provide adequate seismic stability the offshore slope upon which the facility 
built. Extensive design phase pumping tests from in. (305 mm) diameter 
test well, installed shore, confirmed that would technically feasible 
accomplish the required pressure relief. 

was decided eliminate all excess pressure within the aquifer, relative 
tide water levels during dredging. onshore, artesian pressure relief, pumped 
well system was provided lower the piezometric level El. The system 
consisted five in. (305 mm) diameter wells located adjacent the shore, 
shown Fig. The system was designed for operating capacity 
2,500 gpm (158 L/s). Seventy-five (22.9 long well screens were used 
with the wells installed depths 200 (61 below ground surface. 
Two submersible (18.6 kW) pumps, each rated 750 gpm (47 L/s) 
(27.4 head, were provided each well. 

Because the critical nature the pressure relief and the extremely rapid 
which the aquifer would recover pumping were interrupted, several 
backup systems were provided. These consisted two submersible pumps 
each well capable being operated independently simultaneously, multiple 
generators with automatic transfer and start-up capability, dual power feeders, 
high water alarms each well, and emergency connection the Base power 
system. The system was designed for capacity approximately double that 
expected required for condition with some the pumps inoperable. 

piezometric level El. —67 (—20.4 was required for cofferdam 
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basin unwatering, subgrade preparation, and drydock floor construction. Studies 
disclosed that would necessary transfer pressure relief pumping from 
the onshore system second well system the cofferdam achieve this 
drawdown. series twelve in. (305 mm) diameter wells with (12.2 
long screens was designed for operating capacity 5,300 gpm (334 L/s). 
The pressure relief wells were installed within alternate arcs the cofferdam 
after fill compaction, shown Fig. 12. Two submersible (11.1 kW) 
pumps, each rated 300 gpm (18.9 L/s) 140 (42.7 head were provided 
each well. Redundancy and backup systems similar those for the onshore 
wells were also provided. The system was designed for capacity approximately 
50% greater than that expected required. After basin unwatering, 
additional system open, stand-pipe, subgrade level wells was installed through 
the subgrade provide for controlled venting aquifer pressure case 
pumping failure. The well locations are shown Fig. 12. 

The cofferdam and subgrade level wells were subsequently into 
the drydock system. These provide the permanent level artesian pressure 
relief and backup required for drydock service conditions. 


Summary 


using high unit weight cell fill, high interlock strength steel sheet piling 
and extruded wyes, pumped cell fill dewatering system and pumped artesian 
pressure relief system, was possible design cellular cofferdam that would 


permit the in-the-dry construction graving drydock the required deep 
water location. Since most the cofferdam will permanent part the 
completed facility, vibratory probe compaction the fill was considered 
necessary prevent liquefaction and minimize seismically-induced settlements. 
The dense glacial till the site provided excellent foundation soil for the 
cofferdam, terms both bearing capacity and high strength for providing 
restraint the sheet pile tips. 


FIG. Unwatered Cofferdam 
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Unwatering the cofferdam basin was performed without incident. Fig. 
photograph from within the unwatered cofferdam. Drydock floor construction 
took approximately one year after unwatering. Throughout this period, the cell 
fill dewatering and artesian pressure relief systems performed satisfactorily and 
all aspects the cofferdam performance were well within expected limits. 
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CELLULAR COFFERDAM FOR TRIDENT 
PERFORMANCE 


INTRODUCTION 


The construction the drydock for Trident class submarines the Naval 
Submarine Base, Bangor, Bremerton, Washington included steel sheet pile 
cellular cofferdam. order build the drydock in-the-dry the required 
offshore location, was necessary that cofferdam designed accommodate 
into the final drydock configuration. The cofferdam among the first 
constructed the United States, and because the expected interlock tensions, 
the construction required the use high strength steel sheet piles and extruded 
wye connections. 

The purpose this paper to: (1) Provide summary the cofferdam 
construction procedures; (2) provide data instrumentation that was used 
monitor cofferdam performance throughout construction; and (3) compare mea- 
sured cofferdam performance, specifically cell interlock tensions and deforma- 
tions, with published design procedures. 


plan and typical section the cofferdam are shown Figs. and 
respectively. The structure consists cells, each having diameter 75.83 
(23.11 and connecting arcs. Cells 1-20 and their connecting arcs are 
permanent. The remaining cells and arcs were removed the end drydock 
construction. The subgrade excavation levels within the cofferdam basin required 
for drydock floor construction are shown Fig. Project datum, elevation 
the mean lower low water (MLLW). Extreme high water (EHW) 
elevation 14.6 (4.5 and mean tide level elevation 6.4 (1.9 m). Soil 
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conditions the site consisted loose sand and gravel alluvium overlying 
very compact glacial till which was underlain very dense sand and gravel 
stratum (10). 


Site 


The cofferdam configuration required that the sheet piles driven below 
the drydock floor subgrade the basin side the cofferdam. Due expected 
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FIG. 2.—Typical Section 0.305 


hard driving the dense glacial till, was not considered possible drive 
the sheets from original bottom grades the required tip elevations without 
significant risk sheet damage jumped interlocks. Dredging was therefore 
conducted the cells prior cell construction. 

Cell Construction.—High strength, PSX-32 sheet piles and 40° extruded 
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wyes were used, described (10). Cell construciion began 
Cell and initially proceeded toward Cell then continued from Cell toward 
Cell 20. The remaining cells and arcs were then built with closure being made 
Connecting Arc 22-23. 

The cell sheets were set around two (3.1 high rigid, circular steel 
templates approximately 74.7 (22.8 diameter. The templates were 
positioned four H-pile studs driven into the foundation soils. The vertical 
distance between templates varied from (3.1 (9.3 m), with 
the greater spacing being used the deeper cells. The cell design radius was 
achieved the upper template placing wood blocks between the template 
and the sheets. attempt was made block the sheets the lower template. 

Sheet piles were set commencing each the four connecting wyes. Four 
closure points per cell were made this way, one near the center each 
cell segment between wyes. Similarly, closure was made near the central portion 
each connecting arc. Two arc sheets were threaded the wyes prior 
driving the wyes the cell sheets. 

Sheets were driven with MKT 10B3 singie acting hammer having rated 
energy 13,100 ft-lb (17.8 kN-m). Driving was terminated when driving 
resistance blows per inch (25.4 mm) was attained for sheets driven 
pairs. Alternatively, driving was stopped the penetration pair sheets 
exceeded (1.8 below the dredge level and the driving resistance was 
equal blows per inch (25.4 mm). Sheets typically were driven pairs 
penetration least (0.6 with most sheets being driven about 
(1.2 m). Positioning the cell template sheet pile threading, and driving usually 
required ten 10-hr shifts per cell. 


Following completion sheet pile driving cell, the cell was mucked 
out with smooth-edged clamshell bucket. diver inspection was made 
verify that the sheets were not damaged nor out interlock and that essentially 
firm, clean bottom existed prior filling. second diver inspection the 
interlocks was made after filling. Cell fill consisted well-graded, gravelly, 
coarse fine sand having in. mm) maximum particle size and not 
more than 10% weight passing No. 200 sieve (Fig. 3). The fill was placed 
clamshell bucket manner that displaced any remaining dredging sludge 
away from the locations the cell fill dewatering and observation wells which 
were installed later. All fill placed below water was deposited lowering the 
bucket the surface the previously placed fill before opening. Typically, 
the cells were filled prior placement fill within adjacent connecting arcs. 
all cases, however, the soil within the arcs was maintained below 
the fill level the adjacent cells. order prevent tilting other detrimental 
movements individual cells and connecting arcs during filling, the fill surface 
adjacent the sheets was maintained approximately level (within range 
(1.5 m)). 

The fill all permanent cells and connecting arcs was densified minimum 
75% relative density prevent seismically induced liquefaction, minimize 
settlements, and obtain the necessary effective weight required for cofferdam 
stability (10). Fill the temporary cells and arcs was similarly densified 
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provide uniformity soil properties for cofferdam stability. The depth fill 
compacted ranged from (21.3 within portions the eastern and southern 

The contractor used vibratory probe manufactured Toyomenka (America), 
Inc. The probe was in. (559 mm) diameter and approximately 100 (31 
long. Penetration the probe and compaction was achieved combination 
vertical and horizontal vibrations generated vertically mounted vibratory 
pile driver, VM2-25000A (1,450 ft-lb, 1.96 the top the probe, 
and horizontally mounted vibratory pile driver, (360 ft-lb, 0.49 
kN-m), the probe tip. Penetration the probe was aided jetting with 
two in. mm) diameter jet pipes which discharged vertically. The probe 
was vibrated full depth for min. Subsequent compaction was performed 
consecutively raising the probe (1.8 m), holding momentarily, and 
then reinserting (0.9 until the probe was withdrawn the surface. 
compaction proceeded, fill was added placing additional soil around the 
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FIG. 3.—Fill Gradation Cells and 


probe. Penetration the probe (24 depth typically took 5-15 
min. The penetration and compaction cycle one individual probe location 
required about one hour. 

Various center center probe spacings were evaluated during trial compaction 
program Cell order establish procedure compact the remaining 
cells. The next six cells were compacted using probes located about 
(3.7 center center probe spacing. This initial spacing was found 
provide level compaction excess that required. Therefore, nineteen 
probes were used within the remaining cells, all about (4.6 spacings. 
The sequence compaction for each pattern shown Fig. Compaction 
results were measured the standard penetration test, with results correlated 
relative density, accordance with Gibbs and Holtz (4). Typical results 
are shown Fig. Near surface, in-place unit weight tests were also performed 
each the cells. Soil moist unit weights after compaction ranged from 118.2 


pef (1,894 139.5 pef (2,235 kg/m’) with average 130.3 
(2,087 
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Compaction the connecting arcs was conducted following compaction 
the adjacent cells. The number probes within each arc varied depending 
the area but typically ranged from four six. 
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FIG. 4.—Vibratory Probe Patterns 0.305 
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the first four cells compacted, high pore water pressures were observed 
develop within the fill during probe insertion and compaction. Momentary 
excess heads (15 more were observed piezometers located 
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close the probe locations. Increases interlock tension were also observed 
associated with development the compaction pore pressures. 

pore pressure relief system was designed during the construction phase 
and installed the remaining cells limit excess pore pressures and aid the 
compaction facilitating water flow from the soil. Two types pore pressure 
relief systems were utilized: 


gravity system, consisting six in. (51 mm) diameter slotted pipes 
installed for the full depth the cell and located adjacent the sheets. During 
compaction operations, water freely flowed from the pipes and was diverted 
overboard. This procedure was used during compaction the next ten cells. 

system, consisting four in. (203 mm) diameter wells installed 
for the full depth the cell adjacent the wye connections. The wells were 
equipped with submersible pumps which were operated maintain water levels 
the wells about elevation —10 m). This procedure was used during 
compaction the final ten cells. 


After fill compaction, wells were installed within each cell and connecting 
arc dewater the fill order maintain interlock tensions within tolerable 
limits during and subsequent cofferdam basin unwatering (10). addition, 
drilled weep holes were installed and below elevation —49 (—14.9 (10). 
During the construction phase, was established that additional pumps should 
installed the cell observation wells. 

Both wells within each cell were pumped during basin unwatering. Intermittent 
pumping from only one well per cell maintained adequate drawdown thereafter. 
Measurements disclosed that the dewatering system maintained hydraulic gza- 
dients across the cells, ranging from near zero less than (1.5 m). 

During initial cofferdam basin unwatering, the fill water levels were maintained 
below the basin level. Following basin unwatering below elevations 
—45 (—13.7 m), the fill water levels were maintained between elevation —46 
and elevation —61 (—18.6 until the drydock floor was completed. 
The cell and arc water levels were incrementally raised subsequent drydock 
construction proceeded. 


Basin 


The cofferdam was closed tide level about elevation (2.7 m). The 
contractor was permitted lower the basin elevation without restrictions. 
Thereafter, the basin was lowered rate greater than (1.5 per 
period. 

each increment below elevation the basin water level was 
constant for day while the cofferdam performance monitoring instrumenta- 
tion was read, diver inspection was made the interior and exterior perimeter 
the cofferdam, and the data were reviewed evaluate cofferdam performance. 

The basin level was maintained elevation —20 (—6 for three week 
period while the cells were proof tested. The proof loading was obtained 
temporarily raising the water level the cells and connecting arcs elevation 
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Typically, two three cells and the adjacent connecting arcs were proof 
tested one time. The estimated maximum interlock tension under this loading 
was comparable that expected during the final unwatering stage. 


Cofferdam instrumentation consisting vibratory wire strain gages monitor 
interlock tensions, and inclinometers measure sheet pile deflections, was 
installed the locations shown Fig. Also, optical survey measurements 
were made reference points established the tops two sheets each 
cell. 

total 112 strain gages were installed two sheets Cell and six 
sheets Cell Locations the instrumented sheet piles are shown Fig. 
The gages were installed four levels, shown Fig. 7(a), except for 
Sheets 5-4 and 5-6 where they were provided only the lower three and two 
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FIG. 6.—Locations Instrumentation 


levels, respectively. The gages were placed the sheets before sheet pile 
threading and driving levels such that, after driving, the gages were 
approximately elevation —20 m), —40 (—12 m), and —58 
m). This placed the gages near the dredge level and intervals approximately 
one-quarter the cofferdam height above. Additional installation details are shown 
Fig. 8(a). 

Inclinometers were installed the centers the basin and Hood Canal sides 
Cells and 19, and the basin side Cell shown Fig. 
Installation details are shown Figs. 7(b) and 8(b). Horizontal and vertical 
sheet movements for each cell were monitored using optical survey points 


the top two sheets, one each the center the basin and Hood 
cell segments. 


Interlock tensions measured after filling and compaction are shown Figs. 
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and for Cells and respectively. The following comments apply: 


Cell was the cell which trial compaction operations establish 
production probe spacings took place. such, the cell compaction was some- 
what atypical with total probes being made versus the probe 


arrays utilized the other cells. 
LL 
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During tide cycles, interlock tensions varied kip/in. (175 kN/m), 
with the higher tensions occurring low tide. All data presented herein are 
for conditions approximately low tide. 


The limit lines shown Figs. and represent the writers’ estimate the 
probable total range interlock tensions the basin and Hood Canal sides 
the cells. The causes the abnormally low measured tensions Sheet 
5-6 are not known. Cell consistently exhibited higher tensions than Cell 
Reasons for this are not certain but possible factors include: (1) Differing 
compaction patterns; (2) variations permeability and gradation the cell 
fill (Fig. 3); and (3) effects sheet embedment tensions measured the 
dredge level. Cell embedments and near the instrumented sheets ranged 
from about 4.0-7.5 (1.7-2.0 m), whereas Cell embedments were from 
about 2.5-4.5 (0.8-1.4 m). significant differences tensions between 
the sheets the cofferdam basin side and those located adjacent the wyes 
along the common wall were observed Cell except for the data after 
filling sheet 5-5 for the elevation —20 (—6 and elevation —40 (—12 
levels. seen that typically compaction caused increase tension 
from near about kips/in. (350 kN/m). 


Interlock tensions measured subsequent basin unwatering are shown 
Fig. 11. The limit lines shown for Cell are for the sheets the basin side 


(5-1 and 5-2). The limit lines for Cells and indicate relationship between 
the two cells similar that observed prior basin unwatering, with the tensions 
continuing higher Cell For example, the post-unwatered limit line 
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tensions Cell elevation —40 (—12 are 7.2-8.6 kips/in. 
whereas Cell the range from 5.5-7.2 kips/in. 

Changes interlock tension which occurred Cell elevation —58 
(—17 subsequent basin lowering are shown Fig. 11. These changes 
are attributed removal loose fill berm which had accummulated against 
the basin side the cell elevation —49 (—15 m). Sheet 5-2, the tension 
the elevation —60 (—18 upper level strain gages increased 1.65 
kips/in. (290 kN/m) value 10.2 kips/in. (1,790 kN/m). Other lesser 
changes occurred indicated. The relatively small increases observed Sheets 
5-3 and 5-5 suggest that significant resisting frictional forces had developed 
between the sheets and the cell fill along the common walls near the wyes. 
The data also indicate that even relatively loosely placed interior soil berms 
can beneficial limiting interlock tensions. 

Following unwatering, some differences were observed between the Cell 
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interlock tensions the sheets the basin side and those the common 
wall side the wyes. the elevation level, the common wall tensions 
were essentially identical the lower limit the tensions the basin side. 
the elevation —20 (—6 level, the average the common wall tensions 
was about 20% higher than the average the basin side sheets. the elevation 
—40 level, significant difference measured tensions existed 
between the two common wall sheets. However, the average for the common 
wall sheets was essentially equal the average for the basin side sheets. 

When the cofferdam basin and cell fill water levels were lowered, the measured 
interlock tensions the cell sheets the basin side increased, whereas the 
tensions the Hood Canal side decreased. Typical changes are shown Fig. 
12. Lowering the water levels increased the lateral pressure the basin 
side the cell. However, the lowering the cell water level decreased the 
lateral pressure the Hood Canal side. The direction the measured changes 
compatible with the expected behavior. 
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Data obtained during cofferdam proof testing are also included Fig. 12. 
The spikes the tensions correspond the maximum test water level the 
cell fill. The tensions after proof testing, when the cell and pool water levels 
were both elevation —20 (—6.1 m), were somewhat larger than those measured 
for this condition prior the test. This attributed slight increase 
cell diameter under the maximum test water level the cell fill and the 
accompanying small radial expansion the cell fill. The sheets experienced 
only partial elastic shortening following removal the test load because full 
recovery was prevented the cell fill which had previously been displaced 
outwards against the sheets during the proof loading. The data Fig. indicate 
that this prestressing effect was essentially overcome the time the basin 
water level was elevation —30 (—9 m). 


Inclinometer profiles for Cell are shown Fig. 13. These data are typical 
those for the other inclinometers. Since the inclinometer casings were not 
installed until after cell filling, the as-driven sheet profiles and the sheet movements 
during cell filling are not known. The data show the sheet profiles relative 
the design locations which are referenced the cell centerline. 

seen that the sheets were significantly out plumb. Data from the 
nine inclinometers revealed that the instrumented sheets ranged from 12.7-54.5 
in. (323-1,380 mm) out plumb, top bottom. The data consistently showed 
that the tops the sheets were outside the design locations and the bottoms 
were inside. For example, Sheet I-SE was about in. (400 mm) inside 
the dredge level, and about in. (480 mm) outside the top. Compaction 
the fill resulted movement the Hood Canal and Basin sides the 
cell away from the cell centerline, producing somewhat oval shape. The data 
have been plotted the assumption that the sheet pile tips were fixed. 

Because the degree out-of-plumbness measured the inclinometer sheets, 
extensive diver plumb line and dredge line surveys were made each cell. 
These surveys disclosed that the dredge line, the sheets near the wyes were 
typically outside the design radius, whereas those midway between the wyes 
were inside the design radial location. This believed the result the 
cell construction procedures wherein blocking was only used the upper template 
level. The sheets initially set near the wye were positioned with the bottom 
outside the radial line; the remaining sheets had toed complete 
closure. All inclinometers were located sheets that were near the closure 
points where the sheet tips were inside the design radius. 

Movements the Cell inclinometers during unwatering are shown Fig. 
The displacements shown are incremental movements for the various stages 
unwatering from the sheet pile positions which existed when the basin was 
elevation 5.3 (1.6 m). Movements resulting from proof testing the cell 
basin level elevation —20 (—6 are included. The causes the 
movements during proof testing are not fully understood, but are believed 
be, least partially, the result creep movements which occurred during 
the three week period that the basin was elevation —20 (—6 m). 

Total movements measured the top the cell were very small, with the 
movements both inclinometers being less than in. (76 mm) during unwatering 
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elevation —59 (—18 m). The data indicate substantially rigid body rotation 
the cell with slight bulging the basin side. 

Optical survey measurements the movements the top the west side 
cells, function cofferdam basin level, are shown Fig. 15. The data 
shown are for the reference points the basin sides the cells. The largest 
total movement the end unwatering was slightly less than in. (102 mm) 
Cell Movements the east side cells were generally equal less 
than those for the deeper west side cells. The largest incremental movement 
for each unwatering interval occurred from elevation —0.3 —20.0 
m). This explained the manner which the net hydrostatic 
overturning moment the cofferdam developed. For equal changes basin 
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level, the largest incremental moments occurred during the initial stages basin 
lowering. 

Vertical movements the reference points were also checked. With the 
basin level elevation —59 m), movement the reference points 
indicated average heave 0.5 in. (12.7 mm) the Hood Canal side, and 
settlement 0.5 in. (12.7 mm) the basin side. 

Movement data for the top Cells and are shown Fig. 16. Data 
for Cells and represent the approximate limit movements within the 
cofferdam. The data Fig. show that creep type deformations continued 
after the cofferdam was unwatered. The rate movement for Cells and 
decreased with time, whereas the creep rate for Cell remained approximately 
constant during the observation period. Final excavation for, and preparation 
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the drydock floor subgrade occurred during this period. the end this 
time, in. (460 mm) depth compacted fill, and in. (150 mm) concrete 
working slab had been placed. 


Interlock Tensions.—Several methods for estimating maximum sheet pile 
interlock tension are presented the literature. Interlock tensions calculated 
using several the sources for the geometry Cell are summarized 
Table Measured interlock tensions Cells and are shown Table 


Soil parameters and design assumptions recommended the several sources 
vary follows: 


The maximum interlock tension computed above the dredge line heights 
which vary from zero H/3. 

Lateral earth pressure coefficients vary from the active condition, 
0.28, Krynine value, 0.45, with empirical values most often used. 

For the Cell geometric configuration, the effect the connecting arc 
loading the common wall sheets next the wye assumed increase 
tensions from zero high 50% the main cell tension. 


The writers’ interpretations the cofferdam data are follows: 


basin side sheets, the interlock tensions both before basin unwatering, 
(Figs. and 10) and after (Fig. 11), indicate that maximum tensions occurred 
approximately above the dredge line instead the commonly 
used H/4 level, which, for the instrumented cells, approximately 
above the dredge line. 

For the common wall sheets adjacent basin side wyes, the after filling 
and after compaction data indicate that the maximum tension also occurred 
approximately above the dredge level. However, for these sheets 
after unwatering, the point maximum tension higher elevation. This 
believed the result the restraint provided the connecting arc 
fill, which limited deformations the common wall sheets near the bottom 
the cell during basin unwatering. 

Values shown Fig. have been calculated based in-place 
soil densities and the limit lines measured interlock tensions indicated 
Figs. 10, and 11. After filling (Fig. values for the loose fill were 
relatively constant with depth both cells, ranging from low about 0.2 
high about 0.4. After compaction, (Fig. 17(b)), values increased 
with depth and ranged from low about 0.2 elevation high 
Cell 0.46 the elevation —40 (—12 level. Conversely, after basin 
unwatering, (Fig. 17(c)), values decreased with depth, ranging from 0.3-0.45 
the elevation level, and from 0.23-0.36 the elevation —40 (—12 
level. For the three conditions, the values for the after-filling condition 
are believed the most representative values the lateral earth 
pressures within the cells, because the absence prestress effects the 
sheets. The values after compaction are believed have been affected 
residual prestressing the sheets due compaction operations. Residual 
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TABLE Interlock Tensions Following Basin Unwatering Using Various 


Design Procedures 


Assumed 
point 
maximum 
interlock 
Reference tension 
(1) (2) 
TVA (3) 


Lacroix, 
al. (6) 
Navy Design 
Manual (7) 
Swatek (11) H/4 


Terzaghi (12) H/3-H/4 


Corps Engi- 
neers (13) 


dredge level 


pr. 


Lateral 
earth 
pressure 
coefficient, 


Calculated Interlock 
Tension, kips per inch 
wye factor 
(3) (5) (6) 
active 
0.28 
empirical 
0.40 
Krynine 
0.45 
empirical 
0.25-0.35 
empirical 
0.4 
active 
0.28 


Note: kip/in. 175 kN/m; 38° Assumed. 


TABLE 2.—Measured Tensions Following Basin Unwatering 


Average Range 
Measured Values, 
kips per inch 


Lateral 
Point earth 
measured pressure 
interlock coefficient, wye 
Location tension connection Remarks 


(1) (2) (3) (6) 


0.29-0.36 7.4-8.8 Fig. 
(estimated) 
5.5-7.6 
5.8-7.6 
(estimated) 


0.23-0.31 Fig. 


“Magnitude maximum interlock tension based interpretation data. 
insufficient for estimate. 
Note: kip/in. 175 kN/m. 
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prestresses developed result compaction pore pressures, soil displacements 
due probe insertion, and fill densification. The values after unwatering 
were controlled the net effect changes vertical effective stress and 
interlock tension. The pumped dewatering the cell and connecting arc fill 
resulted increased vertical effective stress, with the largest increments occurring 
the lower levels the cells and arcs. the lower levels, the magnitude 
the increased effective stress was sufficient counteract all most 
the prestress effect. The range values for the after-filling and after-unwater- 
ing conditions are almost identical elevation —20 (—6 Cell and 

For the unwatered condition, the back calculated values are within 
the range values recommended the literature. The values Fig. 
show that not constant and varies result construction procedures. 

The data for Cell show that the interlock tensions the common wall 
sheets next the wyes were from zero 20% higher than measured the 
basin side. Thus, the procedures Table utilizing the TVA secant formula 
account for the arc effect appear conservative. Published field data the 
connecting arc effect are available for cellular bulkhead where common 
wall sheet close T-pile connection was instrumented, (14). After filling 
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FIG. 17.—Range Lateral Earth Pressure Coefficients Cells and (a) After Filling 
After Compaction (c) After Basin Unwatering 0.305 


behind the connecting arc was completed, the common wall tensions were 
generally found equal less than the tensions the cell sheets 
front the T-pile connection. Recent studies for another cellular bulkhead 
(9) showed that the common wall tensions close the connection were also 
equal less than those the cell sheets front the connecting member. 


Cofferdam Movements.—Swatek (11) states that cofferdams commonly move 
inward toward the basin amount equal about the cofferdam height 
during unwatering. This would equal about 9.5 in. (240 mm) for this 
cofferdam. Actual movements ranged from 1.5-4.0 in. (38-100 mm), which 
0.16-0.42% the height. The reduced deflections are considered principally 
the result the quality the fill and the compaction operations which increased 
the fill density and shear modulus. These observations are consistent with those 
previously reported Rimstad (8) where movements model sheet pile 
cofferdams were measured. The study found that laboratory models having 
dense soil deflected less than those which the soil was less dense. 

Prior cofferdam unwatering, prediction cell movements was made 
Forrest (2) using computerized procedure patterned after method proposed 
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Brown (1). The shear modulus the fill with depth was based assumed 
75% relative density. The predicted movement, labeled Fig. 16(a), 
within the range the observed inward movements. The calculation assumed 
rigid foundation for the cells, which appears reasonable since the inboard 
sheets were observed move only about 0.5 in. (12 mm) vertically during 
unwatering. 

Many investigators, e.g., Swatek (11), Terzaghi (12), and others, report that 
following unwatering, maximum bulging the cells occurs approximately 
H/4 above the dredge level. From these observations has been concluded 
that the point maximum tension similarly occurred H/4. The inclinometer 
data Figs. and 14, and the interlock tension data, not support these 
observations. The profiles show that the cells did not experience the degree 
bulging that common most cofferdams. believed that the high 
quality the fill and the densification may have contributed this occurrence. 


The following principal conclusions are made: 


High strength sheets and extruded wyes can successfully utilized 
construct steel sheet pile cellular cofferdams suitable for resisting (24 
water depths. 

The data show that substantial interlock tensions developed near the dredge 
level. For the basin side sheets, the maximum interlock tension believed 
have occurred approximately above the dredge level. The 
common assumption that the maximum interlock tension occurs H/4 above 
the dredge line for the basin side sheets may not valid for deep cofferdam 
cells having shallow sheet pile embedments. 

The data indicate that the elevation maximum tension for the commcn 
wall sheets also about above the dredge level prior basin 
unwatering. However, believed that because the restraint provided 
the connecting arc fill, the maximum tensions after unwatering occurred 
greater distance above the dredge line. Insufficient data exist establish 
the specific elevation, and use the distance for the common 
wall tensions cofferdams this type considered prudent. 

Lateral earth pressure coefficients can variable with depth, depending 
construction procedures. Therefore, empirical values should used design 
rather than coefficients based highly theoretical considerations. For compact- 
ed, well graded sand and gravel fill, values the range 0.35-0.40 are 
considered appropriate for computing sheet pile interlock tensions following 
basin unwatering. 

For the 40° extruded wye connection, the data indicate that the common 
wall tension near the wye not significantly higher than the tension the 
basin side sheet piles. allowance 15% for increased tension considered 
reasonable factor. 

The quality and density the cell fill have significant influence 
the magnitude cell movements during cofferdam basin unwatering. 

deep, large diameter cells, two more template levels, including one 
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near the dredge level, should used ensure sheet pile plumbness 
and geometric control. 


The successful construction the cofferdam was largely dependent 
the high bearing capacity and strength the glacial till the site which provided 


nearly rigid foundation and large lateral resistance the embedded sheet 
piles. 


The facility was constructed for the Department the Navy, Commander, 
Naval Facilities Engineering Command through the Officer Charge Con- 
struction, Naval Facilities Engineering Command Contracts, Trident. The design 
was Fay, Spofford Thorndike, Inc., Boston, Massachusetts. The geotechni- 
cal consultant was Haley Aldrich, Inc., Cambridge, Massachusetts. The 
construction contractor was Willamette-General-Manson, Joint Venture 
Willamette-Western Corporation, Portland, Oregon, General Construction 
Company, Seattle, Washington and Manson Construction and Engineering 
Company, Seattle, Washington. The authors wish express their appreciation 
the engineers the staff the Officer Charge Construction, Naval 
Facilities Engineering Command Contracts, Trident for their contributions 
the successful completion the project. 
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COMPRESSIBILITY AND BEARING CAPACITY 
Nabil Ismael," ASCE and Aleksandar ASCE 


INTRODUCTION 


The analysis ultimate loads shallow foundations usually employs the 
well-known Terzaghi equation (6), which contributions the bearing capacity 
from different soil and loading conditions are summed. This analysis based 
solutions developed for the rigid-plastic solid the classical theory 
plasticity. presently known, this solid assumed exhibit deformation 
prior shear failure and plastic flow constant stress after failure. Thus, 
the capabilities theoretical predictions ultimate load are limited relatively 
incompressible soils the general shear failure mode (8). For compressible 
soils, common practice use the available solutions with possible reduction 
the strength parameters for the effects compressibility. 

Compressible soils fail usually either local punching shear. These modes 
are characterized substantial strain prior failure and limited slip surface 
ending the soil mass below the footing failure. The failure mode that 
can expected depends the relative compressibility the soil the particular 
geometrical and loading conditions. Among the principal factors governing the 
failure mode are the soil compressibility, relative density, the foundation size, 
and the presence overburden pressure. 

Because the lack rational methods for analyzing bearing capacity failures 
for compressible soils, Terzaghi proposed the use the bearing capacity equation 
for the general shear failure mode with reduced strength parameters and 
defined 2/3 and tan 2/3 tan Although this proposal simple 
apply, not supported either documented laboratory model tests 
field tests. also suggests jump bearing capacity transition general 
shear failure, phenomenon which, course, does not occur. Furthermore, 
apply this method, the practicing engineer must have data available the 
load-settlement behavior the soil judge the probable failure mode. Such 
data are not always available. 

Proposals this kind may useful practice; however, they are based 
doubtful premise that the effects the relative compressibility soil 
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under different geometrical and loading conditions can accounted for 
adjusting and different solution proposed recently Vesic (8) based 
extension the cavity expansion solutions (5,7) include shallow 
foundations compressible soils. apply this solution, the rigidity index, 
I,, defined the inverse the shear strain the soil, first determined; 
then the compressibility factors are calculated for use the bearing capacity 
formula. These are reduction factors which vary with the angle shearing 
resistance and the parameter, The method emphasizes the importance 
the deformation parameters the soil, the calculation the compress- 
ibility factors. 

This paper examines the effect soil compressibility the ultimate bearing 
capacity shallow foundations means small-scale model tests (4) carried 
out two frictional soils having identical shear strength characteristics, 
but different deformation characteristics, The tests include surface 
tests well tests with initial overburden pressure varying psi 
(103 strip and circular footings. 

Experimental results are compared with the available analytical methods. The 
influence soil compressibility the load-settlement behavior, and bearing 
capacity assessed light test results. 


Son 


The test medium, Chattahoochee sand, subangular medium sand containing 
traces mica. Mixing with 10% (by weight) mica gave homogeneous 


TABLE 1.—Classification Properties the Test Soils 


Property Sand plus percent mica 
(1) 
Coefficient uniformity 
Mean diameter 
Specific gravity 
Maximum density 
Void ratio maximum density 
Minimum density 
Void ratio minimum density 


TABLE 2.—Specifications Number 60, “dry ground” 


Passing number Maximum specifications, percentage 


(2) 


“United States Gypsum Company. 
States Standard Sieve. 
Note: Minimum density 


(1) 
max 
100 max 
200 max 
325 max 
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compressible material. Its physical properties and that Chattahoochee sand 
are shown Table The mica used dry roofing grade having the specifications 
Table 

Fig. the two soils have almost identical grain size distribution curves. 
Triaxial compression tests with volume change measurements were performed 
air dry samples both soils initial relative density 70%. This 
same relative density was achieved vibration the model tests described 


% FINER 


LEGEND 


Mc 


MICA CONTENT 


FIG. Size Distribution Curves 


A single line for all mica contents 


5 


% at failure (psi) 


6.9 


the next section. Because the curvature the Mohr envelope, the secant 
angle shearing resistance, varies with the mean normal stress failure, 
line uniquely determined this variation for the two soils employed 

While Fig. shows identical shear strength characteristics for the two 
cohesionless soils, their deformation characteristics varied considerably with 
the addition mica. For example, the initial modulus, E,, determined from 
triaxial test, o,, versus axial strain, plots decreased nearly half 
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Oo MC «0% 
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adding mica. Table summarizes the strength and initial modulus values 
determined from the triaxial tests. Values Poisson’s ratio were calculated 
means 


v= 


which Poisson’s ratio; the volumetric strain; and incremental 


TABLE 3.—Results Triaxial Consolidated Drained Compression Tests 


Initial Properties 


Mica Density, 
asa pounds pounds Secant pounds 
per- per Rela- per angle, per 
cent- square tive square square 
foot density inch 


(4) 


change. Values calculated from the initial portion the stress-strain curves 
varied narrow range 0.3-0.33. 


Significant parameters the laboratory model footing test program were soil 
compressibility, initial overburden pressure, and foundation shape. Assuming 
that the influence the soil above the base level the footing can replaced 
uniform surchage, the effect overburden pressure was simulated 
air pressure transmitted around the footings means rubber bags similar 
arrangement that reported Beer (3). This assumption the safe 
side and implies that the shearing resistance the overburden soils above 
the foundation level ignored. Three initial overburden pressure equal 
and psi (0, 34, 103 were used, the largest value equivalent 
nearly (6.1 overburden soils and represents transition deep 
conditions Width for most foundations. rectangular strip footing 
measuring 2.5 in. (63 381 mm) and 4.3-in. (109-mm) diameter circular 
footing were employed. Duplicate tests were performed and the average results 
were taken. total tests were completed for each footing the same 
relative density 70%. The footings had rough bases achieved attaching 
sand paper the base. 

Test Equipment.—The program used sand box measuring 


Test 
num- 
ber 
0.70 2.9 47.8 1,210 
93.7 0.69 30.9 38.0 5,000 
0.71 2.8 46.4 585 
87.8 0.71 11.3 40.8 1,350 
85.9 0.65 31.4 38.4 2,360 


~ 
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in. depth (1,143 1,143 457 mm). rigid frame mounted the box provided 
reaction. manually operated hydraulic jack (20-ton capacity, 5.06-in. stroke) 
applied the load. Loads were measured proving ring connected the 
jack and the footing. Figure shows general view prior test with overburden 
pressure. Restraining steel plates attached the box enabled application 
overburden pressure inflating rubber bags connected the steel plates. 

regulator and test gauge with pressure range 0-30 psi (0-207 
supplied air under controlled pressure the rubber bags through quick fittings 
the steel plate. The pressure gauge (Fig. has accuracy +0.25%. 

Density measurements were made each test the sand cone method. 
Penetrometer soundings were also performed for density control. small 0.5-in. 
(12.5-mm) diameter penetrometer was employed for this purpose. 

Experimental Procedure.—To fill the box with sand having the desired relative 
density 70%, equal weights sand were placed, levelled, and vibrated for 
three min form 2-in. (51-mm) layers. V-35A (Syntron Company) vibrator 


3.—General View Prior Test with Overburden Pressure 


mounted 0.25-in. (6-mm) thick steel plate achieved the same relative density 
for both soils using nearly the same vibrating effort. 

tests with overburden pressure, the steel plate equipped with rubber sheet 
was tightened the top the box, and air pressure applied the desired 
level. Slow jacking enabled the footing touch the soil surface and two dial 
gages mounted the footing measured vertical displacements. 

Equal increments stress were applied and maintained for least min 
until all movements had ceased based recorded deflection readings. Each 
loading increment was approximately equal the expected ultimate 
pressure. Tests with high overburden pressure psi (103 required 
additional jacking maintain constant stress. 


Anatysis Discussion Test 


Ultimate Load Criterion.—Two ultimate load criteria were evaluated. The first 


¥ 
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defines the ultimate load the point which the slope the load-settlement 
curve first reaches zero minimum value. Another ultimate load criterion 
defines the ultimate load the point break the load-settlement curve 
log-log plot (3). Both criteria require that the loading test carried 


Pressure (psi) 


Settlement (in) 


5in—+| 


Footing pressure (psi) 
100 150 


Settlement (in) 


FIG. 5.—Pressure-Settlement Diagrams for Strip Footing, psi in. 
mm; psi 6.9 


large displacements, exceeding 25% the foundation size. Because the log-log 
method more consistent was selected define the ultimate bearing capacity 
for these model tests. 
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Angle Shearing Resistance.—The secant angle, determined from triaxial 
tests was used the analyses irrespective the shape the footing. Since 
this angle varies with the mean normal stress along the shearing surface, the 


Footing pressure (psi) 
200 300 


Settlement (in) 


FIG. Settlement Diagrams for Strip Footing, psi in. 25.4 
mm; psi 6.9 


FIG. 7.—Ultimate Load Based Log Load/Log Settlement Plot—Strip Footing, 
psi 


analyses employed Beer’s (2) recommendation that the bearing capacity 


should made using strength characteristics corresponding average mean 
normal stress equal 
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which the ultimate pressure; the initial overburden pressure; and 
the tangent friction angle being constant high stress levels. 

Test Results.—Figs. and show pressure-settlement diagrams for the 
strip footing with initial overburden pressure and psi (0, 34, 103 
respectively. Figs. and show dimensionless log-log plots the 
data with arrows the locations failure loads. Tests circular footings 
are not presented because they gave similar conclusions. 

examination Figs. indicates the remarkable effect compressibility 
the reduction the bearing capacity and increased settlement. The arrows 
the curves show the failure load determined the log-log criterion. 


Q/ayB 


FIG. Load Based Log Load/Log Settlement Plot—Strip Footing, 
psi 


The slope the initial straight line portion the curves greatly reduced 
with increased compressibility. 

Another observation worth noting the change the shape the pressure- 
settlement diagrams. With increased initial overburden pressure, the failure mode 
gradually changes local punching failure regardless the relative compres- 
sibility the soil. 

determine the bearing capacity factors, and the shape factors 
determined Beer (3) from small-scale laboratory testing were used, thus: 
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which €,, are shape factors; width footing; and length 
footing. 


obtain approximate value from tests the surface, the 
bearing capacity equation written the form 


which unit weight the soil; the settlement the footing 


TABLE 4.—Bearing Capacity Factors for Strip Footings 


Bear- 
Unit ing 
weight, capac- 
pounds pounds failure, 
per content, per asa 
square per- square square percent- 
inch centage foot inch age 


(1) (2) (3) (4) 


(b) 38° 


boa = 


Solving for yields 


Solving for gives 


7B 


95.6 24.4 10.2 45.5 164 

87.9 18.1 47.7 107 

95.1 152.1 12.7 40.0 23.3 

87.6 121.7 45.2 40.4 16.9 

96.1 266.2 28.1 38.2 13.0 

190.1 45.4 39.0 9.4 

the moment failure; and depth factor defined 

w 
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determine from Eq. Eqs. and are used with assumed factor 
the form 


Introducing the value, determined this way into Eq. approximate 


determine approximate value, from tests with overburden pressure, 
the bearing capacity equation written 


Introducing into Eq. the values previously determined approximate 


TABLE 5.—Bearing Capacity Factors for Circular Footings 


Bear- 
Unit ing 
weight, capac- 
ity, 
pounds Mica pounds pounds failure, 
per content, per per asa 
square square square percent- 
inch centage foot inch age 


(1) (2) (3) (4) 


(b) 38° 


value, for was obtained. Values and were determined 


from and using expressions and for shape factors. 
similar analysis Beer (3), approximate values the bearing capacity 
factors were subsequently used deduce more exact values since improper 
shape factors were assumed deriving first approximate values. However, 
the first values were found sufficiently accurate this analysis. 

The aforementioned analysis employed the secant angle shearing resistance 
The bearing capacity factors were also obtained using the constant tangent 
friction angle 38°, and the results for both cases are summarized Tables 


(a) 
88.3 12.9 9.2 47.3 
95.5 237.9 25.3 38.7 22.9 
87.5 178.4 45.8 39.4 15.6 
95.7 376.6 20.6 37.3 
87.7 337 42.8 37.4 9.4 
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and for the strip and the circular footings, respectively. comparison 
the theoretical values with experimental results shown Fig. 
The theoretical values are given by: 


Equation represents Terzaghi’s proposal (6) use reduced shear strength 
parameters for soils failing local shear. 


500 


Circuler Strip Ftg. 


© O Pure Sand 
Sand with 10%mica 


Tests by De Beer 


Ttong 2 


FIG. Experimental and Theoretical Values 


Test results Fig. indicate that for surface tests incompressible soil 
the values are equal less than the theoretical values (Eq. 11). 
Surface tests compressible soil with 10% mica give value smaller 
than the theoretical one, but larger than that suggested Terzaghi for local 
shear failure. Thus, flat reduction proposed Terzaghi appears 
conservative for surface tests. With increased overburden pressure 
psi (103 the bearing capacity factor, for relatively compressible 
soils approaches that for relatively incompressible soils, and both become closer 
the Terzaghi values for local shear. This consistent with the fact that 
with increased overburden pressure the mode general shear failure disappears, 
and failure the local punching mode, regardless the soil compressibility. 
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The above data show that the effect soil compressibility more pronounced 
surface tests compared tests with initial overburden pressure. For the 
latter with moderate overburden pressure, the Terzaghi proposal use reduced 
strength parameters agreement with the present measurements. 

The values the bearing capacity factor, obtained from surface tests 
are summarized Tables and for strip and circular footings, respectively. 
The values for the strip footing compressible soil are nearly 70% those 
for incompressible soils, increasing 80% the case the circular footing. 
The value for the strip footing incompressible soil compares well 
with the theoretical value based Eq. and the analysis made Caquot 
and Kerisel (1) leading Eq. shown Tables and values were 
not listed for tests with overburden pressure because these tests were used 
only determine 


TABLE 6.—Summary Actual and Predicted Compressibility Factors 


Mica 
pounds content, 
per square asa 
inch percentage degrees 


(1) (2) (3) 


Note: 


far, test results were compared with conventional plasticity theories. Now 
actual compressibility factors are determined and compared with the analytical 
factors proposed Vesic (8). The compressibility factors, are introduced 
the classical equation obtain the bearing capacity compressible soils: 


Vo Ng be Seabee + YBN, 


The analytical values the compressibility factors are determined from 


exp | + 0.6 *) tanod + 


(a) Strip footing 
45.5 651 1.0 0.85 
47.7 284 0.69 0.52 
40.4 0.73 0.47 
39.0 0.72 0.40 
(b) Circular footing 
47.3 243 0.82 0.65 
39.4 0.68 0.70 
37.4 0.93 0.6 
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which the rigidity index defined as: 


This index, appearing solutions the problem expansion cavities 
infinite solid, associated with the assumed elastic-ideal plastic behavior 
soil. compute €,., the rigidity index, calculated, then introduced 
into Eq. 14. calculate the average mean normal stress the expansion 
zone taken the initial mean normal stress depth, B/2, below the 
base the footing, and the modulus, corresponding this mean normal 


700 


NOTE: width of plate= 2 5in 
length of plate= 15 in 


w 


(psi) 


FIG. 10.—Variation Subgrade Reaction Modulus with Overburden Pressure psi 
6.9 


stress determined from triaxial compression test results. The measured Poisson 
ratio 0.3 was used. Examples the computation procedure determine 
the compressibility factors are given Vesic (8). Actual compressibility factors 
were determined dividing the actual bearing capacity factor for the soil 
those incompressible soil having the same secant angle shearing resistance, 
Thus, for the same overburden pressure and footing shape compressibility 
factors were determined from: 


ve qc 
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For the relatively incompressible sand with mica, the above factors are 
equal unity. summary the results shown Table Examination 
Table reveals that agreement resulted between measured and predicted 
compressibility factors for tests with overburden pressure. However, tests 
performed with initial overburden pressure indicate generally much higher 
compressibility factors than predicted Vesic’s (8) method. This should 
expected since this analysis was proposed for surface near surface footings. 

The results indicate that Vesic’s method analysis (8) can successfully 
determine the compressibility factors for shallow foundations compressible 
soils. For footings with substantial overburden pressure, the method underesti- 
mates the compressibility factors. Terzaghi’s method suitable for cases with 
overburden pressure, and fairly conservative for surface near surface 
footings. These model test findings should confirmed full-scale field tests 
different cohesionless soils. 

Effect Compressibility Modulus Subgrade modulus 
subgrade reaction, K,, can determined the tangent the pressure- 
settlement diagrams given Figs. The variation this modulus with 
overburden pressure shown Fig. 10. The effect compressibility 
more pronounced than the effect the bearing capacity, particularly 
tests with initial overburden pressure. This probably because the bearing 
capacity reached large settlements, long after initial localized failures have 
already started the soil mass. 


experimental study evaluated the effect soil compressibility the 
ultimate bearing capacity shallow foundations. Small-scale model tests were 
performed two frictional soils with identical shear strength characteristics, 
but different deformation characteristics, Based test results the 
following conclusions are reached: 


Increased soil compressibility constant strength causes decrease 
bearing capacity and increase settlement. 

The proposal Terzaghi use reduced shear strength parameters for 
compressible soils agreement with test results with initial overburden 
pressure, but, for surface tests, conservative. 

The method proposed Vesic use compressibility factors the bearing 
capacity formula satisfactory for surface tests, but too conservative high 
overburden stress. 

With increased overburden pressure, the failure mode changes local 
punching failure, regardless the soil compressibility. 


Increased soil compressibility substantially reduces the modulus subgrade 
reaction. 
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PROBABILISTIC SOIL EXPLORATION: 


Tien Wu' ASCE and Kinfun 


INTRODUCTION 


Geotechnical engineers are familiar with the uncertainties associated with the 
evaluation soil properties and inaccuracies analytical methods. Casagrande 
explained the concept these risks his classic paper (6). The use probability 
theories assess risk and reliability designs represents rational approach, 
and recent contributions have added the theoretical basis reliability 
assessment. slope stability analysis, studies have considered spatial variation 
soil properties and inaccuracies analysis and measurement soil strength 
(12,16,18,19). However, important note that most failures geotechnical 
construction were initiated failure modes not anticipated during design (14). 
Often these failure modes were related undetected weaknesses the subsoil. 
Thus, one major source uncertainty lies the interpretation data from 
subsoil exploration emphasized Terzaghi (17) many years ago. Probability 
methods are helpful the evaluation this uncertainty. 

major problem the interpretation data from soil exploration the 
classification soils the observed points and the use this information 
draw inferences about soils unobserved points (4,5). Baecher (4) has 
introduced the concept probabilistic map, which associate probabilities 
with the materials that may exist various locations. Such maps would allow 
obtain discrete distributions significant soil properties such strength 
compressibility, which turn may used estimate risk reliability. 

This paper describes the use probability methods interpret information 
obtained from borehole samples different stages soil exploration program. 
addresses the special problem probable failure caused the presence 
weak layer. illustration, data from the embankment failure Cleveland 
(21) used construct hypothetical case history. Opinions several 
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professional engineers are presented subjective probability and compared 
with the results probability analysis. 


Son 


Fig. shows the plan the site, and two representative cross-sections are 
shown Figs. and Preliminary soil exploration the site the embankment 
began with borings WX7 and which were made prior design. Standard 
penetration tests were made these boreholes, and the representative samples 
recovered the penetration spoon were used for classification and determination 
index properties. Limited numbers in. (51 mm) Shelby tube samples 
also were taken from these boreholes. These were used largely for unconfined 
compression tests and consolidation tests. Prior construction, Borings and 
were made from which almost continuous in. (76 mm) Shelby tube samples 
were obtained. Unconfined compression and consolidation tests were performed 


FIG. 1.—Plan Interchange 


most these samples. The site underlain alternating layers lacustrine 
clay and glacial till. The lacustrine clay consists stratified silt and clay. Where 
the silt and clay layers are thin, the material may called varved clay. 
However, thicker layers silt and clay have also been found (21). The data 
obtained from Boring shown Fig. preliminary information for Section 

March, 1971, slide occurred, the extent which shown Fig. 
Large horizontal displacements were also noted section C-C. Additional borings 
were made part the investigation into the failure mechanism. Fig. shows 
the detailed information from all the boreholes section A-A. 

Stability analysis for sections A-A and B-B was made after the failure. For 
safety factor the calculated shear strength along the horizontal portions 
the failure surfaces the order 550 psf (26 For section 
C-C, the calculated shear strength the order 1,100 psf (53 
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our hypothetical case history address the following question: stability 
analysis, made prior construction, had given the information that layer 
with shear strength 550 psf (26 would lead failure section 
A-A and B-B, how could the soil exploration data analyzed access the 
risk due this mode failure? The approach presented the subsequent 
sections allows one evaluate the uncertainties different stages soil 


exploration program and decide any stage whether additional exploration 
would beneficial. 


surface 


*Sy 


FIG. 2.—Section A-A, Preliminary Information 


200 
* Sy — 


FIG. 3.—Section A-A, Detailed Information 


order make the probability analysis, certain information about the 
characteristics the subsoil must known. the present study, the necessary 
information was derived from borehole records and results tests samples 
from the embankment site and nearby site, shown Fig. This information 
analyzed the next two sections. Because have used information available 
much later date, our case history hypothetical and primarily serves 
illustrate how various sources information can utilized planning soil 
exploration. real-life problem some the information normally would 


Final ground surface 
200 - 660 
190 
40% 2 
200 
170 200 
Final ground surfoce 
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come from the engineer’s previous experience with the soil deposit under 
consideration other deposits similar geologic history, both. 


Properties 


many cases, including the present one, preliminary investigation emphasizes 
classification tests and may contain only few strength tests. Hence, 
often necessary estimate the strength from the index properties. large 
number unconfined compression tests have been performed samples 
obtained from this area and from site Index properties these samples 
are also available. 

Linear regression analysis was performed find the correlation between 
the various index properties and the shear strength measured the unconfined 
compression test. The results are shown Table and indicate that the best 


TABLE 1.—Summary Regression 


Source Degrees Sum Mean F-test Significance Correlation 
variation freedom squares squares probability coefficient 
(1) (2) (3) (4) (5) (6) (7) 


(a) Dependent variable = water content 


Model 903.10 903.10 
Error 873.40 10.92 
Corrected total 1,776.50 


Dependent variable liquid limit 


Model 102.63 102.63 5.59 0.0205 
Error 1,469.09 18.36 
Corrected total 1,571.72 


(c) Dependent variable liquidity index 


Model 3.71 3.71 28.46 0.0001 
Error 10.43 0.13 
Corrected total 14.14 


“Ref. provides details analysis variance. 
Note Independent variable unconfined compression strength. 


correlation between the strength and water content. The variations strength 
and water content primarily reflect the different proportions silt and clay 
the stratified material. also important note that all the samples that 
show low shear strength contain mostly clay with few silt layers. Considering 
only the samples clay, there are total samples the shear strengths 
which range from 300-1,360 psf (14-65 with average 900 
psf (43 


STRATIFICATION 


estimate the probable size given layer, the characteristics the 
must known. Necessary characteristics were obtained from 
undisturbed tube samples and block samples taken from the embankment site 
and are described the following paragraphs. 


GT12 PROBABILISTIC SOIL EXPLORATION 1697 


Transition Matrices.—The variations soil stratification can described 
Markov chains (1,7,8,9). The lacustrine deposit the site composed 
silt and clay layers various thicknesses. typical example the stratification 
shown the cross-section Fig. vertical line drawn across such 
cross-section intersects silt and clay. vertical transition matrix used 
describe the change material type along the vertical line. our study, transition 
counts were made every 0.1 in. (2.5 mm), which was chosen the step length 
the Markov chain. The vertical transition matrix obtained from examination 
the block samples and tube samples given Appendix 


FIG. 4.—Cross-section Cleveland Varved Clay 


FIG. 5.—Random Walk Model Soil Layer 


Next, consider the boundaries between silt and clay, shown Fig. 
and be. The undualtions the boundaries modeled horizontal 
transition matrix. this matrix, the changes the elevation the boundary 
with each step, e.g. from Fig. are denoted states, listed 
Appendix Each step the horizontal direction covers distance in. 
(2.5 mm). Block samples the varved clay were taken and the changes 
elevation the boundaries were determined regular intervals. From these, 
the horizontal transition matrix, Appendix was obtained. Details the 


in. 
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computations are given Ref. 20. The horizontal transition matrix then may 
used generate soil profiles. For clay layer with known thickness, 
(Fig. 5), the horizontal transition martix used trace the 
paths random walks beginning and These paths simulate, respectively, 
the upper and lower boundaries the clay layer. The layer terminates 
distance, (Point c), where the two paths meet. The procedures used obtain 
the transition matrices and generate the boundaries are described detail 
(1,2). 

Five hundred simulations each were used obtain the distributions width 
clay layers with initial thicknesses and in. (2.5, 10, mm). 
The distributions are shown Fig. 

Poisson Lines Model.—The Poisson lines model may used simulate 


initial thickness (cm) 
H=5 


Pal 
Oo 
Cc 


100 200 300 
Distance (meters) 


FIG. 6.—Frequency Distribution 500 Simulations 


vertical cross-section soil deposit (10,15). This done generating random 
lines with density Convex polygons are formed these lines. 
Next, the polygons are stretched stretching ratio: 


which the mean width each lithological unit within the deposit; 
and the mean thickness each lithological unit within the deposit. The 
real plane and the stretched planes are shown Figs. 7(a) and respectively. 
Soil types are randomly assigned the polygons based Pi, the fraction 
the soil mass that belongs type This plane now represents vertical 
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cross-section the soil deposit. The geometry the cross-section represented 
the elongation 


which width the cross-section; and height the cross-section. 
line, ab, drawn across the section shown Fig 7(b), would intersect different 

materials sample points spaced regularly distance apart. The frequency 

transition from material material between two sample points denoted 


FIG. 7.—Poisson Lines: (a) Real Plane; (b) Stretched Plane 


n,,. has been shown (10) that the following relations hold between the 
parameters P,, and 


i=l 


2 


G 


which and the prior probabilities that materials and will found 
any location. 

the present case, and n,, for the subsoil the Cleveland site, were 
obtained from the borehole logs. The stratification revealed the block samples 


| — | 
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was used estimate trial and error. Profiles block samples were 
generated, using different value E,, and the vertical transition matrices for 
these profiles were determined. These were compared with the measured vertical 
transition matrix the block samples. means hypothesis testing, the 
writers identified the generated cross-sections the transition matrices which 
were identical with those the block samples. The values for these 
cross-sections are all larger than 400. find that 400 general agrees 
with the distribution width obtained simulation (Fig. which shows that 
E,, the order 1,000 larger for in. (51 cm). 

Since the soft material composed primarily clay (21), assumed 
that their dimensions have the same characteristics those the clay layers 
described the preceding paragraphs. 


return our hypothetical case history. The characteristics the deposit 
described the two preceeding sections are considered information available 
the engineer from past experience. The first task recognize detect 
the soft material when encountered the borings. The correlation between 
strength and water content used for this purpose. order make use 
the aforementioned information, necessary construct model the 
variation strength the clay layers. consider the variations fluctuations 
random field, the probability the mean strength 550 psf (26 
along the horizontal portion the failure surface section extremely 
small. For the present case, are concerned with the probable presence 
soft clay layers the mean strengths which are less than 550 psf (26 
Then, more conservative model would one that contains clay layers having 
different mean strengths. assume that the samples with strengths between 
300 and 850 psf (14 and the mean which 550 psf (26 
represent the soft clay. All other soils are considered the stiff material. 
The assumption arbitrary or, more precisely, inductive (5), and important 
note the subjective nature the model. real problem, this represents 
the engineer’s choice appropriate model for the failure mode based 
his experience and judgment. After the model has been chosen, proceed 
the deductive process which material can classified. 

Decision Criterion.—In order distinguish soft material from stiff material 
according their water contents, decision criterion needed. Criteria frequently 
used decision theory include the maximum likelihood criterion, Neyman-Pear- 
son criterion, probability error criterion, and Bayes risk criterion (11). 
this study, the maximum likelihood criterion was employed because simple 
and does not involve the prior probability the cost associated with decision 
making. 

Suppose have single observation, which may classified belonging 
class Then, the decision statements are d,: classified belonging 

The likelihood function, defined the probability that 
the random observation will take value, given that the state nature 
This probability evaluated from probability density function which 
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depends certain parameters. Furthermore, the likelihood ratio defined 


The maximum likelihood criterion states that A(z) choose d,; A(z) 
choose 

The water contents are taken normally distributed. The mean and standard 
deviation water contents for the soft material are 30.6% and 1.9%, respectively; 
and those for the stiff material are 23.2% and 5.1%, respectively. Then the 
decision statements can written follows: d,: the sample classified 
soft material (1) with mean water content 30.6% and standard deviation 1.9%; 
d,: the sample classified stiff material (2) with mean water content 23.2% 
and standard deviation 5.1%. 


Thus, sample with water content equal classified soft clay 


(1.9) exp 


This gives 27.5% the dividing line between class and class Applying 
this decision criterion the borings drilled before the slide, identified the 
soil elevation 610 (186 Boring and elevation 560 (168 
Borings WX-8, WX-10, and WX-11 the soft material. 


Size Sort Laver 


Having identified the soft layer the boring log and knowing the characteristics 
stratification, the next question is: What are the sizes the layers? Or, 
what fractions the potential failure surfaces sections AA, BB, and 
are composed this material? Here note again that the low strengths are 
attributed the presence thick layers clay. Thus, the size the soft 
layer can treated the size clay layers with thickness exceeding some 
minimal value that will allow low strength measured. Then the Poisson 
lines model may used derive simple proximity rule. 

Proximity Rule.—Proximity rules are based the assumption that materials 
locations close together are more likely the same class than those 
locations far apart. The probability that the material location class 
given that observation distance from the location reveals material 
class can expressed (4,15) 


which prior probability that class material found; parameter 
obtained from Poisson lines model, and 
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which different classes material. For two class problem 


which and are classes. Eq. shows that location next the observation 
point 0), the conditional probability that the material the same class 
tends !.0. location far from the observation point the conditional 
probability tends the prior probability that the material class 


FIG. 8.—Probability Contours for Soft Layer near Elevation 610 (186 Constructed 


from Preliminary Information. Numbers Indicate Probabilities that Soft Layer may 
Present 


FIG. 9.—Probability Contours for Soft Layer near Elevation 610 (186 Constructed 


from Detailed Information Numbers Indicate Probabilities that Soft Layer may 
Present 


Probability Contours.—The proximity rule then used establish lines 
equal probability that certain material may exist. Such probability contour 
map based the assumption that soil layers extend horizontally. Borehole 
locations are considered observation points. other locations, the unobserved 
points, the conditional probability that material may exist depends whether 
found the observation points. Bayes’ theorem, the posterior probability 
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that the material class given (4) 


which prior probability that class material exists; L(k|i) likelihood 


FIG. 10.—Probability Contours for Soft Layer near Elevation 560 (171 Constructed 


from Preliminary Information. Numbers Indicate Probabilities that Soft Layer may 
Present 


FIG. 11.—Probability Contours for Soft Layer near Elevation 650 (186 Constructed 
from Detailed Information. Numbers Indicate Probabilities that Soft Layer may 
Present 


that the material belongs class observation point the material 
class unexplored location; [k|i(d)] conditional probability function 
that material observation point defined Eqs. 9-11; material 
points. The necessary subsoil properties are and which includes the term 


G2 - 1-77. 
07 06 4 
Z\ 
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E,. value 2,000 was chosen conservative estimate after 
consideration the values given the preceding section. The value 
was determined noting the proportions soft clay encountered 
samples from the preliminary borings. 

First, consider the soft layer near elevation 610 (186 sections 
and BB. Boring the only source which provides preliminary information 
this layer. The probability contours based this information are shown 
Fig. The data from borings T-1-71-T-6-71, P-2-71-P-5-71, and S-1-71-S-4-71, 
which were made after the slide, were considered detailed information and 
used update this plot. The result shown Fig. can seen that 
the region defined the 90% probability contour close the actual slide 
area. 

section CC, the soft layer was found near elevation 560 (171 m). Prior 
design, information was available from borings The probability 
contour map based this information shown Fig. 10. Boring which 
was drilled later date before the construction, provided information. 
After the data from borings became available, the probability 
contour map was updated and shown Fig. 

note that the soft material also may present elevations where 


was not detected the boreholes. any point from borehole, that probability 


The calculations are similar those just described and are not repeated here. 
this paper, the writers concentrate the layers detected elevation 610 
and 560 (186 and 170 illustrate the proposed method. 


One major application probability contour map estimate the shear 
strength the soil certain elevation. The results then can used for 
stability analysis. 

The preliminary probability contour map (Fig. used obtain the 
probabilities that the detected soft layer layers might have length 
along the horizontal portion slip surface (Fig. 12). The average shear 
strength along 


which shear strength; the subscripts and refer soft and stiff 
materials, respectively; and distance ab. Stability calculations have shown 
that the critical slip surface does not change appreciably with changes 
Thus, may taken constant and 


the strengths the soft and stiff materials are taken equal their 
mean values 550 psf and 2,500 psf (25 and 120 respectively, and 
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the variation strength within the soft and stiff layers ignored, then the 
probability failure equals the probability that Then, 0.64, according 
the probability contours Fig. include the contributions strength 
variations within layers, necessary consider the joint distribution 
Zand Calculations using Eq. 22, Appendix II, gives 0.67. see that 
the effect strength variations small and may ignored this case. Thus, 
the preliminary information can provide ample warning possible failure. 
similar calculations, the probability, 1,100 psf] along the horizontal portion 
the slip surface section about 0.40, which also suggests high 
failure probability. 

note that these failure probabilities pertain only the specific failure 
surfaces. Since soft clay layers also may exist other elevations, other potential 
failure surfaces should considered for complete assessment failure 
probability. 


The overall objective this soil exploration reduce the uncertainty 


FIG. 12.—Posterior Distributions 


about the strength. Thus examine the values the mean strength along 
and its coefficient variation. 
The mean and variance are, respectively 


which the bar and tilde over s,, s,, and the mean and standard 
deviation, respectively. the strength variation within the soft and stiff layers 


ignored before, then and represents only the error due inadequate 
sampling 


Posterior 
Consultant 
Data 
Engineer 
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which and sample standard deviations; and and number 
samples. 

For Section AA, the computed mean value, (5), and its coefficient 
variation, V,, are summarized Table (Two-class: Preliminary information). 
For comparison, also show the values and 5,, and the coefficients 
variation obtained from and Eq. (Sample data: Preliminary 
information). The large difference the values between Two-class: 
Preliminary information and sample data: Preliminary information indicates the 
important influence the subsoil stratification the evaluation the preliminary 
information. When the detailed information becomes available, for Two-class: 
Detailed information naturally reduced, and approaches that Sample data: 
Detailed information. Then, additional boreholes would not reduce the uncertainty 
below this level. 


PROFESSIONAL JUDGMENT 


Since the identification probable failure modes and the choice subsoil 
models contain large subjective element, important learn the practicing 
engineer’s view this problem. Experienced engineers usually are aware 
the risk posed undetected soft layers even when cannot evaluated 
numerically. safety factor, appropriate for the risk involved, chosen 
the basis subjective opinion which reflects the engineer’s experience and 
professional judgment. 

first step assess the subjective probability, group geotechnical 
engineers was asked choose the shear strength, s,, and factor safety, 
used the design the slope considered this paper. They were 
also asked estimate the probability failure. This was done part 
the Workshop Reliability, ASCE Specialty Conference Analysis and Design 
Geotechnical Engineering, Austin, Tex., 1974. the first questionnaire, only 
the limited information shown Fig. was given the group. well-known 
consultant with years’ experience chose 1,000 psf (48 F,, 
1.5, and estimated 0.01. 

may consider the consultant’s choice and F,, jointly reflect his 
overall uncertainty because likely that engineers not assign quantitative 
values each source uncertainty. Then the safety factor with respect 
the mean strength approximately 


Assuming that the safety factor has normal distribution, the risk 


F,-1 


0.01, obtain the coefficient variation, 0.25. 

The most conservative reply, Engineer who had eight years experience, 
consists 1,000 psf 1.5, and 0.10. Similar calculations 
were made and the results are shown Table subjective probabilities. 
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the second questionnaire, the group was given the detailed information 
(Fig. obtained after extensive soil exploration. The well-known consultant 
chose 1,100 psf and 700 psf, (53 and for the varved clay 
and silty clay, respectively, F,, 1.5, and estimated 0.001. Engineer 
chose 1,060 and 560 psf (51 and for the varved clay and 
silty clay, respectively, F,, 1.3, and estimated 0.10. Calculations similar 
these were made and the results also are listed Table 

note that the design strengths and safety factors the two responses 
are not very different and yet the two responses differ widely the estimated 
P,. However, when the results the detailed exploration were made available, 
the consultant’s was significantly reduced. the same time, Engineer 
chose reduce F,, and design for the same P,. Both indicate the engineers’ 
reduced uncertainty and reflected the smaller values Consultant: 
Detailed information and Engineer Detailed information, Table The 


TABLE 2.—Shear Strength and Coefficients Variation, Section and 


Preliminary Detailed 
information information Posterior 


(1) (2) (3) (4) 
(a) Probability models 


Two-class 


1,000 psf 


(48 (26 
(26 kN/m?) (26 


(120 (120 


Subjective Probabilities 


(48 (43 

(48 (39 


Consultant 


1,000 psf 


1,000 psf 
(48 


0.06 


(48 


relative reductions are both smaller than that obtained the probability 
model, (Two-class: Preliminary information and Two-class: Detailed information). 

can also study engineer’s changing opinions the design the slope 
the use Bayes’ theorem (3). 


which (a) and prior and posterior distributions, respectively; 
paramater the distribution; L(a) likelihood function the data for 
given value and normalization factor. the present case, the 
prior probability represents the engineer’s choice and estimates after was 
given the preliminary information (Consultant: Preliminary information, Table 
2). The detailed information, which the sample statistics the test results 
from all the boreholes the slide area, the same Case 2(b). using 
Eq. and the relationship for conjugate distribution (3), obtained the posterior 
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distributions, the means and variances which are given column Table 
There almost difference between the two posterior distributions which 
plotted Fig. The engineers’ opinions (Consultant: Detailed information 
and Engineer Detailed information, Table also are shown Fig. 12. 
note again that both the consultant’s and engineer X’s subjective opinions 
are more conservative than the posterior distribution. The reductions 
are smaller than that obtained Bayes’ Theorem. 

realized that answers short questionnaire cannot fully express engineers’ 
opinions and judgment. Nevertheless, encouraging that the engineers’ opinions 
qualitatively are agreement with the model estimates. The larger coefficients 
variation represented the engineers’ opinions are not surprising; the engineers 
may have mind other uncertainties besides the failure mode considered 
this analysis. While detailed comparison unwarranted, the limited results suggest 
that quantitative measures subjective probabilities can obtained and 
incorporated probability models. Research into the engineers’ assessment 
various uncertainties should very rewarding, and the results should advance 
the adoption probability concepts practice. 


Summary 


Probability concepts are applied the interpretation soil exploration data 
site where failure weak layer considered probable. The subsoil 
modeled two-class material with soft clay layers included within stiff 
clay. hypothetical case history constructed illustrate the interpretation 
the data obtained different stages the soil exploration program. The 
process consists the following sequence: Detection and recognition the 
soft material from limited information are based the correlation between 
strength and index properties and decision rule based likelihood. the 
next step, which represents inference, the proximity rule used obtain the 
probability that soft layer exists unexplored locations given that has 
been detected not detected given explored locations. This probability then 
used obtain the probability distribution shear strength and estimate 
the probability failure. 

The model two-class material addresses the problem weak layers 
with unknown dimensions. allows the engineer assess the uncertainties 
any stage exploration program and the reduction the uncertainties 
exploration progresses. The methed useful tool but the engineer first 
must identify the failure mode and choose the subsoil model. This inductive 
process based largely the engineer’s judgment. 

The approach used this example general and may applied other 
problems, such those where the presence particular soil layer threatens 
the safety the structure. The following summarizes the procedure assess 
the risk caused the suspected layer: (1) Calculate the dimension dimensions 
the layer that would result safety factor (2) use the proximity 
rule, Eq. calculate the probability that the layer may have the dimension 
dimensions computed step This probability equal the failure 
probability. The parameters and Eqs. may evaluated from borehole 
records. may necessary assume value for value between 100 
and 500 would seem reasonable for sedimentary materials. plan the soil 
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exploration program, the following additional steps should taken; (3) determine 
the safety factor, the controlling soil property such mean shear strength, 
function the dimension dimensions the layer; (4) use the proximity 
tule, Eq. 12, calculate the probability distribution the dimension 
dimensions. This also gives the probability distribution the safety factor 
the soil property; (5) calculate the variance, coefficient variation, the 
safety factor the soil property; (6) repeat steps new data are obtained 
from boreholes. The progressive decrease the variance with increasing data 
measures the reduction uncertainty. Cost effectiveness requires that the cost 
additional exploration worth the reduction uncertainty. 


The writers are grateful the National Science Foundation for the support 
this research under its Grant ENG 76-20852, and Baecher for his 


review and criticism the manuscript. Ali provided valuable assistance 
with the calculations. 


Vertical transition probability matrix: 


clay silt 
clay 
silt 0.781 


States horizontal transition matrix: 


State Change elevation, millimeters 
less than —8.50 
—8.50 —6.25 
—6.00 —3.75 
—1.25 
—1.00 1.25 
1.50 3.75 
4.00 6.25 
6.50 8.75 
more than 8.75 


Horizontal transition probability matrix: 


0.0 0.238 0.286 0.095 0.048 0.095 0.095 0.048 
0.0 0.130 0.261 0.304 0.174 0.087 0.0 0.0 
0.074 0.012 0.099 0.025 0.025 0.0 
0.0 0.038 0.876 0.063 0.013 
0.0 0.078 0.320 0.431 0.065 0.026 0.013 
0.0 0.029 0.086 0.086 0.114 0.314 
0.0 0.0 0.273 
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When and are random variables 


exp 


inwhich 


The function, (7), obtained from the probability contour map and the 
integration carried out numerically. 
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The following symbols are used this paper: 


distance between sampling points; 
decision statements; 
safety factor; 

fraction slip surface that passes through soft material; 
total length slip surface; 


L(z/i) probability that observing given state 


portion slip surface that passes through soft material; 
n,, frequency transition from material material 
probability that material will found; 
mean shear strength; 
standard deviation shear strength; 
coefficient variation; 
width soil layer; 
observation; 
intensity Poisson lines; and 
standard deviation. 
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LATERAL PILE RESPONSE EARTHQUAKES 
Takaaki and Leland Kraft, Members, ASCE 


INTRODUCTION 


Liquefaction and pore pressure build-up cohesionless soils around piles 
during earthquakes cause the natural frequencies soil-pile systems decrease, 
the vibration amplitudes piles increase, and, extreme events, superstruc- 
tures collapse, due large pile movements and the P-A effects (8). Adequate 
seismic design criteria for pile-supported structures liquefiable sand deposits 
are needed design piles that will provide desired level support. Proven 
analysis methods, however, are not readily available predict the seismic 
performance pile-supported structures saturated sands. 

reliable prediction pile behavior saturated sands requires model 
that can follow pore pressure build-up during earthquake-type irregular loadings, 
and rational procedure evaluate seismic soil-pile interactions effects. 

This paper describes nonlinear numerical method for soil-pile-structure 
interaction saturated sands. The method combines liquefaction model (21) 
with nonlinear, seismic p-y relationships. concept derive nonlinear seismic 
p-y relationships based our previous studies the seismic p-y relationships 
for layered soil conditions with equivalent linear soil parameters (19,20). The 
general validity the method was evaluated analyzing shaking table test 
results model soil-pile-structure systems Maruyama (31). Although the 
method needs further confirmation, represents state-of-the-art technology that 
can used for practical problems assess the potential effects soil nonlinearity 
and pore pressure build-up free-field soil response and seismic soil-pile-struc- 
ture interaction. 


Overview 


Soil-pile interaction can studied using foundation spring methods (51), 
beam-on-Winkler foundation methods (1,2,7,20,27,34,37,39,40), elastic solutions 
(18,23,36,50), finite element methods (3,20,25,44). these methods, the 
beam-on-Winkler foundation method most suited sophisticated design 


Consultant, McClelland Engrs., Inc., Houston, Tex. 

Special Projects Group, McClelland Engrs., Inc., Houston, Tex. 

Note.—Discussion open until May 1982. extend the closing date one month, 
written request must filed with the Manager Technical and Professional Publications, 
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analyses, the method economical and can incorporate complex nonlinear 
soil-pile interaction effects. Several researchers have used the beam-on-Winkler 
foundation method study nonlinear aspects seismic soil-pile interaction 
(1,7,34). These studies showed that the cyclic degradation soil-pile stiffness, 
due soil-pile interaction and free-field response, was important the pile 
response. 

Pile response from the beam-on-Winkler foundation method totally depends 
assumed soil-pile springs. Well-established procedures determine the soil-pile 
springs under seismic loading conditions the nonlinear range, however, are 
not present available. Most previous studies nonlinear, seismic soil-pile 
interaction adopted the lateral load-deflection (p-y) relations obtained from 
limited number pile load tests for pile-head loading conditions (42). 

expand our understanding dynamic soil-pile springs, the writers analyti- 
cally studied the fundamental characteristics soil-pile springs under dynamic 
loading conditions (18,19). These studies demonstrated that the soil-pile springs 
are affected the relative stiffness the pile and soil, pile diameter, pile-head 
fixity condition, and excitation frequency. Also, the soil-pile springs for seismic 
loading conditions can differ from those for pile-head loading conditions that 
simulate wave loadings. Most these factors are not adequately accounted 
for the p-y relations currently used. Therefore, need exists establish 
procedure determine rational, nonlinear p-y relations and study nonlinear 
seismic soil-pile interaction. 


Anatysis 


Analytic typical analytic model, based the beam-on-Winkler 
foundation method for this study, shown Fig. pile-supported structure 
idealized single-pile system. The superstructure represented 
series lumped masses, connected beams, that possess lateral degrees-of- 
freedom only. Flexural rigidity the beam right above the pile head may 
adjusted simulate rotational constraint the pile cap. The pile discretized 
into segments, and each pile segment connected free-field soil through 
discrete soil-pile elements spring and dashpot) (Fig. 2). These discrete 
elements represent: (1) Spatial variation soil properties; (2) nonlinear stress- 
strain behavior soils; (3) pore pressure build-up the free field and around 
pile; and (4) three-dimensional soil-pile interaction. 

Seismic excitation the model provided free-field soil displacements. 
Free-field soil response obtained separate analysis nonlinear site 
response method, SRANG, based effective stresses coupled with pore pressure 
buildup (21). 

Linear p-y Relationship.—For pile heterogeneous and linear elastic soil 
layers, Kagawa and Kraft (19) showed that the soil reaction, some depth 
unit pile length can approximated 


p=E,d,y + 2p, B(V, + 


which effective Young’s modulus soil; the average soil-pile 
spring coefficient invariant with depth; the pile displacement relative 
the free field; the mass density soil; the equivalent width 
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the pile; and and the compression and shear wave velocities the 
free-field. The soil properties Eq. vary with depth. The soil-pile element 
that represents Eq. shown Fig. 

The soil-pile spring coefficient, for elastic soil conditions, 
can related the local pile flexibility factor shown 
Fig. Since depends loading conditions, obtained for seismic 
loading are shown Fig. 

for heterogeneous, elastic soil conditions may approximated those 
for homogeneous, elastic soil profiles when the local pile flexibility factor 
evaluated using appropriate average for the site (19). The detailed 
procedures determine for heterogeneous, elastic soil conditions are described 
Kagawa and Kraft (19). 

The velocity-dependent term Eq. represents the three-dimensional wave 
energy dissipation from the pile using dashpots that absorb plane and waves. 

Nonlinear p-y with Pore Pressure Free Field.—Nonlinear Response: 
incorporate the nonlinear stress-strain behavior soils into the p-y relationship 
Eq. used nonlinear stress-strain relationship, schematically shown 
Fig. The backbone curve defined 


which shear stress; shear modulus low shear strain; 
shear strain; and reference shear strain The maximum shear 
stress, sin sin derived considering the passive failure 
the soil element adjacent pile. the effective overburden pressure 
and the internal friction angle. The function, may determined 
conventional cyclic triaxial tests simulating the stress system around pile and 
can modeled by, for example, the hyperbolic model (10,24), the Ramberg-Os- 
good model (38,41,43), and the Martin-Davidenkov model (28). Unloading and 
reloading curves may constructed using the Masing rule (32). High shear 
strain concentration occurs the vicinity pile due soil-pile interaction, 
and the soil strain used Eq. for the evaluation soil stiffness around 
pile must representative average this variation. The average shear 


which Poisson’s ratio; and factor that depends the shape 
loaded area, distribution soil modulus, and relative stiffness pile 
soil, among other factors. estimate the factor, for typical offshore flexible 
piles, used closed-form solution for pile homogeneous, elastic soil 
layer (18). Our analysis showed that factor typically ranged from 2.0-3.0. 
This range values agreement with those used Matlock (33) for 
the lateral response pile under wave loading and Skempton (48) for footing 
settlements. 

Eq. can used with Poisson’s ratio compute average modulus, 
E,. This modulus then can combined with Eq. generate the dynamic 
p-y relationship nonlinear soils, which 
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This p-y relationship may vary along the pile length account for vertical 
variations the soil properties. 

Pore Pressure Generation Model.—Soil-pile interaction governed the 
effective stress state the soil the immediate vicinity the pile. The change 
effective stress the soil around the pile can differ from the change 
the free-field. The change total mean stress due the lateral pressure, 
and the change pore pressure, result the change the effective mean 
stress, AG,,. The change the normalized effective stress 


the induced excess pore pressure; and the change total mean 
stress. The backbone curve Eq. can modified, shown Eq. 
account for the change effective stress: 


The change relative considered small and neglected. The 
degraded soil modulus Eq. then may used Eq. account for 
the soil nonlinearity effects the p-y relationship. 

Since the occurrence the Niigata and Alaska earthquakes, considerable 
effort has been devoted clarification the fundamental mechanism 
liquefaction (4,5,12,13,14,15,26,29,35,45,46,47,54,55,56). Although models are 
available predict pore pressure response, the methods generally involve many 
parameters from various types laboratory tests, and the basic rules evaluate 
the pore pressure response against earthquake-type irregular loadings are not 
clearly defined some models. 

The model used here was simple and computationally efficient (21). According 
the model, the pore pressure increase associated with shear stress increase 


(t2/8,) 7 
(+, 
which pore pressure increase; shear stress; and the potential 


function for pore pressure build-up that determined for given density 
soil and confining pressure. given 


which and model parameters; n(N) the function’’ that 
affected only the number times the soil sheared; and 
function representing the change the pore pressure generation capability 
the soil caused increase the pore pressure ratio, u/G,. 
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All model parameters can determined from conventional, cyclic, undrained 
simple shear triaxial tests. The model requires only limited number 
laboratory tests and only few model parameters. Eq. can numerically 
integrated obtain pore pressure increases for irregular stress variations. Details 
the model are presented elsewhere (21), but typical comparison measured 
and predicted pore pressure responses from the model shown Fig. 

the field, drainage occurs certain extent, and the pore pressure 
redistributes with time. Seed, Martin, and Lysmer (45) studied the liquefaction 
characteristics idealized uniform sand deposit determine the influence 
drainage the number cycles liquefaction. The study indicated that 
pore pressure dissipation had small effects the pore pressure response for 
medium fine sand deposits with the coefficient permeability less than 
about 1.5 cm/s. Similar conclusions may derived from the study 
Finn, Lee, and Martin (5). For soil conditions with permeabilities higher 
than this value, the pore pressure build-up appreciably affected the dissipation 
effects, and the proposed method fails. Under these conditions, however, the 
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FIG. 5.—Measured and Predicted Pore Pressure for Niigata Sand and Irregular Loading, 
Niigata Earthquake W-E Component 


pore pressure generated the free-field and due the soil-pile interaction 
dissipates rather quickly, and the impact pore pressure build-up pile response 
expected small. The effects pore pressure dissipation around pile 
will assessed later section, this dissipation can important the 
pile response. 


With Eq. the relationship between and reduces 


which represents transformation between the p,-y relationship and the stress- 
strain relationship the soil. The soil-pile spring thus derived varies according 
the nonlinear stress-strain properties the soil. 

Pore Pressure Build-up due Soil-Pile Interaction.—During seismic event, 
the soils the vicinity pile may develop pore pressure excess that 
the free-field, due solely the soil-pile interaction. The pore pressure gradient 


TIME, SECONDS 


1720 DECEMBER 1981 GTi2 


generated the soil-pile interaction may high, and rapid redistribution 
the pore pressure may take place when the loading rate slow and the permeability 
the soil high. Under these conditions, the pile response mainly affected 
the degradation soil-pile springs due the pore pressure build-up 
the free field. 

Martin, Lam, and Tsai (30) presented the results radial dissipation pore 
pressure idealized soil-pile system. According their results, the rate 
dissipation depends the pile diameter, the permeability and compressibility 
sand, and the loading rate. For typical values the coefficient compressibil- 
ity (4.2 and permeability cm/s) for well-drained 
sand, high pore pressure gradients around pile still can maintained for 
large diameter piles loaded period ranging from 3-7 sec. Thus, the effects 
the pore pressure, developed due the soil-pile interaction pile response, 
cannot neglected since they have the potential impact pile response. The 
pore pressure generation due the soil-pile interaction was obtained our 
analytic model assuming that the representative soil shear stress around 
pile for the pore pressure assessment can approximated p,/(2 
Eq. 

The effects the assumption pore pressure generation due soil-pile 
interaction, and the rate dissipation pile response, were evaluated 
studying two cases: (1) The soil-pile springs degraded only due the pore 
pressure build-up the free-field; and (2) the soil-pile springs degraded due 
the pore pressure build-up caused the soil-pile interaction. 

Analysis Procedure.—Seismic response analysis soil-pile-structure system 
the proposed method consists two steps: (1) Nonlinear free-field analysis; 
and (2) nonlinear soil-pile-structure analysis. Although any nonlinear site response 
method based effective stresses coupled with pore pressure response may 
used accomplish the first task, nonlinear site response method SRANG 
was developed. Nonlinear stress-strain relations were approximated the 
hyperbolic model, and the pore pressure response was obtained using the model 
presented the writers (21). 

With the computed free-field soil motion, nonlinear soil-pile-structure 
interaction analysis carried out the time domain. each time step, free-field 
soil response read from tape, and equations motion soil-pile-structure 
system are established assuming the appropriate soil-pile spring coefficient, 
Since depends the deflected shape the pile and the relative stiffness 
between the soil and pile, iteration loop provided select appropriate 
value for each time step. This iteration usually, however, not required 
varies slowly with time and the time interval small enough follow 
this variation. After the equations motion are solved, system properties such 
soil constants, pore pressure around the pile, and the soil-pile spring coefficient 
are updated for the next time step. The procedure was automated into computer 
program NONSPS. Time-step integration for the programs SRANG and NONSPS 
was performed the unconditionally stable a-method (11). 


The proposed nonlinear seismic p-y relationships involve several simplifying 
assumptions include the effect soil nonlinearity and pore pressure build-up. 
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basic assumption related the scaling relationship between p-y relations 
and stress-strain behavior soil the concept average soil strain around 
pile. Redistribution with time the pore pressure generated the soil-pile 
interaction not incorporated into the relationships. The effects these 
assumptions pile response will examined following. 

Performance and validity the proposed method should tested using 
field observational data seismic response pile-supported structures. Avail- 


TABLE 1.—p-y Elastic and Plastic Ranges 


kilograms per square foot thousands pounds per foot 
Depth/B Reference Equation Reference Equation 
(1) (2) (3) 


DEPTH/B 3.0 

s 39° 
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REESE, COX, AND 
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FIG. 6.—Predicted and Measured p-y Relations 


able field data, however, are limited seismic events low intensity where 
the responses pile-supported structures essentially were within the linear 
Thus, field data are not available examine the performance nonlinear 
seismic soil-pile-structure interaction methods. The proposed nonlinear p-y 
relations were first compared with the static p-y relations Reese, Cox, and 
Koop (42), and the assumptions and performance the seismic soil-pile 
interaction method presented herein were evaluated comparing computed 
results with tests model soil-pile-structure system Maruyama (31). 
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Evaluation Nonlinear p-y.—Nonlinear p-y relationships, using Eq. were 
compared with those Reese, Cox, and Koop (42). make the comparison, 
the writers considered the test conditions the Mustang Island (42): Earth 
pressure coefficient rest 0.4; angle interal friction 39°; coefficient 
subgrade reaction 125 pci (34,328 effective unit soil weight 
pef (10.5 and pile diameter ft. The shear stress-strain relationships 
the sand the depth the p-y evaluation were based maximum shear 
modulus, G,,,,, obtained from (10) 


2¢ 0.5 
The void ratio 0.66, used Eq. 10, was estimated from the measured water 
content the site. Using these parameters, constructed p-y relationships. 
was 2.5 and was 1.2, which are typical values for piles layered soil 
conditions, with soil modulus parabolically increasing with depth, loaded 
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FIG. 7.—Soil-Pile-Structure Model for Shaking Table Tests, after Maruyama (31) 


the pile head. The computed slope the p-y curves the origin, and the 
ultimate values from the two methods, are compared Table The results 
are reasonable agreement. 

generate the values function the shear modulus reduction 
curve Eq. was represented the simple hyperbolic model 


max 


The p-y relation Eq. compared with that Reese, Cox, and Koop 
(42) Fig. for depth pile diameter ratio Our p-y relation agrees 
reasonably well with that Reese, Cox, and Koop (42) for small and large 
y/B values. Since the agreement for intermediate values depends the 
characteristics the soil model used and easily can improved using 
soil model that involves more parameters than the hyperbolic model Eq. 
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11, can conclude that Eq. may used predict nonlinear p-y relationships. 
Further refinement the results was not attempted here more studies are 
needed compare the performance the proposed method the light 
field observational data. For this reason, the present study based the 
simple hyperbolic model Eq. 11. Although the preceding comparisons are 
limited the p-y relationship for pile-head loading conditions, the proposed 
p-y should also perform well for seismic loading conditions. 

Shaking Table Tests Model Soil-Pile Structure System.—The shaking table 
tests, which were carried out the Ministry Construction Japan, were 
reported Maruyama (31) and Hakuno, Iwasaki, and Tatsuoka (9). The steel 
piles the model were 726 long, wide and had rectangular cross 
section. Two different thicknesses and 4.5 mm) model piles were used. 
Concentrated masses 4.59 kgf for the pile and 0.62 kgf for the 4.5 
pile were attached the top the model piles. 

The piles were placed beds pluvial deposited Toyoura sand placed 
relative densities between and 50% 1,265 495 799 steel 
box supported shaking table. eliminate the constraining effects the 
side walls, they were designed rotate freely with soil (Fig. 7). 

Acceleration, pore pressure, and dynamic earth pressures were measured 
the tests. Instrumentation lay-out for the tests shown Fig. Acceleration 
pick-ups were equipped with light-weight stabilizers maintain their positions 
during shaking. 

Characteristics Test Sand.—The Toyoura sand had uniformity constant 
1.46, maximum void ratio 0.96, minimum void ratio 0.64, and 
permeability coefficient 0.0329 cm/sec. The duration shaking was short 
and the pore pressure gradient was low for the small model. Our analysis, 
based Darcy’s flow theory, indicated that the pore pressure movements 
were small during the tests and that measured pore pressure responses were 
not affected pore pressure redistributions most test results. The hyperbolic 
stress-strain model was used fit the estimated stress-strain behavior the 
test sand. 

Liquefaction characteristics Toyoura sand were provided Tatsuoka (52). 
The tests were performed confining pressures 0.5 and 1.0 kgf/cm’, 50% 
relative density, and isotropic consolidation. The sand layers analyzed this 
study had average overburden pressure 0.04 and 31% relative 
densities, and anisotropic consolidation. The test results Tatsuoka (52) were 
modified for the difference confining pressures, relative densities, and 
consolidation conditions determine the liquefaction characteristics the sand 
layers the model tests. 

Analysis Model Test were made predict the response 
the model piles and compare the predicted and measured responses. Although 
the real sand layer was three-dimensional, our analysis was made two 
dimensions, and the friction between the soil and the walls parallel the shaking 
direction was ignored. The rotating wall was simulated series beam 
elements. The writers used the nonlinear plane-strain program NONPLS 
compute the pore pressure response the sand layer without the presence 
the piles. Fig. compares the computed and measured pore pressure responses 


for typical case. Predicted pore pressure responses agree very well with test 
results. 
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Behavior model piles saturated sand layer was then predicted using 
the free-field responses the sand layer input for the NONSPS analysis. 
The measured fundamental resonance frequency the model piles was about 
before shaking. The soil-pile springs softened, and the measured funda- 
mental resonance frequency the soil-pile-structure system decreased with 
continued shaking about Hz, due the pore pressure build-up supporting 
soils. The frequency the sinusoidal table input motion was Hz. Therefore, 
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FIG. 8.—Measured and Predicted Pore Pressure Responses Sand Layer 


the model piles passed through the resonant conditions, and the soil-pile-structure 
system exhibited relatively large responses the transient stage. The measured 
and predicted pile responses Fig. support this phenomenon. Our prediction 
good agreement with the measured data about three cycles which 
the pile response controlled the elastic p-y. This indicates that our p-y 
relations small values are reasonable. Between three seven cycles, however, 
the predicted pile response underestimated the measured results. This difference 
consistent with that observed the computed and measured resonance behavior 
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the sand layer about Hz, and may due the possible existence 
high frequency components actual input table motions. After cycles 
loading, the upper one third the sand layer liquefied and piles started 
move with small constraint soil motion. The predicted peak accelerations 
vary erratically after the sand layer starts liquefy. The mean values the 
histories, however, agree with the observed data. The fluctuation may due 
the lack viscous damping the pile provided the liquefied sands 
our analysis, and the comparison should improved incorporating this 
effect into our study. 

The predicted model pile responses include the effects pore pressure build-up 
the free-field sand layer, and due soil-pile interaction. the upper portion 
the free field was brought into liquefaction few cycles loading, the 
pore pressure generated the soil-pile interaction had only small effect 
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FIG. 9.—Measured and Predicted Pile Response Liquefying Sand Layer 


the pile responses. This was confirmed separate analysis that included 
only the pore pressure build-up the free field. Thus, the validity the pore 
pressure model the proposed soil-pile springs could not assessed using 
the model test data. Pore pressure build-up the soil around pile, due 
the soil-pile interaction, however can have large effect pile response, 
demonstrated following. Further studies are needed quantify this aspect 
soil-pile interaction. 

Significance radiation damping pile response was evaluated for cases 
with and without the damping term Eq. The radiation damping affected 
the pile response only early stage loading. 


The seismic response typical 400 (122 water depth offshore structure 
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was computed with the proposed method demonstrate the impact pile 
response pore pressure build-up the soil around pile. The structure 
was represented single vertical pile system the properties which are 
given Kagawa and Kraft (20). The soils the site consist uniform 
sand depth 400 (122 m). Three soil densities (D, 55%, 65%, and 75%) 
were used the analysis. Shear moduli the sand site were computed using 
the formula Iwasaki, Tatsuoka, and Takagi (17). Liquefaction strengths 
the sand used here are shown Fig. 10. The soil motion depth 400 
(122 was based the deconvoluted Centro N-S component 200 
(61 sand site with 85% relative density. The maximum acceleration 
the base motion was 0.173 and the duration shaking was 10.24 sec. 

The writers first performed nonlinear free-field analysis using the program 
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FIG. 10.—Liquefaction Strengths for Example Analysis 
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FIG. 11.—Response Hypothetical Off-shore Sand Site 0.305 


SRANG. Fig. represents the acceleration and displacement response spectra 
obtained from this analysis for the site with 55%. Large pore pressures 
developed only the upper 15-20 (4.58-6.1 m). pointed out other 
investigators (6,28,49), nonlinear site response methods predict ample high 
frequency responses and large low frequency displacement responses, 
compared with the equivalent linear method SHAKE. The results from SHAKE 
analysis are also shown Fig. 11. The difference the high frequency response 
characteristics the motions from the two methods has relatively small impact 
the response the structure studied here with natural frequency 0.41 
Hz. The difference the displacement spectra low frequencies, however, 
can have large effect the response tall, flexible structures such 
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offshore pile-supported structures deep water. 

The free-field motion obtained for the sand site was then used input 
NONSPS analysis. Two separate NONSPS analyses were made, the first 
only considering the effect pore pressure build-up the free field, and 
the second also including the effect pore pressure build-up around the pile 
due soil-pile interaction. 

the pore pressure develops around the pile due the soil-pile interaction, 
liquefaction may extend about pile-diameter depth, and large moment 
and shear can transmitted pile-diameter depth. the other 
hand, shear and moment profiles for the three soil densities were similar 
the soil-pile springs only degraded due the pore pressure build-up the 
free-field. Thus, our analysis indicated that the pile response was affected more 
the pore pressure build-up around the pile than that the free field. 
Fig. shows typical comparison for 65% relative density between the 
shear and moment profiles for these two conditions. The pile shear normalized 
the weight the superstructure and the pile moment divided the 
weight the superstructure and the pile moment divided the weight 
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FIG. 12.—Effects Pore Pressure around Pile Pile Shear and Moment 


times the effective height the superstructure. The effective height was 
determined using the first mode vibration shape the structure. The maximum 
shear and moment decrease rather rapidly with depth when the soil-pile springs 
degraded due solely the pore pressure the free-field. the pore pressure 
develops around the pile due the soil-pile interaction, the maximum shear 
along the pile does not occur the same time, due the progressive softening 
the soil-pile springs. Therefore, the variation the maximum pile shear 
somewhat irregular. The same phenomenon observed for the maximum 
pile moment. 


The numerical method presented this paper can used assess seismic 
response characteristics off-shore well onshore pile-supported structures 
sited saturated sand deposits. The primary feature the proposed method 
consistent approach determine nonlinear seismic p-y relationships. The 
pore pressure build-up effects the free-field around pile are included 
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the proposed p-y relationships. The proposed method reproduced the static p-y 
relations under the pile-head loading condition Reese, Cox, and Koop (42) 
with reasonable degree agreement. 

Performance and general validity the method were examined through 
analysis shaking table tests model soil-pile-structure systems. The pore 
pressure generation model the proposed p-y could not fully tested from 
this analysis, due the lack pore pressure measurement data the vicinity 
the model piles. General agreement between the predicted and measured 
pile responses during the liquefaction process, however, encouraging. The 
proposed method provides state-of-the-art technique assess the potential 
impact pore pressure build-up around pile pile response. 

Using the proposed method, the writers analyzed seismic response characteris- 
tics idealized, off-shore, pile-supported structure sand site. Results 
indicated that the pile response can significantly affected pore pressure 
develops around pile due the soil-pile interaction. Shear and moment within 
pile dramatically increased, due the pore pressure build-up around the 
pile. The study identified the potential importance the pore-pressure build-up 
around the pile, due the soil-pile interaction, for evaluating the lateral pile 
responses during earthquakes. 

More studies, however, are needed clarify the pore pressure response 
characteristics sands around pile under dynamic loading. Large stress 
concentrations and high pore pressure gradients around pile may invalidate 
use the pore pressure model that based the simple shear condition. 
Pore pressure dissipation the radial well vertical directions along the 
pile can delay the pore pressure build-up, due the soil-pile interaction. The 
proposed p-y may overpredict the pore pressure build-up around the pile. Pore 
pressure measurements around full-scale large diameter piles under dynamic 
loads are needed clarify these points. With such field data available, the 
validity the proposed method and the significance the pore pressure build-up 
due the soil-pile interaction pile response can fully examined. 
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The following symbols are used this paper: 


model parameter; 

effective width pile; 

void ratio; 

Young’s modulus soil; 

flexural rigidity pile; 

pore pressure function; 

shear modulus reduction function; 

scaling factor between pile displacement and soil strain; 
shear modulus low strain; 

earth pressure coefficient rest; 

local pile flexibility factor 
memory function; 

soil reaction unit pile length; 

pore pressure ratio (Au 
radius pile; 

compression and shear wave velocities soil; 
pile displacement relative free field; 

model parameter; 

shear strain and reference shear strain; 

average soil-pile spring coefficient; 

Poisson’s ratio soil; 

mass density soil; 

effective mean principal stress; 

effective vertical stress and initial effective vertical stress; 
shear stress; 

potential function for pore pressure generation; 
angle internal friction; 

pore pressure increase; 

change total mean stress; and 

change effective mean stress. 


A = 
B = 
e = 
f = 
n(N) 
=> 
r => 
y = 
B = 
yandy, 
v => 
Cn = 
G,ando,, 
T = 
= 
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TECHNICAL NOTES 


SEISMIC DISPLACEMENT ANALYSIS EARTH 


the limit equilibrium method seismic design earth dams, the computation 
displacements during earthquakes play important part. Since Newmark 
(1) presented his technique rigid block sliding plane surface, designers 
have tried use compute the displacement the crest during earthquake. 
The sliding block technique gives displacements along fictitious plane which 
not properly defined. The definition given Sarma (3) also arbitrary. 

There are only two kinds surfaces along which rigid body movement 
possible. These are plane and circular arc. All other surfaces, such 
log spiral nonlinear one, cannot accommodate the rigid body mechanics 
sliding. the case plane surface, the sliding purely translational; 
the case circular arc, the sliding purely rotational about the center 
the circle. rigid body movement applied, then for all other surfaces 
equivalent plane circular arc may found. 

Most limit equilibrium methods stability analysis use circular arc slip surfaces, 
which provide simple computational process search for the critical surface. 
Therefore, displacements computed the basis rotation about the center 
the circle are readily applicable. 

Let assume rigid body top circular arc slip surface ABC, shown 
Fig. which being acted upon earthquake load which 
function time. also assume that sometime the past slide has 
been initiated, i.e. that was greater than which the critical 
acceleration. Let the center gravity the body ABCA. the 
center the circle with radius The forces acting the body are its weight, 
the earthquake load, KW, the normal forces, and the shear forces, 
the instant the center gravity moves through angle from 
which the position the center gravity the beginning 
the movement. 

can write the equation motion the sliding body 


which the mass moment inertia; and the distance from 


Civ. Engrg. Dept., Imperial College, London, England, SW7. 

Note.—Discussion open until May 1982. extend the closing date one month, 
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ASCE. Manuscript was submitted for review for possible publication April 15, 1981. 
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this position the sliding body, can find value which satisfies 


also can find value which was operative the beginning 
the slip position G,, which gives 


The value the critical acceleration, computed from the stability analysis 


FIG. 1.—Rotational Sliding Surface 


slopes, Sarma (2). the shear stress the slip surface the beginning 
the slip. 


Assumption Case 1.—We assume that movement, small that 
which case Eq. may written 


Solution Eq. closed form, when function time, not possible 
and numerical solutions may devised. However, special case K(t) may 
studied. 

K(t) represents rectangular pulse that 
for 


for 


The initial conditions are that 
fe) 
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then, for small Eq. may solved with initial conditions given Eq. 
obtain 


K., 
K., 


co 


The displacement, x,,, along the arc the circle is, therefore 


K,, 


2 
The solution Eq. the same form was obtained for sliding 
block plane surface (3), which 


inwhich 

The value for simple geometry shown Fig. This shows 
that approximately equal sin and, for all practical purposes, will 
have value between 0.5 and 

Looking back into the equation motion, Eq. note that the same 
solution would have been arrived the assumption were made that for all 
practical purposes 


With such assumption, statement may made that displacements 
computed the basis sliding block plane horizontal surface are 
approximately equal the displacement along the arc the circular sliding 
surface when the critical acceleration remains unchanged. 

Assumption Case 2.—We assume that during sliding the shear stresses 
the base the sliding surface remain constant and are equal those the 
beginning the slide. this case, Eq. will read 


sin 
sin 
Therefore, Eq. becomes 
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Solution Eq. with initial conditions given Eq. and for rectangular 
pulse, defined Eq. yields 


for 
Eq. shows that there are two possibilities. the pulse duration, 


which occurs during the pulse. the duration, small compared with 
Eq. 


= —K., + v + 2K,,9 ro 


C=(Rd/r* sing, 


3 
% sin’y 
¥-sinycos® 


f = radius of gyration 


FIG. 2.—Variations Quantity, Rd/r’, for Simple Geometry 


case 


case 2,S=S 


To ne, 


FIG. 3.—Angular Displacement for Rotational Sliding Against Duration Pulse 


Keo 
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This analysis shows that for rotating along circular arc, the shear strength 
along the slip surface remains constant, there maximum rotation 
2(K,, which cannot exceeded. But this general large 
and therefore not practical use. For example, for (K,, 0.1, 
and for radius curve, the maximum possible displacement along the 
arc x,, which very large. 

However, this also shows that even when the static factor safety equal 
one, i.e. zero, earthquake acceleration 0.1 will give maximum 
displacement instead infinity, which will obtained the case 
plane sliding surfaces with circular surfaces when the critical acceleration 
remains constant. Fig. shows comparison the displacements obtained 
the two cases. This shows that for small angular displacements, the two 
cases give identical results. The case when the critical acceleration actually 
deteriorates with displacement not included this analysis. 


Newmark, M., Earthquakes Dams and Geotechnique, 
Vol. 15, No. 1965, pp. 139-160. 

Sarma, K., Analysis Embankments and Geotechnique, Vol. 
23, No. 1973, pp. 423-433. 

Sarma, K., Stability Earth Dams and Geotechnique, 
Vol. 25, No. 1975, pp. 743-761. 
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DISCUSSION 


encouraging see the growth and application nonlinear analysis 
seismic response studies and the authors are commended for their interesting 
and informative paper. 

The question low strain damping problem inherent hysteretically 
damped nonlinear soil models, since very small strains they predict elastic 
behavior and attenuation motion. For this reason, prediction residual 
Strains soil masses difficult when using such models. The inclusion 
viscous damping linear model under steady state cyclic loading results 
the maximum stress and maximum strain being nonconcurrent, which 
variance with experimental data. The inclusion viscous damping Iwan 
model will result widening the hysteresis loops the midstrain range 
when harmonically cycled fixed amplitude. 

Experimental work the University Auckland has been carried out 
measure soil damping strains down (15). free vibration torsion 
test apparatus employed using almost friction free air bearing (16,17). 
The data recorded and analyzed using on-line computer for high resolution 
low strain data and accurate post processing. This work established that 
energy dissipation still occurring strains this magnitude. Equivalent 
viscous damping factors about have typically beeen measured. Fig. 
shows some the data gathered from dynamic torsion test. important 
note that this work does not prove the existence viscous damping 
low strains, but establishes that energy absorbed during cyclic loading 
very low strains (one microstrain). 

One reason for using viscous damping reduce the noise,”’ 
the form high frequency oscillations. has been the writer’s experience 
that such high frequency motion often not problem. This numerical noise 
usually occurs outside the frequency band interest cases involving high 
energy seismic events. This likely the case for soft sites, since soil-structure 
interaction will increase the fundamental period the structure being studied. 

Degradation additional sophistication introduced simulate the reduction 
soil stiffness with strain history. study degradation cohesive soils 
was undertaken Bacchus (18), using cyclic triaxial test equipment. Similar 
results were obtained and linear relationship was found between log and 


log Using similar curve fitting procedures the results for normally consoli- 
dated clay were expressed 


1980, Chan-Feng Tsai, Ignatius Lam, and Geoffrey Martin (Proc. 
Paper 15708). 
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1744 DECEMBER 1981 GT12 


which the octahedral stress, initial value and function 
the strain amplitude. 

The study used series tests with variation strain amplitude and 
measured excess pore pressures and the change rheologic properties. 

has been the writers experience that degradation not great importance 
seismic response studies, except very soft soils the near source region 
high magnitude earthquake. The numbers cycles loading reaching 
high strains, say 0.5%, usually limited. Degradation possibly important 
feature severe storm loading off-shore structures where the number 
cycles different order magnitude. 

The transmitting boundary described Joyner and Chen (3) significant 
advance nonlinear methods. The writer has obtained similar results using 
this procedure, but note should made that the methodology based 
the propagation plane shear waves through elastic medium. assumes 
constant shear wave velocity for the material below the transmitting boundary. 


Sreer Strain 


FIG. 14.—Variation Damping with Shear Strain 


This not really the case most seismic events engineering significance, 
thus the procedure approximation most cases. 


15. Mead, L., Soil Report No. 194, School Engineering, 
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Taylor, W., and Parton, M., Testing Soils,’’ Proceedings the 
International Conference Soil Mechanics and Foundation Engineering, Vol. 
1973, 425. 
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Discussion Kurt and Herbert Klapperich* 


The suggestion the paper appears helpful and applicable for practical 
purposes with respect suitable position grouted zone around tunnels 
and subways. For calculations with Finite-Element-Method the constitutive laws 
and their parameters are essential importance with considerable influence, 
e.g., the deformations (1). 

This holds for soils and injections whereby the injections are always time 
dependent two ways: (1) Dependency the deformation the onset 
loading; and (2) the well-known creep. Both influences should taken into 
consideration calculating the soil grouting. this way the factor mentioned 
Table can increased well decreased the case late respectively 
early onset loading after the injection (2). 

The creep deformations are also dependent the onset loading and are 
the same magnitude the initial deformations after excavation. For the 
design procedure the calculation creep deformation important especially 
before tunneling performed and should considered more strongly. 


the paper, the figures for design procedure were only determined 
linear-elastic FE-calculations. interest know whether the stiffness 
the onset loading decrease stiffness was used order consider 
the subsequent creep deformations directly. 

For both assumptions would interest for practical purposes 
far the ratio deformations after certain duration loading (even 
one day) differ from the determined ones. 


16. Miiller-Kirchenbauer, H., 
Boden unter Dauerlast,’’ Baugrundtagung, Berlin, West Germany, 1978. 

17. Richter, T., der FlieBregel bei der Berechnung eines Streifenfundamentes 
auf rolligem Boden mit der Methode der finiten Die Bautechnik, No. 
Berlin, West Germany, 1978. 


“September, 1980, Tan and Clough (Proc. Paper 15716). 

Grundbauinstitut, Civ. Engrg. Dept., Technical Univ. Berlin, Berlin, 
West Germany. 

“Research Asst., Grundbauinstitut, Civ. Engrg. Dept., Tech. Univ. Berlin, Berlin, 
West Germany. 
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Anchors, Boris S. Browzin, GTS:May 81:629:16275 


Collapsible Soils 

Construction of Large Canal on Collapsing Soils, Paul 
C. Knodel, GT1:Jan 81:79:15992 

Design Considerations for Collapsible Soils, Samuel 
P. Clemence and Albert O. Finbarr, GT3:Mar 81: 
305:16106 


Collars 

Deflection of Monopiles with Collar Caissons, 
Massimo Magni and Alex P. Michalopoulos, TN: 
GTS:May 81:667:16221 


Compaction 

Densification of Loose Deposits by Pounding, Robert 
G. Lukas, GT4:Apr 80:435:15376 

Shear Strength of Rockfill, Nick Barton and Bjorn 
Kjaernsli, GT7:July 81:873:16374 


Compaction (Soils) 

Cellular Cofferdam for Trident Drydock: 
Performance, Max D. Sorota, Edward B. Kinner 
and Mark X. Haley, GT12:Dec 81:1657:16733 


Complexity 

Nitrogen-Cycle Model for Aquatic Systems: Analysis, 
Tavit O. Najarian and Jay L. Taft, GT12:Dec 81: 
1141:16703 


Compressibility (Soils) 

Compressibility and Bearing Capacity, Nabil F. 
Ismael and Aleksandar S. Vesic, GT12:Dec 81: 
1677:16718 


Compressible Soils 

Stresses in Soil around Vertical Compressible Piles, 
K. S. Sankaran, N. R. Krishnaswamy and B. K. 
Sharas Chandra, TN: GT1:Jan 81:107:15956 


Compression Waves 

Drilled-Shaft Integrity by Wave Propagation Method, 
Thomas M. Hearne, Jr., Kenneth H. Stokoe, II and 
Lyndon C. Reese, GT10:Oct 81:1327:16582 


Computer Analysis 

Coupled Heat and Water Flows Around Buried 
Cables, Omar N. Abdel-Hadi and James K. 
Mitchell, GT11:Nov 81:1461:16634 

Laterally Loaded Pile Design, Robert L. Sogge, GT9: 
Sept 81:1179:16510 


Computerized Simulation 
Computer Simulation of Creep of Clay, Roland Pusch 
and Paul Feltham, GT1:Jan 81:95:15998 


Computer Programs 

Nonlinear Seismic Response of Soft Clay Sites, Ram 
D. Singh, Ricardo Debry, Earl H. Doyle and Izzat 
Idriss, GT9:Sept 

Random Surface Generation in Stability Analysis, 
Ronald A. Siegel, William D. Kovacs and Charles 
Lovell, TN: GT7:July 81:996:16360 


Computers 
Laterally Loaded Pile Design, Robert L. Sogge, GT9: 
Sept 81:1179:16510 


Cone Penetrometers 

Cone Penetration in Soil Profiling, Mohsen M. 
Baligh, Vitoon Vivatrat and Charles C. Ladd, GT4: 
Apr 80:447:15377 


Consolidation (Soils) 

Hydroconsolidation Potential of Palouse Loess, Gary 
Lobdell, GT6:June 81:733:16309 

Undrained Settlement of Plastic and Organic Clays, 
Roger Foott and Charles C. Ladd, GT8:Aug 81: 
1079:16421 


1749 


CYCLIC LOADS 


Consolidation Tests (Soils) 

In-Situ Volume-Change Properties by Electro- 
Osmosis—Evaluation, James K. Mitchell and 
Sunirmal Banerjee, GT4:Apr 80:367:15371 


Construction Equipment 
Construction Vibrations: State-of-the-Art, John F. 
Wiss, GT2:Feb 81:167:16030 


Construction Materials 

Geotechnical Considerations for Construction in 
Saudi Arabia, Issa Oweis and John Bowman, GT3: 
Mar 81:319:16092 


Construction Methods 
There Were Giants on the Earth in Those Days, 
George F. Sowers, GT4:Apr 81:383:16190 


Construction Planning 

Geotechnical Considerations for Construction in 
Saudi Arabia, Issa Oweis and John Bowman, GT3: 
Mar 81:319:16092 


Construction Procedure 

Densification of Loose Deposits by Pounding, Robert 
G. Lukas, GT4:Apr 80:435:15376 

Geotechnical Considerations for Construction in 
Saudi Arabia, Issa Oweis and John Bowman, GT3: 
Mar 81:319:16092 

Ground Movements Caused by Braced Excavations, 
Thomas D. O'Rourke, GT9:Sept 81:1159:16511 


Cores (Dams) 

Cracking, Leakage, and Erosion of Earth Dam 
Materials, Erik |. Hjeldnes and Bhagwat V.K. 
Lavania, GT2:Feb 80:117:15220 


Cracking 

Cracking, Leakage, and Erosion of Earth Dam 
Materials, Erik I. Hjeldnes and Bhagwat V.K. 
Lavania, GT2:Feb 80:117:15220 


Cracking (Fracturing) 

Laboratory Study of Hydraulic Fracturing, Gary W. 
Jaworski, James M. Duncan and H. Bolton Seed, 
GT6:June 81:713:16287 


Creep 

Computer Simulation of Creep of Clay, Roland Pusch 
and Paul Feltham, GT1:Jan 81:95:15998 

Time-Dependant Behavior of Solution Caverns in 
Salt, Jameshid Ghaboussi, Afred J. Hendron, Jr. 
and Randall E. Ranken, GT10:Oct 81:1379:16597 

Undrained Settlement of Plastic and Organic Clays, 
Roger Foott and Charles C. Ladd, GT8:Aug 81: 
1079:16421 


Critical Depth 

Pressuremeter Tests at Very Shallow Depth, Jean- 
Louis Briaud and Donald H. Shields, GT8:Aug 81: 
1023:16416 


Cycles 
Cyclic Shear Resistance of Noncohesive Soils, 
Edward G. Prater, TN: GT1:Jan 80:111:15100 


Cyclic Loads 

Cyclic Axial Response of Single File, Harry G. 
Poulos, GT1:Jan 81:41:15979 

Cyclic Strengths Compared for Two Sampling 
Techniques, Carlos Espana, Ronald C. Chaney and 
Dennis Duffy, GTS:May 81:563:16234 

Dilation Angle and Liquefaction Potential, Yoginder 
P. Vaid, Peter M. Byrne and John M.O. Hughes, 
TN: GT7:July 81:1003:16360 

Earthquake-Induced Liquefaction Near Lake 
Amatitlan, Guatemala, H. Bolton Seed, Ignacio 
Arango, Clarence K. Chan, Alberto Gomez-Masso 
and Rebecca Grant Ascoli, GT4:Apr 81:501:16212 

Permanent Displacements Due to Cyclic Wave 
Loading, W. Allen Marr, Jr. and John T. Chrisitian, 
GT8:Aug 81:1129:16474 


CYCLIC LOADS 


Pile Load Tests: Cyclic Loads and Varying Load 
Rates, Leland M. Kraft, William R. Cox and 
Edward A. Verner, GT1:Jan 81:1:16000 

Soi! Failure Modes in Undrained Cyclic Loading, 
Ernest T. Selig and Ching S. Chang, GTS:May 81: 
$39:16238 


Cyclic Strength 

Cyclic Strengths Compared for Two Sampling 
Techniques, Carlos Espana, Ronald C. Chaney and 
Dennis Duffy, GTS:May 81:563:16234 


Dam Failure 

Dam Failure in Benefit/Cost Analysis, Gregory 
Baecher, Marie-Elisabeth Pate and Richard de 
Neufville, TN: GT1:Jan 80:101:15100 


Damping 

Embedded Foundations Under Vertical Vibration, A. 
Sridharan, M. V. Nagendra and C. Chinnaswamy, 
TN: GT10:Oct 


Dams 

Laboratory Study of Hydraulic Fracturing, Gary W. 
Jaworski, James M. Duncan and H. Bolton Seed, 
GT6:June 81:713:16287 


Dams (Earth) 

Comparative Study of Dynamic Response of Earth 
Dam, Ahmed M. Abdel-Ghaffar and Ronald F. 
Scott, GT3:Mar 81:271:16097 

Cracking, Leakage, and Erosion of Earth Dam 
Materials, Erik |. Hjeldnes and Bhagwat V.K. 
Lavania, GT2:Feb 80:117:15220 

Dynamic FEM Model of Oroville Dam, John 
Vrymoed, GT8:Aug 81:1057:16464 

Seepage From Freewater Above Impermeable 
Tailings, Lewis T. Isaacs and Bruce Hunt, GT11: 
Nov 81:1563:16652 

Stability Analysis of Embankments and Slopes, 
Sarada K. Sarma, GT12:Dec 79:1511:15068 

Vibration Tests of Full-Scale Earth Dam, Ahmed M. 
Abdel-Ghaffar and Ronald F. Scott, GT3:Mar 81: 
241:16096 


Dams (Rockfill) 

Northfield Mountain Pumped Storage Project: 
Performance, William F. Swiger, Philip A. Wild and 
Thomas J. Lamb, GT6:June 80:673:15506 


Decomposition 

Ignition Test for Soil Organic-Content Measurement, 
Abdul-Amir W.N. Al-Kafaji and Orlando B. 
Andersland, GT4:Apr 81:465:16174 


Deflection 

Microgravity Surveys for Evaluation of Elevation 
Changes Due to Reservoir Impoundment, Dwain K. 
Butler, TN: GT3:Mar 81:355:16073 


Deformation 

Permanent Displacements Due to Cyclic Wave 
Loading, W. Allen Marr, Jr. and John T. Chrisitian, 
GT8:Aug 81:1129:16474 

The Steady State of Deformation, Steve J. Poulos, 
GTS:May 81:553:16241 

Stochastic FEM In Settlement Predictions, Gregory 
B. Beacher and Thomas S. Ingra, GT4:Apr 81:449: 
16179 

Time-Dependent Deformation Behavior of Clays, 
Edward Kavazanjian, Jr. and James K. Mitchell, 
GT6:June 80:61 1:15488 

Undrained Settlement of Plastic and Organic Clays, 
Roger Foott and Charles C. Ladd, GT8:Aug 81: 
1079:16421 


Deformation Analysis 

Strength and Deformability of Highly Fractured 
Rock, Jerome M. Raphael and Richard E. 
Goodman, GT11:Nov 79:1285:14988 


1981 


Deposits 

Statistical Analysis of Marine Clay Deposits, Hin 
Fatt Cheong and R. V. Subrahmanyam, TN: GT2: 
Feb 81:221:16009 


Design 

Design Considerations for Collapsible Soils, Samuel 
P. Clemence and Albert O. Finbarr, GT3:Mar 81: 
305:16106 


Design Criteria 

Current USA Practices: Slurry Wall Specifications, 
Richard A. Millet and Jean-Yves Perez, GT8:Aug 
81:1041:16458 

New Design Correlations for Piles in Sand, Harry M. 
Coyle and Reno R. Castello, GT7:July 81:965: 
16379 


Design Practices 

Design Considerations for Collapsible Soils, Samuel 
P. Clemence and Albert O. Finbarr, GT3:Mar 81: 
305:16106 


Dielectrics 

Relationship of Electrical Dispersion to Soil 
Properties, Scott S. Smith and K. Arulanandan, 
GTS5:May 81:591:16242 


Displacements 

Permanent Displacements Due to Cyclic Wave 
Loading, W. Allen Marr, Jr. and John T. Chrisitian, 
GT8:Aug 81:1129:16474 


Dunes 

Measurement of Porosity in Natural Sand Desposits, 
S. A. Denekamp and Y. Tsur-Lavie, GT4:Apr 81: 
439:16173 


Dynamic Characteristics 

Comparative Study of Dynamic Response of Earth 
Dam, Ahmed M. Abdel-Ghaffar and Ronald F. 
Scott, GT3:Mar 81:271:16097 

Vibration Tests of Full-Scale Earth Dam, Ahmed M. 
Abdel-Ghaffar and Ronald F. Scott, GT3:Mar 81: 
241:16096 


Dynamic Models 
Dynamic FEM Model of Oroville Dam, John 
Vrymoed, GT8:Aug 81:1057:16464 


Dynamic Properties 

Dynamic Properties of Soils from In-Situ Tests, 
Shamsher Prakash and Vijay Kumar Puri, GT7:July 
81:943:16366 

Field and Laboratory Determination of Soil Moduli, 
William F. Marcuson, III and Joseph R. Curro, Jr., 
GT10:Oct 81:1269:16591 


Dynamic Response 

Comparative Study of Dynamic Response of Earth 
Dam, Ahmed M. Abdel-Ghaffar and Ronald F. 
Scott, GT3:Mar 81:271:16097 

Nonlinear Lateral Dynamic Stiffness of Piles, 
Demosthenes C. Angelides and Jose M. Roesset, 
GTI11:Nov 81:1443:16635 


Dynamics 

Construction Vibrations: State-of-the-Art, John F. 
Wiss, GT2:Feb 81:167:16030 

Nonlinear Lateral Dynamic Stiffness of Piles, 
Demosthenes C. Angelides and Jose M. Roesset, 
GTI1I1:Nov 81:1443:16635 


Dynamic Stability 

Vibration Tests of Full-Scale Earth Dam, Ahmed M. 
Abdel-Ghaffar and Ronald F. Scott, GT3:Mar 81: 
241:16096 


Dynamic Tests 

Comparative Study of Dynamic Response of Earth 
Dam, Ahmed M. Abdel-Ghaffar and Ronald F. 
Scott, GT3:Mar 81:271:16097 


1750 


SUBJECT INDEX 


Strength and Deformability of Highly Fractured 
Rock, Jerome M. Raphael and Richard E. 
Goodman, GT11:Nov 79:1285:14988 


Earth Dam Performance 
Line of Seepage in Earth Dams on Inclined Ledge, 
Yang H. Huang, TN: GTS:May 81:662:16221 


Earth Dams 

Dynamic FEM Model of Oroville Dam, 
Vrymoed, GT8:Aug 81:1057:16464 

Line of Seepage in Earth Dams on Inclined Ledge, 
Yang H. Huang, TN: GTS:May 81:662:16221 

Longitudinal Vibrations of Embankment 
George Gazetas, GT1:Jan 81:21:15980 

Seismic Displacement Analysis of Earth Dams, 
Sarada K. Sarma, TN: GT12:Dec 81:1735:16688 


John 


Dams, 


Earth Dam Seepage 
Line of Seepage in Earth Dams on Inclined Ledge, 
Yang H. Huang, TN: GTS:May 81:662:16221 


Earthfill 
Densification of Loose Deposits by Pounding, Robert 
G. Lukas, GT4:Apr 80:435:15376 


Earth Pressure 

In Situ Tests by Flat Dilatometer, Silvano Marchetti, 
GT3:Mar 80:299:15290 

Laterally Loaded Pile Design, Robert L. Sogge, GT9: 
Sept 81:1179:16510 

Principal Stress Rotation: A Missing Parameter, J. 
Robin F. Arthur, Ken S. Chua, Treve Dunstan and 
Juan |. Rodriguez del C., GT4:Apr 80:419:15375 


Earthquake Damage 
Seismic Displacement Analysis of Earth Dams, 
Sarada K. Sarma, TN: GT12:Dec 81:1735:16688 


Earthquake Engineering 
Evaluating Pile Performance During Earthquakes, 
Issa S. Oweis, TN: GTS:May 81:678:16221 


Earthquakes 
Dynamic FEM Model 
Vrymoed, GT8:Aug 81:1057:16464 


of Oroville Dam, John 

Earthquake-Induced Liquefaction Near Lake 
Amatitlan, Guatemala, H. Bolton Seed, Ignacio 
Arango, Clarence K. Chan, Alberto Gomez-Masso 
and Rebecca Grant Ascoli, GT4:Apr 81:501:16212 

Evaluating Pile Performance During Earthquakes, 
Issa S. Oweis, TN: GTS:May 81:678:16221 

Lateral Pile Response During Earthquakes, Takaaki 
Kagawa and Leland M. Kraft, Jr.. GT12:Dec 81: 
1713:16735 

Liquefaction Analysis for Multidirectional Shaking, 
Jamshid Ghaboussi and S. Umit Dikmen, GT5:May 
81:605:16243 

Longitudinal Vibrations of Embankment 
George Gazetas, GT1:Jan 81:21:15980 

Probabilistic Evaluation of Loads, Wilson H. Tang, 
GT3:Mar 81:287:16100 

Stability Analysis of Embankments and Slopes, 
Sarada K. Sarma, GT12:Dec 79:1511:15068 

Statistical Study of Uniform Cycles in Earthquakes, 
Achintya Haldar and Wilson H. Tang, GTS:May 
81:577:16239 


Earth Tides 

Microgravity Surveys for Evaluation of Elevation 
Changes Due to Reservoir Impoundment, Dwain K. 
Butler, TN: GT3:Mar 81:355:16073 


Dams, 


Elasticity 

Laboratory Tests on Model Piled Raft Foundations, 
Terence J. Wiesner and Peter T. Brown, GT7:July 
80:767:15576 


Electrical Properties 

Relationship of Electrical Dispersion to Soil 
Properties, Scott S. Smith and K. Arulanandan, 
GTS5:May 81:591:16242 


1751 


FABRICS 


Electroosmosis 
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New Design Correlations for Piles in Sand, Harry M. 
Coyle and Reno R. Castello, GT7:July 81:965: 
16379 


Pile Foundation Design 

New Design Correlations for Piles in Sand, Harry M. 
Coyle and Reno R. Castello, GT7:July 81:965: 
16379 


Pile Foundation Performance 
Evaluating Pile Performance During Earthquakes, 
Issa S. Oweis, TN: GTS:May 81:678:16221 


Pile Friction 

Friction Capacity of Piles Driven into Clay, Leland 
M. Kraft, John A. Focht, Jr. and Srinath F. 
Amerasinghe, GT11:Nov 81:1521:16663 


Pile Loading Tests 
Piles Subjected to Torsion, M. F. Randolph, GTS8: 
Aug 81:1095:16424 


Pile Load Tests (Cyclic Loading) 

Pile Load Tests: Cyclic Loads and Varying Load 
Rates, Leland M. Kraft, William R. Cox and 
Edward A. Verner, GT1:Jan 81:1:16000 


Piles 

Laboratory Tests on Model Piled Raft Foundations, 
Terence J. Wiesner and Peter T. Brown, GT7:July 
80:767:15576 

Laterally Loaded Pile Design, Robert L. Sogge, GT9: 
Sept 81:1179:16510 

Repeated Loading on Single Piles in Sand, Sin-Fatt 
cae and Thomas H. Hanna, GT2:Feb 80:171: 
1$222 

Stresses in Soil around Vertical Compressible Piles, 

K. S. Sankaran, N. R. Krishnaswamy and B. K. 

Sharas Chandra, TN: GT1:Jan 81:107:15956 


Piles (Foundations) 


Piles Subjected to Torsion, M. F. Randolph, GT8: 
Aug 81:1095:16424 


Pile Settlement 


Cyclic Axial Response of Single File, Harry G. 


Poulos, GT1:Jan 81:41:15979 
Theoretical t-z Curves, Leland M. Kraft, Jr., Richard 


P. Ray and Takaaki Kagawa, GT1I1:Nov 81:1543: 
16653 


Pipelines 
Response of Buried Structures to Traveling Waves, 
Richard N. Nwang and John Lysmer, GT2:Feb 81: 
183:16052 


Piping (Erosion) 

Permeability and Piping in Fractured Rocks, Richard 
E. Goodman and Panchanatham N. Sundaram, 
GTS5:May 80:485:15433 


Pneumatic Tired Rollers 
Construction of Large Canal on Collapsing Soils, Paul 
C. Knodel, GT1:Jan 81:79:15992 


Pondweeds 


Construction of Large Canal on Collapsing Soils, Paul 
C. Knodel, GT1:Jan 81:79:15992 


Pore Water Pressure 


Cone Penetration in Soil Profiling, Mohsen M. 


Baligh, Vitoon Vivatrat and Charles C. Ladd, GT4: 
Apr 80:447:15377 

Lateral Pile Response During Earthquakes, Takaaki 
Kagawa and Leland M. Kraft, Jr.. GT12:Dec 81: 
1713:16735 

Modeling the Liquefaction Process, Takaaki Kagawa 
and Leland M. Kraft, Jr., GT12:Dec 81:1593:16709 


Porosity 
Measurement of Porosity in Natural Sand Desposits, 


S. A. Denekamp and Y. Tsur-Lavie, GT4:Apr 81: 
439:16173 
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1981 


Portland Cements 
Soil Stabilization with Preheater Fines, Sudipta S. 
Bandyopadhyay, TN: GTS:May 81:654:16221 


Predicting 

Theoretical t-z Curves, Leland M. Kraft, Jr., Richard 
P. Ray and Takaaki Kagawa, GT11:Nov 81:1543: 
16653 


Preheating 
Soil Stabilization with Preheater Fines, Sudipta S. 
Bandyopadhyay, TN: GTS:May 81:654:16221 


Pressure 

Self-Boring Pressuremeter Study of San Francisco 
Bay Mud, G. Wayne Clough and Gordon M. 
Denby, GT1:Jan 80:45:15148 


Pressure Measurement 

Pressuremeter Tests at Very Shallow Depth, Jean- 
Louis Briaud and Donald H. Shields, GT8:Aug 81: 
1023:16416 


Probability Theory 

Probabilistic Evaluation of Loads, Wilson H. Tang, 
GT3:Mar 81:287:16100 

Stochastic FEM In Settlement Predictions, Gregory 
Beacher and Thomas Ingra, GT4:Apr 81:449: 
16179 


Pumped Storage 

Northfield Mountain Pumped Storage Project: 
Performance, William F. Swiger, Philip A. Wild and 
Thomas J. Lamb, GT6:June 80:673:15506 


Rayleigh Waves 

Seismic Techniques in the Laboratory, Richard D. 
Woods and Robert Henke, GT10:Oct 81:1309: 

16599 


Reinforced Earth 

Field Measurements of an Earth Support System, C. 
K. Shen, S. Bang, K. M. Romstad, L. Kulchin and 
J. S. DeNatale, GT12:Dec 81:1625:16734 

Ground Movement Analysis of Earth Support 
System, C. K. Shen, S. Bang and L. R. Herrmann, 
GT12:Dec 81:1609: 16732 

Stability of Loaded Footings on Reinforced Soil, Joe 

O. Akinmusuru and Jones A. Akinbolade, GT6: 

June 81:819:16320 


Reinforced Roads 

Geotextile-Reinforced Unpaved Road Design, Jean- 
Pierre Giroud and Laure Noiray, GT9:Sept 81: 
1233:16489 


Relative Density 

Estimation of SPT-N and Relative Density, Michael 
A. Fardis and Daniele Veneziano, GT10:Oct 81: 
1345:16590 

Measurement of Porosity in Natural Sand Desposits, 
S. A. Denekamp and Y. Tsur-Lavie, GT4:Apr 81: 
439:16173 


Repeated Loading 

Repeated Loading on Single Piles in Sand, Sin-Fatt 
Chan and Thomas H. Hanna, GT2:Feb 80:171: 
15222 


Reservoirs 

Microgravity Surveys for Evaluation of Elevation 
Changes Due to Reservoir Impoundment, Dwain K. 
Butler, TN: GT3:Mar 81:355:16073 


Residual Shear Strength 
The Steady State of Deformation, Steve J. Poulos, 
81:553:16241 


Residual Stress 
New Design Correlations for Piles in Sand, Harry M. 


Coyle and Reno R. Castello, GT7:July 81:965: 
16379 


SUBJECT INDEX 


Responses 

Theoretical t-z Curves, Leland M. Kraft, Jr., Richard 
P. Ray and Takaaki Kagawa, GT11:Nov 81:1543: 
16653 


Retaining Walls 
Total Lateral Surcharge Pressure Due to Strip Load, 
Ramon Jarquio, TN: GT10:Oct 81:1424:16530 


Rheology 

Road on Peat: Observations and Design, Edward T. 
Hanrahan and Martin G. Rogers, GT10:Oct 81: 
1403:16588 


Road Construction 

Road on Peat: Observations and Design, Edward T. 
Hanrahan and Martin G. Rogers, GT10:Oct 81: 
1403:16588 


Road Design 

Geotextile-Reinforced Unpaved Road Design, Jean- 
Pierre Giroud and Laure Noiray, GT9:Sept 81: 
1233:16489 

Road on Peat: Observations and Design, Edward T. 
Hanrahan and Martin G. Rogers, GT10:Oct 81: 
1403:16588 


Rock Analysis 

Miami Limestone Foundation Design and 
Construction, Thomas J. Kaderabek and Richard T. 
Reynolds, GT7:July 81:859:16377 


Rockfill Dam Design 
Shear Strength of Rockfill, Nick Barton and Bjorn 
Kjaernsli, GT7:July 81:873:16374 


Rock Mechanics 

Northfield Mountain Pumped Storage Project: 
Performance, William F. Swiger, Philip A. Wild and 
Thomas J. Lamb, GT6:June 80:673:15506 

Strength and Deformability of Highly Fractured 
Rock, Jerome M. Raphael and Richard E. 
Good> an, GT11:Nov 79:1285:14988 


Rock Properties 

Strength and Deformability of Highly Fractured 
Rock, Jerome M. Raphael and Richard E. 
Goodman, GT11:Nov 79:1285:14988 


Rock Salt 

Time-Dependant Behavior of Solution Caverns in 
Salt, Jameshid Ghaboussi, Afred J. Hendron, Jr. 
and Randall E. Ranken, GT10:Oct 81:1379:16597 


Rock Strength 

Strength and Deformability of Highly Fractured 
Rock, Jerome M. Raphael and Richard E. 
Goodman, GT11:Nov 79:1285:14988 


Rock Tests 

Miami Limestone Foundation Design and 
Construction, Thomas J. Kaderabek and Richard T. 
Reynolds, GT7:July 81:859:16377 


Roughness 
Rock Joints: Roughness-Shear Strength Relationship, 
Ralf Peek, TN: GTS:May 81:672:16221 


Safety 

Safety of a Constructed Facility: Geotechnical 
Aspects, T. William Lambe, W. Allen Marr and 
Franciso Silva, GT3:Mar 81:339:16107 


Safety Engineering 

Safety of a Constructed Facility: Geotechnical 
Aspects, T. William Lambe, W. Allen Marr and 
Franciso Silva, GT3:Mar 81:339:16107 


Safety Factor 

Limit Plasticity Approach to Slope Stability Program, 
Ryszard J. Izbicki, TN: GT2:Feb 81:228:16009 

Line of Seepage in Earth Dams on Inclined Ledge, 
Yang H. Huang, TN: GTS:May 81:662:16221 


SANDY SOILS 


Random Surface Generation in Stability Analysis, 
Ronald A. Siegel, William D. Kovacs and Charles 
W. Lovell, TN: GT7:July 81:996:16360 


Salt Domes 

Time-Dependant Behavior of Solution Caverns in 
Salt, Jameshid Ghaboussi, Afred J. Hendron, Jr. 
and Randall E. Ranken, GT10:Oct 81:1379:16597 


Sampling 

Cyclic Strengths Compared for Two Sampling 
Techniques, Carlos Espana, Ronald C. Chaney and 
Dennis Duffy, GTS:May 81:563:16234 


Sand 

Dilation Angle and Liquefaction Potential, Yoginder 
P. Vaid, Peter M. Byrne and John M.O. Hughes, 
TN: GT7:July 81:1003:16360 

Earthquake-Induced Liquefaction Near Lake 
Amatitlan, Guatemala, H. Bolton Seed, Ignacio 
Arango, Clarence K. Chan, Alberto Gomez-Masso 
and Rebecca Grant Ascoli, GT4:Apr 81:501:16212 

Measurement of Porosity in Natural Sand Desposits, 
S. A. Denekamp and Y. Tsur-Lavie, GT4:Apr 81: 
439:16173 

New Design Correlations for Piles in Sand, Harry M. 
Coyle and Reno R. Castello, GT7:July 81:965: 
16379 

Soil Failure Modes in Undrained Cyclic Loading, 
Ernest T. Selig and Ching S. Chang, GTS:May 81: 
539:16238 

Stability of Loaded Footings on Reinforced Soil, Joe 
O. Akinmusuru and Jones A. Akinbolade, GT6: 
June 81:819:16320 


Sands 

Anisotropic Strength of Cohesionless 
Masanobu Oda, GT9:Sept 81:1219:16491 

Cemented Sands Under Static Loading, G. Wayne 
Clough, Nicholas Sitar, Robert C. Bachus and 
Nader Shafii Rad, GT6:June 81:799:16319 

Compressibility and Bearing Capacity, Nabil F. 
Ismael and Aleksandar S. Vesic, GT12:Dec 81: 
1677:16718 

Cyclic Strengths Compared for Two Sampling 
Techniques, Carlos Espana, Ronald C. Chaney and 
Dennis Duffy, GTS:May 81:563:16234 

Estimation of SPT-N and Relative Density, Michael 
A. Fardis and Daniele Veneziano, GT10:Oct 81: 
1345:16590 

Lateral Pile Response During Earthquakes, Takaaki 
Kagawa and Leland M. Kraft, Jr.. GT12:Dec 81: 
1713:16735 

Modeling the Liquefaction Process, Takaaki Kagawa 
and Leland M. Kraft, Jr., GT12:Dec 81:1593:16709 

Permanent Displacements Due to Cyclic Wave 
Loading, W. Allen Marr, Jr. and John T. Chrisitian, 
GT8:Aug 81:1129:16474 

Principal Stress Rotation: A Missing Parameter, J. 
Robin F. Arthur, Ken S. Chua, Treve Dunstan and 
Juan I. Rodriguez del C., GT4:Apr 80:419:15375 

Resorcinolic Grout for Injecting Sandy Foundations, 
Arvind V. Shroff and Dhananjay L. Shah, TN: 
GT10:Oct 80:1153:15724 

Static Shear and Liquefaction Potential, Yoginder P. 
Vaid and W. D. Liam Finn, GT10:Oct 79:1233: 
14909 

Statistical Analysis of Sand Liquefaction, Michael N. 
Fardis and Daniele Veneziano, GT10:Oct 81:1361: 
16604 


Sand Waves 

Sand Waves in Lower Cook Inlet, Alaska, Arshud A. 
Mahmood, Clarence J. Ehlers and Blase A. 
Cilweck, GT10:Oct 81:1293:16583 


Sandy Soils 

Earthquake-Induced Liquefaction Near Lake 
Amatitlan, Guatemala, H. Bolton Seed, Ignacio 
Arango, Clarence K. Chan, Alberto Gomez-Masso 
and Rebecca Grant Ascoli, GT4:Apr 81:501:16212 

Statistical Analysis of Sand Liquefaction, Michael N. 
Fardis and Daniele Veneziano, GT10:Oct 81:1361: 
16604 


Sands, 


SAUDI ARABIA 


Saudi Arabia 

Geotechnical Considerations for Construction in 
Saudi Arabia, Issa Oweis and John Bowman, GT3: 
Mar 81:319:16092 


Seepage 

Acoustic Emission Monitoring of Seepage, Robert M. 
Koerner, W. Martin McCabe and Luis F. 
Baldivieso, TN: GT4:Apr 81:521:16156 

Centrifuge Modeling of Coal Waste Embankments, 
Mosaid M. Al-Hussaini, Deborah J. Goodings, 
Andrew N. Schofield and Frank C. Townsend, 
GT4:Apr 81:481:16180 

Line of Seepage in Earth Dams on Inclined Ledge, 
Yang H. Huang, TN: GTS:May 81:662:16221 

Seepage From  Freewater Above Impermeable 
Tailings, Lewis T. Isaacs and Bruce Hunt, GT11: 
Nov 81:1563:16652 


Seepage Control 
Soil-Bentonite Slurry Trench Cutoffs, David J. 
D'Appolonia, GT4:Apr 80:399:15372 


Seepage Theory 

Acoustic Emission Monitoring of Seepage, Robert M. 
Koerner, W. Martin McCabe and Luis F. 
Baldivieso, TN: GT4:Apr 81:521:16156 


Seismic Design 

Response of Buried Structures to Traveling Waves, 
Richard N. Nwang and John Lysmer, GT2:Feb 81: 
183:16052 

Seismic Displacement Analysis of Earth Dams, 
Sarada K. Sarma, TN: GT12:Dec 81:1735:16688 


Seismic Stability 

Nonlinear Seismic Response of Soft Clay Sites, Ram 
D. Singh, Ricardo Debry, Earl H. Doyle and Izzat 
M. Idriss, GT9:Sept 81:1201:16493 

Seismic Displacement Analysis of Earth Dams, 
Sarada K. Sarma, TN: GT12:Dec 81:1735:16688 


Seismic Tests 

Seismic Techniques in the Laboratory, Richard D. 
Woods and Robert Henke, GT10:Oct 81:1309: 
16599 


Seismic Waves 

Instrumentation Corrections to Wave Velocity Data, 
John R. Hall, Jr., John T. Christian, Eduardo 
Kausel and James L. Wolfgang, TN: GT10:Oct 81: 
1419:16530 


Seismographs 

Instrumentation Corrections to Wave Velocity Data, 
John R. Hall, Jr., John T. Christian, Eduardo 
Kausel and James L. Wolfgang, TN: GT10:Oct 81: 
1419:16530 


Sensitivity 

Nitrogen-Cycle Model for Aquatic Systems: Analysis, 
Tavit O. Najarian and Jay L. Taft, GT12:Dec 81: 
1141:16703 


Settlement (Structural) 

Laboratory Tests on Model Piled Raft Foundations, 
Terence J. Wiesner and Peter T. Brown, GT7:July 
80:767:15576 

Settlement Comparison Used in Tank-Failure Study, 
Roy A. Bell and Jun Iwakiri, GT2:Feb 80:153: 
15219 

Stochastic FEM In Settlement Predictions, Gregory 
B. Beacher and Thomas S. Ingra, GT4:Apr 81:449: 
16179 


Tolerable Settlement of Buildings, Harvey E. Wahls, 
GT11:Nov 81:1489:16628 


Settlement Analysis 

Undrained Settlement of Plastic and Organic Clays, 
Roger Foott and Charles C. Ladd, GT8:Aug 81: 
1079:16421 


1981 


Shadows 
Rock Joints: Roughness-Shear Strength Relationship, 
Ralf Peek, TN: GTS:May 81:672:16221 


Shafts 

Drilled-Shaft Integrity by Wave Propagation Method, 
Thomas M. Hearne, Jr., Kenneth H. Stokoe, II and 
Lyndon C. Reese, GT10:Oct 81:1327:16582 

Repeated Loading on Single Piles in Sand, Sin-Fatt 
Chan and Thomas H. Hanna, GT2:Feb 80:171: 
15222 


Shaking 

Liquefaction Analysis for Multidirectional Shaking, 
Jamshid Ghaboussi and Umit Dikmen, 
81:605: 16243 


Shallow Foundations 

Experimental Study on Footings in Layered Soil, 
Adel M. Hanna, GT8:Aug 81:1113:16447 

Miami Limestone Foundation Design and 
Construction, Thomas J. Kaderabek and Richard T. 
Reynolds, GT7:July 81:859:16377 


Shearing 
Slip Circles and Critical Shear Planes, Eric Spencer, 
GT7:July 81:929:16369 


Shear Modulus 

Field and Laboratory Determination of Soil Moduli, 
William F. Marcuson, III and Joseph R. Curro, Jr., 
GT10:Oct 81:1269:16591 

Seismic Techniques in the Laboratory, Richard D. 
Woods and Robert Henke, GT10:Oct 81:1309: 
16599 

Undrained Settlement of Plastic and Organic Clays, 
Roger Foott and Charles C. Ladd, GT8:Aug 81: 
1079:16421 


Shear Planes 
Slip Circles and Critical Shear Planes, Eric Spencer, 
GT7:July 81:929:16369 


Shear Resistance 
Cyclic Shear Resistance of Noncohesive Soils, 
Edward G. Prater, TN: GT1:Jan 80:111:15100 


Shear Strain 

Cyclic Axial Response of Single File, Harry G. 
Poulos, GT1:Jan 81:41:15979 

Evaluating Pile Performance During Earthquakes, 
Issa S. Oweis, TN: GTS:May 81:678:16221 


Shear Strength 

Effect of Wetting and Drying on Shear Strength, 
Mehter Mohamed Allam and Asuri Sridharan, 
GT4:Apr 81:421:16178 

Shear Strength of Rockfill, Nick Barton and Bjorn 
Kjaernsli, GT7:July 81:873:16374 


Shear Strength (Rock) 
Rock Joints: Roughness-Shear Strength Relationship, 
Ralf Peek, TN: GTS:May 81:672:16221 


Shear Strength (Soils) 

Clay Anisotropy and Braced Wall Behavior, G. 
Wayne Clough and Lawrence A. Hansen, GT7:July 
81:893:16391 

Probabilistic Soil Exploration: Case History, Tien H. 
Wu and Kinfun Wong, GT12:Dec 81:1693:16764 


Shear Stress 

Evaluating Pile Performance During Earthquakes, 
Issa S. Oweis, TN: GTS:May 81:678:16221 

Modeling the Liquefaction Process, Takaaki Kagawa 
and Leland M. Kraft, Jr., GT12:Dec 81:1593:16709 


Shear Tests 

Dilation Angle and Liquefaction Potential, Yoginder 
P. Vaid, Peter M. Byrne and John M.O. Hughes, 
TN: GT7:July 81:1003:16360 

Principal Stress Rotation: A Missing Parameter, J. 
Robin F. Arthur, Ken S. Chua, Treve Dunstan and 
Juan I. Rodriguez del C., GT4:Apr 80:419:15375 


SUBJECT INDEX 


Shear Strength of Rockfill, Nick Barton and Bjorn 
Kjaernsli, GT7:July 81:873:16374 

Statistical Analysis of Sand Liquefaction, Michael N. 
caone and Daniele Veneziano, GT10:Oct 81:1361: 


Shrinkage 

Methodology for Foundations on Expansive Clays, 
Michael W. O'Neill and Nader Poormoayed, GT12: 
Dec 80:1345:15949 


Simple Shear Tests 

Static Shear and Liquefaction Potential, Yoginder P. 
Vaid and W. D. Liam Finn, GT10:Oct 79:1233: 
14909 


Sinkage 


Anchored Bulkheads: Horizontal and Sloping 


Anchors, Boris S. Browzin, GTS:May 81:629:16275 


Site Evaluation 

Nonlinear Seismic Response of Soft Clay Sites, Ram 
D. Singh, Ricardo Debry, Earl H. Doyle and Izzat 
M. Idriss, GT9:Sept 81:1201:16493 


Site Preparation (Construction) 
Densification of Loose Deposits by Pounding, Robert 
G. Lukas, GT4:Apr 80:435:15376 


Slabs 

Methodology for Foundations on Expansive Clays, 
Michael W. O'Neill and Nader Poormoayed, GT12: 
Dec 80:1345:15949 


Slipping 

Extreme-Value Problems of Limiting Equilibrium, 
Michael Garber and Rafael Baker, GT10:Oct 79: 
1155:14901 


Slopes 
Limit Plasticity Approach to Slope Stability Program, 
Ryszard J. Izbicki, TN: GT2:Feb 81:228:16009 


Slope Stability 

Centrifuge Modeling of Coal Waste Embankments, 
Mosaid M. Al-Hussaini, Deborah J. Goodings, 
Andrew N. Schofield and Frank C. Townsend, 
GT4:Apr 81:481:16180 

Drum Centrifuge Studies of Over Consolidated 
Slopes, Richard J. Fragaszy and James A. Cheney, 
GT7:July 81:843:16362 

Extreme-Value Problems of Limiting Equilibrium, 
Michael Garber and Rafael Baker, GT10:Oct 79: 
1155:14901 

Random Surface Generation in Stability Analysis, 
Ronald A. Siegel, William D. Kovacs and Charles 
W. Lovell, TN: GT7:July 81:996:16360 

Slip Circles and Critical Shear Planes, Eric Spencer, 
GT7:July 81:929:16369 

Stability Analysis of Embankments and Slopes, 
Sarada K. Sarma, GT12:Dec 79:1511:15068 


Slope Stability Analysis 

Random Surface Generation in Stability Analysis, 
Ronald A. Siegel, William D. Kovacs and Charles 
W. Lovell, TN: GT7:July 81:996:16360 

Slip Circles and Critical Shear Planes, Eric Spencer, 
GT7:July 81:929:16369 


Slurry Excavation 

Current USA Practices: Slurry Wall Specifications, 
Richard A. Millet and Jean-Yves Perez, GT8:Aug 
81:1041:16458 


Slurry Trenches 

Current USA Practices: Slurry Wall Specifications, 
Richard A. Millet and Jean-Yves Perez, GT8:Aug 
81:1041:16458 

Soil-Bentonite Slurry Trench Cutoffs, 
D’Appolonia, GT4:Apr 80:399:15372 


Sociological Factors 
There Were Giants on the Earth in Those Days, 
George F. Sowers, GT4:Apr 81:383:16190 


David J. 
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SOIL MECHANICS 


Soft Soils 

Nonlinear Seismic Response of Soft Clay Sites, Ram 
D. Singh, Ricardo Debry, Earl H. Doyle and Izzat 
M. Idriss, GT9:Sept 81:1201:16493 

Self-Boring Pressuremeter Study of San Francisco 
Bay Mud, G. Wayne Clough and Gordon M. 
Denby, GT1:Jan 80:45:15148 


Soil Analysis 
Optimal Estimators for Soil Properties, Gregory B. 
Baecher, TN: GTS:May 81:649:16221 


Soil Cement 

Relationship of Electrical Dispersion to Soil 
Properties, Scott S. Smith and K. Arulanandan, 
GTS:May 81:591:16242 


Soil Compression Tests 

At-Rest Lateral Pressure of Peat Soils, Tuncer B. 
Edil and Abdulmohsin W. Dhowian, GT2:Feb 81: 
201:16063 


Soil Consolidation Tests 
Hydroconsolidation Potential of Palouse Loess, Gary 
T. Lobdell, GT6:June 81:733:16309 


Soil Density 

Estimation of SPT-N and Relative Density, Michael 
A. Fardis and Daniele Veneziano, GT10:Oct 81: 
1345:16590 


Soil Dynamics 

Dynamic Properties of Soils from In-Situ Tests, 
Shamsher Prakash and Vijay Kumar Puri, GT7:July 
81:943:16366 

Ignition Test for Soil Organic-Content Measurement, 
Abdul-Amir W.N. Al-Kafaji and Orlando B. 
Andersland, GT4:Apr 81:465:16174 

Liquefaction Analysis for Multidirectional Shaking, 
Jamshid Ghaboussi and Umit Dikmen, 
81:605:16243 

Prediction of Swelling Potential for Natural Soils, 
Sudipta S. Bandyopadhyay, TN: GT5:May 81:658: 
16221 

Seismic Techniques in the Laboratory, Richard D. 
Woods and Robert Henke, GT10:Oct 81:1309: 
16599 

Soil Failure Modes in Undrained Cyclic Loading, 
Ernest T. Selig and Ching S. Chang, GTS:May 81: 
539:16238 

Undrained Settlement of Plastic and Organic Clays, 
Roger Foott and Charles C. Ladd, GT8:Aug 81: 
1079:16421 


Soil Engineering 

Relationship of Electrical Dispersion to Soil 
Properties, Scott S. Smith and K. Arulanandan, 
GTS5:May 81:591:16242 


Soil Investigations 

Cone Penetration in Soil Profiling, Mohsen M. 
Baligh, Vitoon Vivatrat and Charles C. Ladd, GT4: 
Apr 80:447:15377 

Probabilistic Soil Exploration: Case History, Tien H. 
Wu and Kinfun Wong, GT12:Dec 81:1693:16764 


Soil Liquefaction 

Static Shear and Liquefaction Potential, Yoginder P. 
Vaid and W. D. Liam Finn, GT10:Oct 79:1233: 
14909 


Soil Mechanics 

Cemented Sands Under Static Loading, G. Wayne 
Clough, Nicholas Sitar, Robert C. Bachus and 
Nader Shafii Rad, GT6:June 81:799:16319 

Effect of Wetting and Drying on Shear Strength, 
Mehter Mohamed Allam and Asuri Sridharan, 
GT4:Apr 81:421:16178 

Embedded Anchor Response to uplift Loading, 
Andreas Andreadis, Roger Harvey and Eldon 
Burley, GT1:Jan 81:59:15985 

Experimental Study on Footings in Layered Soil, 
Adel M. Hanna, GT8:Aug 81:1113:16447 


SOIL MECHANICS 


Geotechnical Considerations for Construction in 
Saudi Arabia, Issa Oweis and John Bowman, GT3: 
Mar 81:319:16092 

Ignition Test for Soil Organic-Content Measurement, 
Abdul-Amir W.N. Al-Kafaji and Orlando B. 
Andersland, GT4:Apr 81:465:16174 

Laterally Loaded Pile Design, Robert L. Sogge, GT9: 
Sept 81:1179:16510 

Seepage From Freewater Above Impermeable 
Tailings, Lewis T. Isaacs and Bruce Hunt, GT11: 
Nov 81:1563:16652 

Slip Circles and Critical Shear Planes, Eric Spencer, 
GT7:July 81:929:16369 

Soil-Steel Structure Response to Live Loads, Baidar 
Bakht, GT6:June 81:779:16316 


Soil Moisture 

Coupled Heat and Water Flows Around Buried 
Cables, Omar N. Abdel-Hadi and James K. 
Mitchell, GT11:Nov 81:1461:16634 


Soil Profiles 
Optimal Estimators for Soil Properties, Gregory B. 
Baecher, TN: GTS:May 81:649:16221 


Soil Properties 
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$1:1041:16458 


Stability 

Anchored Bulkheads: Horizontal and Sloping 
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Instability of Amuay Cliffside, T. William Lambe, 
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GTS:May 81:684:16221 

Optimal Estimators for Soil Properties, Gregory B. 
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Mitchell, GT11:Nov 81:1461:16634 


Tensile Stress 

Laboratory Study of Hydraulic Fracturing, Gary W. 
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Unpaved Roads 
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Slip Circles and Critical Shear Planes, GT7:July 81: 
929:16369 


Sridharan, A. 

Embedded Foundations Under Vertical Vibration, 
with M. V. Nagendra and C. Chinnaswamy, TN: 
GT10:Oct 81:1429:16530 


See M. Krishna Murthy, TN: GT10:Oct 80:1143: 
15724 


Sridharan, Asuri 
See Mehter Mohamed Allam, GT4:Apr 81:421:16178 


Sridharan m116361, Asuri 


See Mysore V. Nagendra, TN: GT7:July 81:989: 
16360 


Stokoe, Kenneth H., II 


See Thomas M. Hearne, Jr., GT10:Oct 81:1327: 
16582 


Subrahmanyam, R. V. 
See Hin Fatt Cheong, TN: GT2:Feb 81:221:16009 


Sundaram, Panchanatham 
See Richard E. Goodman, GT5:May 80:485:15433 


Swiger, William F. 
Northfield Mountain 


Pumped Storage 
Performance, with Philip A. Wild and Thomas J. 
Lamb, GT6:June 80:673:15506 

Dsc: Claude A. Fetzer, GT7:July 81:1020:16356 

Clo: GT7:July 81:1021:16356 

TR:145:731 


Project: 


Symes, Matthew 
See Atonio Gens, GTS:May 81:702:16216 


Taft, Jay L. 
See Tavit O. Najarian, GT12:Dec 81:1141:16703 


Tan, David Y. 
Ground Control for 


Shallow Tunnels by Soil 


Grouting, with G. Wayne Clough, GT9:Sept 80: 
1037:15716 


TAN 


Dsc: Kurt K. Borchert and Herbert Klapperich, 
GT12:Dec 81:1745 sOASCE:16683 


Tang, Wilson H. 

Probabilistic Evaluation of Loads, GT3:Mar 81:287: 
16100 

See Achintya Haldar, GTS:May 81:577:16239 


Tokimatsu, Kohji 
See Yoshiaki Yoshimi, GT2:Feb 81:237:16003 


Townsend, Frank 
See Mosaid M. Al-Hussaini, GT4:Apr 81:481:16180 


Tsur-Lavie, 
See S. A. Denekamp, GT4:Apr 81:439:16173 


Vaid, Yoginder P. 

Dilation Angle and Liquefaction Potential, with Peter 
M. Byrne and John M.O. Hughes, TN: GT7:July 
81:1003:16360 

Static Shear and Liquefaction Potential, with W. D. 
Liam Finn, GT10:Oct 79:1233:14909 

Dsc: Yoshiaki Yoshimi and Kohji Tokimatsu, 
GT2:Feb 81:237:16003 
TR:145:555 


Veneziano, Daniele 
See Michael A. Fardis, GT10:Oct 81:1345:16590 
See Michael N. Fardis, GT10:Oct 81:1361:16604 


Verner, Edward A. 
See Leland M. Kraft, GTl:Jan 81:1:16000 


Vesic, Aleksandar 
See Nabil F. Ismael, GT12:Dec 81:1677:16718 


Vivatrat, Vitoon 
See Mohsen Baligh, GT4:Apr 80:447:15377 


Vrymoed, John 
Dynamic FEM Model of Oroville Dam, GT8:Aug 
81:1057:16464 


Wahls, Harvey 
Tolerable Settlement of Buildings, GT11:Nov 81: 
1489: 16628 


1770 


Wiesner, Terence 
Laboratory Tests on Model Piled Raft Foundations, 
with Peter T. Brown, GT7:July 80:767:15576 
Dsc: Joe O. Akinmusuru, GT8:Aug 81:1153: 
16408; Eldon Burley, Roger C. Harvey, Bhas 
Kar Nath and Lawrence A. Wood, GT8:Aug 81: 
1155:16408 
Clo: GT8:Aug 81:1157:16408 


Wild, Philip A. 
See William F. Swiger, GT6:June 80:673:15506 


Wiss, John F. 
Construction Vibrations: State-of-the-Art, GT2:Feb 
81:167:16030 


Wolfgang, James 
See John R. Hall, Jr., TN: GT10:Oct 81:1419:16530 


Wong, Kinfun 
See Tien H. Wu, GT12:Dec 81:1693:16764 


Wood, David 

Dsc, GTS:May 81:704:16216, of “Principal Stress 
Rotation: A Missing Parameter”, GT4:Apr 80:419: 
15375 


Wood, Lawrence 
See Eldon Burley, GT8:Aug 81:1155:16408 


Woods, Richard 
Seismic Techniques in the Laboratory, with Robert 
Henke, GT10:Oct 81:1309:16599 


Wu, Tien 
Probabilistic Soil Exploration: Case History, with 
Kinfun Wong, GT12:Dec 81:1693:16764 


Yoshimi, Yoshiaki 

Dsc, with Kohji Tokimatsu, GT2:Feb 81:237:16003, 
of “Static Shear and Liquefaction Potential”, GT10: 
Oct 79:1233:14909 


Znidarcic, Dobroslav 
Finite Strain Consolidation: Test Conditions, with 
Robert L. Schiffman, TN: GTS:May 81:684:16221 
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